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Spillways  at  lake  outlets  are  used  to  reduce  water- 
level  fluctuations  and  promote  year-round  lake  access. 
Such  hydraulic  alteration  may  cause  accelerated 
sedimentation  in  the  lake.   The  impact  of  a  1967  spillway 
was  quantified  in  eutrophic  Newnan's  Lake  (Florida)  using  a 
model  and  experimental  work. 

A  transect  of  sedimentary  profiles,  dated  with  210Pb 
and  137Cs  by  y-ray  spectroscopy,  showed  threefold  increases 
in  accumulation  rates  of  organic  matter,  total  Kjeldahl 
nitrogen  (TKN) ,  and  total  phosphorus  (TP)  1200  m  lakeward 
of  the  spillway  since  its  construction.   Concentrations  of 
TKN  and  TP  increased  3.5  and  2.4  times,  respectively,  in 
sediments  deposited  since  1967.   These  increases  were 
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progressively  less  at  stations  farther  from  the  spillway. 
Postspillway  accumulation  of  TP  was  focused  toward  the  dam, 
while  recent  TKN  deposition  was  similar  lakewide. 

A  90-day  spillway  removal  flushed  60  tons  (dry  weight) 
of  sediment  from  the  lake.   Discharge  concentrations  of 
particulate  organic  matter,  TKN  and  TP  were  highest  during 
the  first  month  of  the  drawdown  under  adequate  hydraulic 
head.   Storms  stirred  the  water  column  and  promoted 
flushing  of  resuspended  matter.   Field  and  laboratory  tests 
did  not  show  net  oxidative  removal  of  organic  matter  from 
exposed  lake  bottom.   Consolidated  sediments  remained 
moderately  firm  after  reflooding.   Lowering  the  water  depth 
from  62  to  30  cm  in  the  littoral  zone  removed  6.08 
g  m~2day-1  of  organic  littoral  sediment,  likely  due  to 
increased  wind-wave  action  on  this  substrate.   Littoral 
sediments  with  low  bulk  density  eroded  fastest  and  bulk 
density  of  remaining  substrate  increased  by  an  average  of 
250%.   Redistribution  of  this  material  to  deeper  portions 
of  the  lake  was  not  demonstrated. 

Results  from  experiments,  literature  and  theory  were 
synthesized  using  feedback  dynamics.   The  resulting 
simulation  model  was  most  sensitive  to  changes  in  sediment 
resuspension  and  actions  of  bacteria  and  benthic 
heterotrophs  on  phosphorus  dynamics  in  surface  sediments. 
Research  aimed  at  quantifying  these  processes  will  enhance 
confidence  in  the  model.   Separate  impacts  of  the  spillway 
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(decreased  flushing,  littoral  zone  vegetation,  and  sediment 
resuspension)  had  minor  implications  for  model  behavior. 
Combined,  however,  they  increased  seston,  water-column 
phosphorus,  and  sedimentation  by  19,  39,  and  27%, 
respectively. 

In  conclusion,  after  25  years,  intended  benefits  of 
the  spillway  (lake  access,  navigation)  have  been  replaced 
with  lake  management  "costs"  such  as  increased  turbidity, 
primary  production,  and  sedimentation. 
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CHAPTER  1 
INTRODUCTION 


Effects  of  spillways 

Spillways  have  been  built  at  the  outlets  of  many  large 
shallow  lakes  in  Florida  during  the  past  50  years  (Bishop 
1967;  Davis  1973).   The  majority  of  these  water  control 
structures  were  designed  to  prevent  low  lake-stages  during 
the  dry  winter  season  characteristic  of  Florida's  climate. 
Year-round  access  to  these  lakes  was  promoted,  and  their 
use  for  irrigation,  recreation,  and  sportfishing  maximized. 
The  structures  varied  from  simple  earthen  dams  to  stoplog 
weirs  and  concrete  spillways  with  or  without  removable 
boards . 

The  hydraulic  modification  produced  by  these  dams  may 
produce  short-term  benefits  for  lake  access  and  navigation, 
but  long-term  conseguences  for  lake  ecosystems  are  poorly 
understood.   Stabilized  lake  levels  may  accelerate  the 
accumulation  of  sediment  and  nutrient-rich  detritus  in  the 
lake.   This  can  occur  in  several  ways. 

First,  the  design  of  most  spillways  eliminates 
discharge  of  bottom-water.   This  layer  of  water  is  rich  in 
solids  due  to  sediment  resuspended  by  wind-wave  induced 


2 
currents  (Sheng  and  Lick  1979;  Paul  et  al.  1982) .   Outflow 
is  restricted  to  less  turbid  surface  water  flowing  over  the 
crest  of  the  spillway.   During  extended  periods  of  low 
rainfall  when  outflow  of  the  lake  is  confined  to  evapo- 
transpiration  and  losses  to  groundwater,  suspended  solids 
are  not  exported. 

Second,  decreased  water-level  fluctuations  reduce 
exposure  of  littoral  sediments  to  air,  thus  preventing 
seasonal  consolidation  of  this  flocculent  substrate  and 
reducing  oxidative  removal  of  organic  matter.   These 
effects  may  be  particularly  pronounced  in  shallow  lakes 
with  vast  littoral  areas.   Such  consolidation  has  been 
documented  when  the  water  level  in  a  lake  was  artificially 
lowered  for  several  months.   Organic  and  nutrient-rich 
sediments  from  Lake  Apopka  (Florida) ,  with  an  initial  water 
content  in  excess  of  80%,  consolidated  40-50%  after 
exposure  to  rain  and  sun  for  5  months  (Fox  et  al.  1977) . 
Significant  consolidation  of  the  top  5-10  cm  of  exposed 
littoral  substrate  was  noted  during  the  drawdown  of  Fox 
Lake  (McKinney  and  Coleman  1980) .   Oxidation  of  exposed 
lake  sediments  has  been  suggested  (Harris  and  Marshall 
1963;  Wegener  and  Williams  1974),  but  experimental  evidence 
is  controversial.   No  significant  oxidation  was  noted 
during  drying  of  Lake  Apopka  sediments  (Fox  et  al.  1977) , 
while  Holcomb  et  al.  (1975)  reported  a  considerable 
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reduction  in  organic  matter  content  of  exposed  sediments  in 
Lake  Tohopekaliga  (Florida) . 

Third,  stabilized  lake  levels  may  reduce  wetland 
habitat  around  the  lake  (Keddy  and  Reznicek  1985) .   Many 
wetland  plants  require  seasonal  fluctuations  in  water  level 
for  growth  and  reproduction  (Gosselink  and  Turner  1978; 
Bedinger  1979) .   Constant  water  levels  result  in  a  gradual 
decline  of  such  wetland  taxa  and  may  alter  the  wetland 's 
ability  to  filter  nutrients  and  sediments  carried  into  the 
lake  in  runoff  from  the  watershed. 

Increased  accumulation  of  flocculent  sediment  and 
nutrient-rich  detritus  may  significantly  affect  lake 
ecosystems.   Unconsolidated  sediments  are  easily 
resuspended  (Lam  and  Jacquet  1976;  Sheng  and  Lick  1979) 
resulting  in  increased  turbidity  in  the  water  column  and 
reduced  light  penetration.   This  may  change  plant 
communities  and  reduce  submerged  vegetation.   Periodic 
resuspension  of  sediment  and  associated  pore-fluid  also  may 
lead  to  increased  availability  of  nutrients  for  algal 
utilization  in  the  trophogenic  layer  (Holdren  and  Armstrong 
1980;  Pollman  1983).   Organic  and  inorganic  suspended 
matter  exert  considerable  biological  and  chemical  oxygen 
demand  (Hargrave  1969;  James  1974)  which  stresses 
heterotrophic  communities  in  the  lake  and  promotes  nutrient 
release  from  the  sediments  (Mortimer  1971;  Theis  and  McCabe 
1978;  Stauffer  1981).   Deposits  of  flocculent  sediment  in 


4 
the  littoral  zone  eliminate  firm  substrate  for  rooted 
macrophytes,  thereby  reducing  refuge  and  feeding  habitat 
for  fish.   Such  accumulations  may  also  reduce  the 
availability  of  preferred  nesting  habitat  for  centrarchids , 
such  as  the  largemouth  bass  (Micropterus  salmoides)  and 
black  crappie  (Pomoxis  nigromaculatus )  (Bruno  1984) ,  which 
rely  on  firm  substrate  for  deposition  of  eggs  (Eddy  and 
Underhill  1978) . 

Related  Studies 

In  spite  of  these  plausible  effects,  little  is  known 
about  the  long-term  response  of  a  lake  to  damming  of  its 
outlet.   Much  work  has  focused  on  short-term  responses 
following  drastic  water-level  manipulations  such  as  the 
filling  of  a  reservoir  by  damming  of  a  stream  or  the 
drawdown  of  an  existing  impoundment.   Many  of  these 
investigations  are  reviewed  in  Baxter  (1977)  and  Cooke  et 
al.  (1986) .   Those  applicable  to  this  work  are  referred  to 
in  Chapters  3  and  4. 

Long-term  investigations  have  concentrated  on  changes 
in  littoral  flora  following  water-level  stabilization  and 
the  possible  effects  of  such  a  reduction  in  high  water-low 
water  perturbation  on  the  ontogeny  of  aguatic  ecosystems. 
Symoens  et  al.  (1988)  reported  a  reduction  in  species 
diversity  in  the  littoral  zone  of  eutrophic  Lake  Virelles 
(Belgium)  due  to  water-level  stabilization  by  an  outflow 
control  structure.   Plants  characteristic  of  the 


periodically  exposed  shores  disappeared.   These  included 
sedges,  rushes,  and  members  of  the  buckwheat  family. 
Fluctuating  water-levels  increased  the  diversity  of 
vegetation  types  and  plant  species  along  the  shoreline  of 
the  Great  Lakes  (Keddy  and  Reznicek  1985) .   Low  water 
periods  allowed  many  mud-flat  annuals  and  emergent  marsh 
species  to  regenerate  from  buried  seeds.   Water-level 
stabilization  in  Lake  Tohopekaliga  (Florida)  eliminated  a 
considerable  portion  of  the  vegetated  flood  plain  (Holcomb 
and  Wegener  1971) . 

Effects  of  stabilizing  water-levels  on  the  ontogeny  of 
aquatic  ecosystems  have  been  described  for  wetlands 
(Gottgens  and  Montague  1988;  Gunderson  1989;  Walters  et  al. 
1992),  aquaculture  systems  (Odum  1971),  and  floodplain 
vegetation  (Goodrick  and  Milleson  1974) .   Freshwater 
marshes,  such  as  the  Florida  Everglades,  appear  to  be 
maintained  at  an  early  successional  stage  by  seasonal 
fluctuations  in  water  levels  and  other  pulsing  factors 
(e.g.  fire,  grazing) .   Aerobic  decomposition  of  accumulated 
organic  matter  may  be  promoted  when  sediments  are  exposed 
during  the  dry  season.   Nutrients,  released  upon 
reflooding,  support  a  wet-season  bloom  in  productivity.   A 
fluctuating  water  regime  may  be  a  dominant  force  in 
preventing  organic  matter  build-up  and  maintaining  an 
aquatic  ecosystem  in  a  state  of  pulsed  stability  (Odum 
1971)  . 
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Research  Methods:  Coupling  Limnology  with  Modeling 
Two  research  approaches  are  convenient  in  an 
investigation  of  long-term  changes  in  a  lake  following 
spillway  installation.   Both  pa leo limnology  and  ecosystem 
modeling  work  with  flexible  time-frames  and  allow 
investigation  of  long-term  and  slow-acting  processes. 
Paleolimnology  is  particularly  useful  when  a  detailed 
historical  data  base  is  absent. 

Paleolimnologists  assume  that  lake  sediments 
accumulate  in  an  orderly  fashion  through  time.   These 
sediments  contain  a  relatively  stable,  historical  record  of 
past  conditions  of  the  lake-watershed  ecosystem.   This 
record  may  be  in  the  form  of  physical,  chemical,  and/or 
biological  signals.   The  timing  of  these  signals  reguires 
dating  of  levels  in  sediment  cores  which  can  be 
accomplished,  for  instance,  by  measuring  the  content  of 
radioactive  carbon  (14C)  ,  lead-210  (210Pb) ,  or  cesium-137 
(137Cs) .   The  age  of  recent  sediments  (up  to  120  years)  can 
be  estimated  by  the  content  of  short-lived  isotopes  (e.g. 
210Pb)  (Eakins  and  Morrison  1978;  Krishnaswami  and  Lai 
1978)  .   In  this  work,  paleolimnological  technigues  are 
applied  to  guantify  the  impact  of  a  1967  spillway  on 
sediment  composition  and  net  rates  of  sediment  and  nutrient 
accumulation  in  Newnan's  Lake  (Chapter  2). 

Knowledge  of  past  lake  conditions  based  on 
interpretations  of  the  sedimentary  record  is  not  only 
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useful  to  evaluate  rates  of  change  in  the  lake  (such  as 
loss  of  volume)  ,  but  it  is  also  helpful  in  the  development 
of  lake  management  models,  because  it  can  assist  in 
reconstructing  the  reasons  for  past  changes  in  the  lake. 
Models  of  different  hypotheses  of  the  causes  of  past  lake 
responses  (such  as  an  increase  in  the  rate  of 
sedimentation)  to  perturbation  (such  as  the  construction  of 
a  dam)  may  be  used  to  predict  the  effect  of  future 
management  actions.   One  of  the  objectives  of  the 
simulation  model  presented  here  (Chapter  5)  is  to  evaluate 
the  long-term  effect  of  spillway  construction  in  the  outlet 
of  Newnan's  Lake  on  variables  such  as  sedimentation  rate, 
lake  depth,  and  concentrations  of  nutrients  and  suspended 
matter  in  the  water  column.   In  addition,  the  model  can 
provide  options  for  lake  management  by  testing  the  effects 
of  alternative  management  strategies  (such  as  a  modified 
dam  or  a  drawdown)  on  intended  uses  of  the  lake. 

The  scheduling  of  a  short-term,  experimental  drawdown 
in  Newnan's  Lake  in  the  spring  of  1989  provided  an 
opportunity  to  investigate  the  impact  of  this  management 
strategy.   This  investigation  focused  on  the  effect  of  the 
drawdown  on  littoral  and  profundal  lake  sediments  for  two 
reasons.   First,  the  accumulation  of  flocculent  and 
nutrient-rich  sediments  in  the  lake  may  have  significant 
implications  for  the  lake  ecosystem  (as  outlined  earlier) . 
This  material  may  be  removed  during  a  drawdown  by  flushing 
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of  suspended  matter  and/or  oxidation  of  exposed  littoral 
substrate.   Both  these  processes  were  quantified  in  this 
work  (Chapter  3).   Second,  by  focusing  on  the  sediments, 
the  effectiveness  of  the  drawdown  to  accomplish  its  main 
objective  may  be  determined.   The  main  objective  was  to 
improve  littoral  habitat  for  fish  growth  and  recruitment 
(Krummrich,  Florida  Game  and  Fresh  Water  Fish  Commission, 
pers.  comm.).   This  improvement  may  occur  when 
consolidation,  and/or  removal  of  flocculent  littoral 
substrate  promotes  the  establishment  of  littoral  plant 
communities  which  may,  in  turn,  provide  refuge  and  feeding 
habitat  for  fish.   Firm  littoral  substrate  also  provides 
preferred  nesting  habitat  for  species  of  the  sunfish  family 
(Centrarchidae)  (Bruno  1984),  economically  important 
sportfish  in  Newnan's  Lake.   This  study  measured  both 
consolidation  of  exposed  littoral  substrate,  and  erosion 
and  subsequent  "sloughing"  of  this  material  to  deeper  parts 
of  the  lake  (Chapter  4) . 

Finally,  studying  processes  in  the  lake  with  and 
without  the  spillway  complements  the  modeling  exercise.   It 
assists  in  determining  coefficient  values  for  modeling  the 
effect  of  the  spillway  on  the  lake,  and  it  allows  field 
testing  of  the  model  by  comparing  the  results  of  the 
drawdown  study  with  the  results  of  the  same  change  made  in 
the  model.   As  such,  the  limnological  work  assists  in  model 
development  and  analysis.   By  forcing  the  quantification  of 
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processes  in  the  lake,  the  modeling  exercise  may  reveal 
gaps  that  exist  in  our  understanding  of  this  ecosystem  and 
identify  necessary  field  or  laboratory  investigations. 

Study  Area 

Newnan's  Lake  is  located  8  km  east  of  Gainesville, 
Florida,  within  the  Orange  Creek  Basin  (Figure  1-1) .   The 
geomorphology  of  this  basin  was  first  described  by  Pirkle 
and  Brooks  (1959) .   It  is  dominated  by  the  Hawthorne 
formation,  a  marine  deposit  of  Miocene  age  consisting  of 
phosphatic  sands,  clays,  and  limestone.   It  is  relatively 
impermeable  compared  with  the  underlying  limestone  and  acts 
as  a  confining  layer.   Surface  flow,  subsurface  flow 
through  sands  and  clayey-sands,  and  direct  rainfall  are  the 
major  sources  of  water  for  the  lake  (Canf ield  1981) .   The 
drainage  area  north  of  the  lake  supplies  surface  water 
inflow  via  Hatchet  Creek,  Little  Hatchet  Creek,  and  several 
smaller  streams.   The  low  mineral  content  of  the  lake  water 
suggests  a  small  input  from  groundwater.   Loss  of  water 
from  the  lake  is  dominated  by  evapo-transpiration  in 
Florida's  warm  climate.   Of  the  130  cm  of  rainfall  received 
by  the  Orange  Creek  watershed  annually  (Adkins  1991) ,  only 
10%  leaves  as  surface  drainage  (Clark  et  al.  1964) . 

Land-use  in  the  basin  is  dominated  by  forest  (Table  1- 
1) ,  including  large  tracts  of  commercial  pine  (Pinus 
elliotti)  plantations.   Extensive  areas  of  undeveloped  pine 
flatwoods  also  occur,  with  riverine  hardwood-swamp  and 
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Figure  1-1.   Drainage  map  of  the  Orange  Creek  Basin. 
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cypress    (Taxodium  distichum)    along  the  water  courses.      A 
prominent  fringe  of  cypress  trees  borders  the   lake. 


Table   1-1.      Land-use  data  for 
the  Newnan's  Lake  watershed. 
Values   in  percent  of  total 
watershed  area. 

Urban  8 . 6 

Forest   1)  75.6 


Source:    Huber  et   al.    1982. 
)  Includes   areas   dominated  by 

cypress  trees. 


Agriculture  7 . 7 

Open  water  8.0 

Wetland  0.1 


Newnan's  Lake  has  a  surface  area  of  approximately  3000 
ha  and  a  maximum  and  mean  depth  of  3.6  m  and  1.5  m, 
respectively  (Nordlie  1976)  (Table  1-2) .   It  is  a  flat 
bottom  basin  with  a  few  depressions  where  the  water  depth 
exceeds  2.5  m  at  average  lake  stage  (Figure  1-2).   The  lake 
has  a  single  major  surface  water  outlet,  Prairie  Creek, 
which  drains  to  the  south  (Figure  1-1)  through  an  area  with 
extensive  marsh  and  bottomland  hardwood  communities.   Prior 
to  the  construction  of  Camp's  Canal  (in  the  early  1930s) 
Prairie  Creek  discharged  into  Payne's  Prairie  (Figure  1-1) . 
This  bypass  canal  was  constructed  around  the  eastern  margin 
of  the  prairie  to  control  flooding  and  manage  the  land  for 
cattle  grazing  (Gottgens  and  Montague  1988) .   The  canal 
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13 
directs  the  water  around  the  prairie  on  a  sluggish  course 
to  the  River  Styx  and  subsequently  to  Orange  Lake  (Figure 
1-1) .   The  prairie  has  been  managed  as  a  state  park  and 
marsh  preserve  since  1970  and  a  limited  amount  of  water 
from  Prairie  Creek  has  been  diverted  back  into  the 
preserve. 

Table  1-2.   Morphometry  of  Newnan's  Lake. 



Surface  Area  3  042  ha 

Maximum  Depth 

Mean  Depth 

Development  of  Shoreline 

Drainage  Basin  Area 

Lake  Volume 

Lake  Detention  Time 

The  outflow  from  Newnan's  Lake  through  Prairie  Creek 
has  been  regulated  with  a  spillway  since  1967.   This  50  m 
long  structure  consists  of  removable  flashboards  down  to 
the  channel  bed.   The  elevation  of  the  top  of  the  spillway 
boards  prevents  surface  discharge  when  the  water  level  in 
the  lake  drops  below  20.1  m  (MSL) .   During  this  study,  one 
of  the  top  boards  was  missing,  lowering  the  elevation  of  a 
1.2  m  section  of  the  top  of  the  spillway  boards  to  19.9  m 
(MSL) . 

A  comparison  of  pre-  and  postspillway  stage-duration 
curves  indicates  a  reduction  in  the  range  of  water-level 
fluctuations  of  approximately  32  cm  or  30%  (Skoglund  1990; 


14 
Adkins  1991) .   This  reduction  is  more  pronounced  if  severe 
prespillway  droughts  are  included,  during  which  tio  water 
level  records  were  kept  (1954-56  and  1962-64) .   In  this 
shallow  lake,  with  a  very  small  elevation  gradient  in  the 
littoral  zone,  such  reduced  water-level  fluctuations  may 
significantly  affect  seasonal  exposure  of  littoral 
sediments  and  the  area  of  vegetated  littoral  habitat.   The 
spillway  has  raised  the  average  stage  of  the  lake  by  13  cm, 
i.e.  it  has  increased  the  average  water  depth  in  the  lake 
by  9%. 

Newnan's  Lake  may  be  classified  as  an  eutrophic,  soft 
water  lake  (Canfield  1981) .   Water  quality  data  show  high 
nutrient  concentrations,  high  and  variable  chlorophyll  a 
values,  and  low  Secchi  disk  transparency  (Table  1-3) .   High 
color  (often  in  excess  of  150  Pt  units)  occurs  as  a  result 
of  input  of  water  rich  in  dissolved  humic  materials  from 
the  surrounding  flatwood  and  cypress  communities. 

Temperature  and  dissolved  oxygen  data  for  Newnan's 
Lake  were  also  typical  of  shallow,  productive  systems. 
Vertical  changes  in  water-column  temperature  rarely 
exceeded  3°C,  regardless  of  season  (Crisman  1986a). 
Consequently,  thermal  stratification  rarely,  if  ever, 
occurred.   Daytime  surface  water  during  summer  and  early 
fall  was  often  supersaturated,  while  the  water  column  in 
the  deeper  portions  of  the  lake  below  3  m  approached  anoxia 
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16 
as  a  result  of  intense  decomposition  in  surficial  sediments 
(Crisman  1986a) . 

Using  the  data  sources  from  Table  1-3 ,  the  ratio  of 
water-column  N:P  (in  grams)  varied  from  12.4  to  27.6.   This 
ratio  is  approximately  7:1  in  aquatic  plant  material 
(Vallentyne  1970) .   Therefore,  these  water  quality  data 
indicate  that  phosphorus  is  generally  less  available  than 
nitrogen  relative  to  their  content  in  aquatic  plant 
material,  and  that  primary  production  in  Newnan's  Lake  is 
likely  limited  by  phosphorus. 

The  lake  had  dense  growths  of  filamentous  algae 
dominated  by  blue-greens  (Cyanophyta)  such  as  Spirulina 
sp. ,  Anabaena  sp. ,  and  Aphanizomenon  sp. (Nordlie  1976; 
Crisman  1986a) .   At  times,  macrophytes  dominated  by  exotics 
such  as  water  hyacinth  (Eichhornia  crassipes)  and  hydrilla 
(Hydrilla  verticillata)  were  abundant.   Since  the  early 
1970s,  periodic  applications  of  herbicides  have  been  made 
to  restrict  the  spread  of  these  exotics.   The  lake  bottom 
was  covered  with  a  homogeneous  layer  of  highly  flocculent, 
organic  sediment  (Holly  1976) .   A  review  of  existing 
literature  and  an  inventory  of  the  existing  data  base  for 
Newnan's  Lake  is  given  in  Gottgens  and  Montague  (1987a), 
and  in  its  accompanying  categorized  bibliography  (Gottgens 
and  Montague  1987b) . 


CHAPTER  2 

QUANTITATIVE  IMPACTS  OF  LAKE  LEVEL  STABILIZATION 

ON  MATERIAL  TRANSFER  BETWEEN  WATER  AND  SEDIMENT 


In  Chapter  1,  it  was  hypothesized  that  spillways 
accelerate  the  accumulation  of  sediment  and  nutrient-rich 
detritus  in  a  lake  and  that  such  increased  trapping  of 
material  may  have  considerable  long-term  impact  on  the  lake 
ecosystem.   Little  is  known,  however,  about  increases  in 
the  rate  of  sedimentation  caused  by  damming  of  a  lake 
outlet.   Measurement  of  this  rate  is  essential  for 
determining  loss  in  lake  volume  and  guantifying  mineral 
cycling  in  lakes,  including  the  modeling  of  nutrient 
concentrations  in  water  (Dillon  and  Rigler  1974: 
Vollenweider  1976) .   It  is  also  needed  for  the  development 
of  a  sediment  or  pollutant  budget. 

Material  accumulation  rates  can  be  calculated  by 
determining  geochronology  in  the  sediment  profile  using 
radioactive  decay  of  fallout  ("unsupported")  lead-210 
(210Pb)  following  burial  in  sediments  (Appleby  and  Oldfield 
1978)  .   Matching  the  occurrence  of  cesium-137  (137Cs)  in 
the  profile  with  the  onset  of  widespread  atmospheric 
nuclear  testing  in  the  early  1950s  (Ritchie  et  al.  1973) 
provides  an  additional  age-marker.   These  paleolimnological 
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techniques  have  been  used  widely  to  detect  changes  in 
sediment  accumulation  rates  caused  by  urban  development  in 
a  watershed  (Smeltzer  and  Swain  1985) ,  major  storm  events 
(Robbins  et  al.  1978),  and  deforestation  (Oldfield  et  al. 
1980) .   They  are  particularly  helpful  in  gaining  insight 
into  past  conditions  in  lakes  when  a  long-term  data  base  is 
absent.   Here,  they  were  applied  to  quantify  the  impact  of 
a  1967  spillway  on  sediment  composition  and  net 
accumulation  rates  of  organic  matter,  nitrogen,  and 
phosphorus  in  Newnan's  Lake. 

Materials  and  Methods 
Collection  of  Sediment  Cores 

Core  locations  were  arranged  in  a  transect  of 
increasing  distance  lakeward  from  the  spillway  (Figure  2- 
1) .   Cores  8,  6,  and  1  were  taken,  respectively,  1200, 
3000,  and  4500  m  lakeward  of  the  structure  at  a  water 
depths  of  115,  150,  and  150  cm.   The  cores  were  collected 
from  the  lake  in  the  spring  of  1989  using  a  Livingstone 
piston  corer  (Livingstone  1955)  equipped  with  4 . 1  cm 
diameter  cellulose  butyrate  tubes.   The  piston  was 
positioned  at  the  lower  end  of  the  coring  tube  which  was 
then  lowered  to  the  level  at  which  sampling  was  to 
commence.   During  this  operation  the  piston  cable  was  payed 
out  freely,  but  once  the  sampling  level  was  reached  the 
cable  was  secured  to  the  boat  to  prevent  the  piston  from 
moving  any  farther  down.   The  tube  was  driven  into  the 
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-{w?}- wayside  Park  with  Boatramp 
H  F  f-Fishcamp 


Figure  2-1.   Insert  map  of  Florida  showing  location  of 
study  area  and  map  of  Newnan's  Lake  with  core  sites 
indicated  (filled  circles) . 
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sediments  by  pushing  on  the  extension  rod  while  the  piston 
prevented  the  sample  from  being  compressed  or  lost  upon 
retrieval.   After  the  apparatus  was  hoisted  to  the  surface, 
appropriate  caution  was  exercised  to  preserve  the  sediment- 
water  interface.   Cores  were  kept  refrigerated  until  they 
were  sectioned.   The  cores  were  sectioned  in  1  cm 
intervals,  which  were  stored  in  plastic  bags  at  4°C. 
Radio-isotope  Analyses 

Concentrations  of  210Pb  and  137Cs  were  measured  by 
direct  Y~assaY  using  a  coaxial  N-type,  intrinsic-germanium 
detector  (Princeton  Gamma  Tech)  .   The  counting  system  used 
for  spectral  analysis  (Figure  2-2)  is  located  in  the 
University  of  Florida's  Department  of  Environmental 
Engineering  Sciences'  Low  Background  Counting  Room.   This 
room  was  designed  to  reduce  the  level  of  background 
radiation  interference  during  sample  counting.   An  outer 
shield  (0.95  cm  steel),  main  shield  (10.1  cm  lead),  and  an 
inner  lining  (0.05  cm  cadmium  +  0.15  cm  copper)  were  used 
to  reduce  background  radiation  at  the  germanium  detector. 
The  detector  has  a  38%  efficiency  (at  1332.5  keV)  and 
operates  at  a  temperature  near  that  of  liquid  nitrogen 
(-196°C) ,  which  is  provided  by  a  high-vacuum  cryostat  dewar 
system  (capacity  30  1).   This  thin-window  (beryllium), 
standard  detector  has  a  relatively  low  efficiency  resulting 
in  the  need  for  long  count  times.   This  type  of  detector, 
however,  counts  over  a  large  range  of  Y~energies,  is 
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22 
suitable  for  many  different  sample  configurations,  and 
therefore  more  generally  useful.   It  also  detects 
efficiently  other  gammas  from  the  238u  decay  series  and  can 
thus  be  used  to  determine  supported  and  unsupported  levels 
of  210Pb  simultaneously  (Appleby  et  al.  1986;  Nagy  1988). 
The  electronics  include  a  preamplifier  (RG11B/C,  Princeton 
Gamma  Tech,  Inc.),  amplifier  (TC  242,  Tennelec) ,  bias 
supply  (5  kV,  TC  950,  Tennelec),  power  supply  (TC  909, 
Tennelec) ,  and  transformer  (Sola) .   A  Zenith  (Z159) 
computer  with  multi-channel  analyzer  ("Maestro"  ADCAM  100, 
EG  &  G  Ortec)  and  mathematical  spreadsheet  (Quattro, 
Borland  Inc.)  was  used  for  data  conversion  and  analysis.   A 
new  procedure  developed  for  Y~raY  spectroscopy  using  the 
Department's  recently  acquired  P-type,  intrinsic  germanium 
well  detector  is  given  in  Appendix  A. 

Samples  for  isotope  analysis  were  dried  at  95 °C  for  24 
hours,  pulverized  by  mortar  and  pestle,  weighed,  and  placed 
in  small  plastic  petri-dishes  (#1006  Falcon) .   Core 
sections  were  combined  (up  to  4  cm)  to  obtain  an  adequate 
sample  weight  (generally  more  than  1  gram) .   Petri  dishes 
were  sealed  with  plastic  cement  and  left  for  14  days  to 
equilibrate  radon  (222Rn)  with  radium  (226Ra) .   Counting 
times  varied  from  14  to  45  hours  depending  on  sample 
weight;  small  samples  needed  longer  counting  times  to 
minimize  uncertainty.   Blanks  were  counted  for  every  two 
samples  to  determine  background  radiation.   Standards  were 
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run  with  the  same  frequency  to  track  efficiency  (counts/y) 
and  to  calculate  a  226Ra  conversion  factor  (pCi/cps) . 
Sample  spectra  were  analyzed  for  activity  in  the  46.5  kev 
(210Pb)  and  662  kev  (137Cs)  peaks.   Activities  at  295  kev 
(214Pb)  ,  352  kev  (214Pb)  ,  and  609  kev  (214Bi)  representing 
uranium  series  peaks  were  used  to  compute  supported  levels 
of  210Pb. 

Calculation  of  210Pb  dates  followed  the  constant  rate 
of  supply  model  (Goldberg  1963;  Appleby  and  Oldfield  1978) 
which  is  able  to  quantify  changing  sediment  accumulation 
rates.   In  this  model,  the  cumulative  residual  unsupported 
210Pb,  At,  beneath  sediments  of  age  t  varies  according  to 


At  =  A0  e"kt  (2-1) 

where  A0  is  the  total  residual  unsupported  210Pb  (pCi/cm2) 
and  k  is  the  210Pb  radioactive  decay  constant.   At  and  A0 
are  calculated  by  numerical  integration  of  the  210Pb 
profile.   The  age  of  sediments  of  depth  x  is  then  given  by 


1     An 
f  |  in  _£  (2.2) 


The  sedimentation  rate  (r)  can  be  calculated  directly 
(Appleby  &  Oldfield  1978): 
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r=  2±  (2-3) 


where  C  is  the  concentration  of  unsupported  210Pb  (pCi/g) 
in  the  sediment  layer  of  interest.   This  model  appears 
applicable  to  Florida  lakes,  particularly  because  210Pb 
residuals  match  both  the  known  atmospheric  flux  of  this 
isotope  as  well  as  the  residuals  of  nearby  cores  (Binford 
and  Brenner  1986;  Gottgens  and  Crisman  1992)  . 
Bulk  Density,  Nutrient,  and  Carbon  Analyses 

Bulk  density  (mg/cm3)  and  organic  matter  content 
(mg/mg)  were  determined  on  1  cm3  subsamples  for  as  many  as 
45  intervals  per  core.   Bulk  density  was  measured  by  drying 
of  subsamples  at  95 °C  for  24  hours,  cooling  under 
dessication  and  weighing.   Organic  matter  content  was 
estimated  by  weight  loss  on  ignition  at  550 °C  for  1  hour 
followed  by  rehydration  with  distilled  water  and  re-drying 
at  95°C  for  24  hours  (method  209D,  A.P.H.A.  1985). 
Nitrogen  was  measured  as  total  Kjeldahl  nitrogen  (TKN) 
using  a  Technicon-II  semi-automated  manifold  after 
digestion  following  Bremner  and  Mulvaney  (1982)  ,  but 
modified  to  exclude  selenium  as  a  catalyst.   The  digestate 
was  also  used  for  total  phosphorus  (TP)  determinations. 
Liberated  orthophosphate  was  determined  with  the  ascorbic 
acid  method  (method  424F,  A.P.H.A.  1985)  using  a  Milton  Roy 
spectrophotometer  (Model  20)  with  a  1  cm  light  path. 
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Between  15  and  20  intervals  per  core  were  analyzed  for 
these  nutrients.   Total  carbon  in  the  sediments  was 
measured  by  combustion  of  the  sample  at  900°C  followed  by 
titration  of  the  evolved  C02  (Lee  and  Macalady  1989)  using 
a  model  5011  coulometer  (Coulometrics,  Inc.)-   Total 
organic  carbon  was  measured  using  a  combustion  temperature 
of  500°C. 

Results  and  Discussion 
Results  of  paleolimnological  analyses  may  be  presented 
in  units  of  concentration  or  as  rates  of  accumulation. 
Concentration,  expressed  as  a  relative  measure  of  sediment 
composition  (e.g.  mg/g  dry  weight) ,  is  the  conventional  way 
of  expressing  sediment  stratigraphy  (Shapiro  et  al.  1971; 
Pennington  1973;  Griffiths  and  Edmondson  1975;  Kramer  et 
al.  1991) .   Such  data,  however,  are  vulnerable  to 
variations  in  sedimentation  of  other  components  in  the 
profile.   For  example,  increased  deposition  of  allogenic 
inorganic  material  lowers  the  concentration  of  organic 
matter  in  the  sediment.   Effects  of  such  dilution  can  be 
eliminated  from  the  calculations  using  ratios  of  components 
in  the  sedimentary  matrix.   Accumulation  rates,  on  the 
other  hand,  are  normalized  to  time  thus  avoiding  the 
problem  of  co-variance  among  different  sedimentary 
components.   Detailed  sediment  dating  with  acceptable 
analytical  precision,  necessary  for  this  procedure,  can  now 
be  done  with  direct  210Pb  low-background  y-assay  (Appleby 
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et  al.  1986;  Gottgens  and  Crisman  1992) .   This  technique 
eliminates  uncertainties  associated  with  226Ra  supported 
activity  and  provides  additional  independent  data  markers 
in  the  profile  via  a  simultaneous  assay  for  other  y- 
emitting  isotopes,  including  137Cs.   Nonetheless, 
accumulation  rate  calculations  are  numerically  sensitive  to 
small  errors  in  210Pb  measurements.   The  uncertainty 
associated  with  this  may  be  reduced  by  averaging 
accumulation  rates  over  longer  periods  of  time.   In  this 
study,  both  concentration  and  accumulation  rate  data  are 
given  to  provide  complementary  information  about  the  recent 
history  of  the  lake. 

Cores  consisted  mostly  of  homogeneous  brown  material 
with  organic  matter  content  generally  greater  than  50%  and 
bulk  density  less  than  70  mg/cm3  throughout  the  top  30-40 
cm  (Gottgens  and  Crisman  1992) .   Abrupt  increases  in  the 
organic  matter,  TKN,  and  TP  profiles  of  core  8  at  33  cm 
depth  coincided  with  a  210Pb  determined  sediment  age  of 
approximately  1967,  i.e.  the  time  that  the  lake  outflow  was 
dammed  (Figure  2-3a) .   Concentrations  of  TKN  and  TP 
increased  3.5  and  2.4  times,  respectively,  in  sediments 
deposited  since  1967.   This  signal  was  progressively  less 
at  stations  farther  from  the  spillway  (Figures  2-3b  and  c)  . 
Increased  nutrient  concentrations  at  a  210Pb  date  of 
approximately  1952-56  (Figure  2-3a,  b,  and  c)  matched  the 
driest  year  (1954)  in  north  Florida's  recent  history.   This 
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Figure  2-3.   Core  sites  8  (a),  6  (b)  ,  and  1  (c) ,  Newnan's 
Lake:  Profiles  for  organic  matter,  total  phosphorus,  and 

total  Kjeldahl  nitrogen.   Spillway  installation  ( )  and 

first  occurrence  of  *37Cs  marker  (***)  indicated. 
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may  have  resulted  from  reflooding  of  exposed  littoral 
sediments  and  subsequent  release  of  nutrients  to  the  water 
column. 

The  first  occurrence  of  detectable  137Cs,  marking  the 
onset  of  widespread  atmospheric  nuclear  testing  in  the 
early  1950s,  matched  well  with  the  determined  210Pb 
chronology  (Figures  2-3a,  b,  and  c) .   This  agreement  is 
remarkable  considering  the  shallow  character  of  the  lake 
and  the  flocculence  of  the  bottom  substrate  which  makes 
these  upper  sediments  vulnerable  to  physical  disturbance. 
Eutrophic  systems  such  as  Newnan's  Lake,  however, 
accumulate  sediment  rapidly,  so  that  these  periodic 
physical  disturbances  likely  affect  short  time-intervals. 
The  bulk  density  profile  for  core  8  (Figure  2-4) 
demonstrated  both  the  fluid  nature  of  the  deposits  and  the 
sudden,  fivefold  decrease  in  bulk  density  after  spillway 
installation.   Bulk  density  in  the  top  33  cm  of  sediment 
(e.g.  24  years)  was  less  than  60  mg/cm3,  and  average  net 
accumulation  is  1.3  cm  per  year.   In  most  sediment 
profiles,  bulk  density  increases  with  depth  due  to 
compaction.   Yet,  the  profile  illustrates  a  significant 
transition  at  the  time  the  lake  outlet  was  dammed. 

Unconsolidated  sediments  may  be  resuspended  easily  by 
wind-induced  currents  in  shallow  systems  such  as  Newnan's 
Lake.   The  resulting  organic  turbidity  increases  sediment 
oxygen  demand  (Bowman  and  Delf ino  1980) ,  particularly  in 
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Figure  2-4.   Bulk  density  profile  for  core  site  8,  Newnan's 
Lake. 
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subtropical  lakes  where  high  temperatures  stimulate 
bacterial  respiration  in  organic  substrate  (McDonnell  and 
Hall  1969) .   In  addition,  resuspended  sediments  may 
contribute  significant  guantities  of  nutrients  to  the 
overlying  water  column  through  desorptive  processes  or 
entrained  pore  fluid  (Bengtsson  1975;  Cooke  et  al.  1977) . 
Pollman  (1983)  reported  a  more  than  twofold  increase  in 
mid-lake  orthophosphorus  concentrations  under  a  moderate 
wind  event  (approximately  8.9  m/s)  in  a  central  Florida 
lake  with  a  mean  depth  and  sediment  composition  similar  to 
Newnan's  Lake.   Pronounced  increases  in  the  nitrogen  and 
phosphorus  content  of  surficial  sediments  deposited  since 
spillway  installation  may  enhance  internal  nutrient  cycling 
and,  therefore,  change  plant  communities  and  lake 
productivity.   For  example,  using  data  from  Pollman 's 
(1983)  sediment-dispersion  nutrient-release  model, 
resuspension  of  the  upper  1  cm  of  flocculent  sediment  at 
station  8  would  cause  a  70%  increase  in  the  phosphorus 
concentration  of  the  overlying  water  through  entrainment  of 
interstitial  water  and  desorption  from  suspended  sediments. 
Carbon: nitrogen: phosphorus  ratios  (C:N:P)  demonstrated 
both  temporal  and  spatial  changes  in  the  composition  of  the 
deposits  (Table  2-1) .   Measurement  of  total  carbon  (TC)  and 
total  organic  carbon  (TOC)  demonstrated  that  almost  all 
sedimentary  C  is  organic  (TOC  =  0.996  x  TC;  standard  error 
=  0.003).   Postspillway  deposits  were  low  in  P  compared 
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with  material  deposited  prior  to  1967.   This  reduction  in  P 
content  was  most  pronounced  relative  to  N  content  at 
stations  6  and  1  (Table  2-1) .   While  many  factors  could  be 
considered  in  interpreting  this  difference,  including 
changes  in  fractionation  of  sedimentary  P  and  variations  in 
Fe  content  of  the  sediments  (Engstrom  and  Wright  1984) , 
such  information  was  not  collected  for  this  work.   The 
reduced  P  content  of  recent  sediments  may  be  a  reflection 
of  a  change  toward  P-limitation  of  primary  production  in 
the  lake. 


Table  2-1.   Total  carbon  (C)  to  total 
Kjeldahl  nitrogen  (N)  to  total  phosphorus 
(P)  ratios  (weight /weight)  for  cores  8, 
6,  and  1  pre-  and  postspillway 
installation,  Newnan's  Lake  (FL) . 


Core 

8 

6 

1 

Post  C:P 

82 

91 

91 

Pre  C : P 

40 

44 

43 

Post  C:N 

13 

10 

8 

Pre  C:N 

7 

6 

6 

Post  N:P 

6 

10 

11 

Pre  N : P 

6 

7 

7 

Spatially  within  the  lake,  the  amount  of  P  in  the 
sediments  relative  to  C  and  N  was  higher  in  "downstream" 
areas  (e.g.  site  8) .   This  may  result  from  rapid 


32 
entrainment  of  P  and  reduced  N:P  ratios  during  windy 
conditions  in  shallow  systems  (Hamilton  and  Mitchell  1988) 
followed  by  downstream  transport  and  density-dependent 
settling  of  the  fine  particulates  with  which  P  is 
associated  (Engstrom  and  Wright  1984) .   High  C:N  ratios  in 
core  8  deposits  may  reflect  a  somewhat  higher  degree  of 
decomposition  of  bottom  material  in  the  lower  portion  of 
the  lake.   The  low  C:N  ratios  lakewide  (between  6  and  13) 
imply  an  autochthonous  origin  of  the  organic  matter  in  the 
sediment  matrix  (cf .  Wetzel  1983) .   The  contribution  of 
non-Kjeldahl  nitrogen  to  this  ratio  is  small  in  anaerobic 
sediments  (Reddy  and  Patrick  1984) . 

Net  accumulation  rates  of  sediment,  organic  matter, 
and  nutrients  have  increased  greatly  close  to  the  spillway 
since  its  construction  (Figure  2-5) .   This  increase  was 
progressively  less  at  sites  farther  from  the  spillway, 
where  accumulation  rates  corresponded  well  with  mid- lake 
values  reported  for  lakes  of  similar  trophic  state  (Deevey 
et  al.  1986:  Binford  and  Brenner  1986).   The  highest  rates 
were  found  consistently  at  the  station  most  distant  from 
significant  material  inflow  in  the  northern  portion  of  the 
lake.   A  45%  higher  sedimentation  rate  at  station  8  amounts 
to  an  additional  11  cm  of  material  accumulation  and  a  9% 
reduction  in  water  depth  in  this  shallow  system  during  the 
24  years  since  dam  construction. 
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A  pronounced  increase  in  the  accumulation  rate  of  P  at 
station  8  was  observed  (3-3.5  times  the  measured  rate  at 
stations  6  and  1) ,  while  recent  N  deposition  is  similar 
among  the  three  stations  (Figure  2-5) .   This  "focusing"  of 
P  in  sediments  in  the  lower  portion  of  the  lake  parallels 
the  decreased  C:P  and  N:P  ratios  at  this  downstream 
station.   The  higher  rate  of  entrainment  of  P  (compared  to 
N)  by  wind-induced  wave  action  (Maceina  and  Soballe  1990) 
likely  contributes  to  P  focusing  in  this  shallow,  exposed 
lake.   Its  short  water  residence  time  and  obstructed 
outflow  promote  horizontal  transport  and  subsequent 
accumulation  close  to  the  dam. 

Other  factors  may  have  contributed  to  increased 
sedimentation,  but  their  effect  would  have  resulted  in 
changed  sedimentation  rates  lake-wide,  rather  than 
concentrated  toward  the  spillway.   These  factors  include 
herbicide  treatment  of  aquatic  macrophytes  and  added 
nutrient  runoff  from  the  watershed.   Between  1972  and  1977 
a  total  of  130  ha  (i.e.  4%)  of  lake  surface  area  was 
treated  chemically  to  reduce  the  abundance  of  Eichhornia 
crassipes  (water  hyacinth)  and  Hvdrilla  verticillata 
(hydrilla) .   In  November  of  1988  an  additional  110  ha 
received  chemical  treatment  (Hinkle,  Florida  Department  of 
Natural  Resources,  pers.  comm.).   Analysis  of  aerial 
photographs  (U.S.  Geological  Survey  1966,  1988)  did  not 
show  significant  land-use  changes  in  the  basin  during  the 
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last  22  years.   The  number  of  residences  less  than  1  km 
from  the  Newnan's  Lake  shoreline  increased  from  16.0  to 
19.2  units/km2  between  1966  and  1988. 

Summary  and  Conclusions 
The  Newnan's  Lake  spillway,  intended  to  stabilize  lake 
stage,  changed  material  transfer  between  water  and 
sediment.   This  change,  operating  over  a  period  of  decades, 
can  be  quantified  using  paleolimnological  techniques.   A 
transect  of  210Pb  dated  profiles  perpendicular  to  the  dam 
at  the  lake  outlet  demonstrated  a  distinct  impact  from  its 
installation  in  1967  on  net  rate  of  material  accumulation 
and  sediment  composition.   While  net  accumulation  rate  of 
bulk  sediment  was  only  45%  higher,  resulting  in  a  9% 
reduction  in  water  depth,  rates  for  organic  matter  and 
nutrients  tripled  compared  to  prespillway  conditions.   P 
accumulation  was  particularly  enhanced  close  to  the 
spillway  compared  with  other  stations.   If  allowed  to 
persist,  increased  material  transfer  between  water  and 
sediments  may  produce  loss  of  lake  volume  and  changes  in 
benthic  habitat.   The  threefold  increase  in  the  trapping  of 
nutrients  may  change  plant  communities  and  lake 
productivity . 


CHAPTER  3 

REMOVAL  OF  PARTICULATE  ORGANIC  MATTER/ NUTRIENTS  AND 

OXIDATION/ CONSOLIDATION  OF  EXPOSED  LITTORAL  SUBSTRATE 

DURING  A  SHORT-TERM  DRAWDOWN 


A  three-month  gravity  drawdown  of  the  lake  was 
initiated  on  April  24,  1989  by  removing  all  spillway 
flashboards.   The  objectives  of  this  drawdown  were  to 
improve  fish  growth  and  recruitment  and  to  flush  nutrient- 
rich  detritus  from  the  lake.   Accumulation  of  flocculent 
and  nutrient-rich  sediments  in  the  lake  (Chapter  2)  and  the 
presence  of  this  flocculent  substrate  in  the  littoral  zone 
may  reduce  the  availability  of  preferred  nesting  habitat 
for  centrarchids ,  such  as  the  largemouth  bass  (Micropterus 
salmoides)  and  black  crappie  (Pomoxis  nigromaculatus )     (Eddy 
and  Underhill  1978;  Bruno  1984).   Exposing  the  littoral 
zone  lake-bottom  during  drawdown  may  allow  oxidation  and 
consolidation  of  this  substrate,  which  may  improve  this 
habitat  for  future  fish  spawning.   Seed  germination  and 
firm  substrate  promote  the  establishment  of  littoral  plant 
communities,  which  may  provide  refuge  and  feeding  habitat 
and  promote  fish  recruitment. 

Research  on  other  Florida  lakes  has  guantified  effects 
of  drawdowns  on  fish  populations  (Wegener  and  Williams 
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1974)  ,  aquatic  invertebrates  (Wegener  et  al.  1974) ,  and 
littoral  plant  communities  (Holcomb  and  Wegener  1971; 
Goodrick  and  Milleson  1974;  Hestand  and  Carter  1975). 
However,  no  information  is  available  on  the  impact  of 
drawdown  practices  on  the  flushing  of  organic  matter  and 
nutrients  from  Florida  lakes.   The  objectives  of  this  part 
of  the  study  were:   (1)  To  measure  the  flushing  of  organic 
matter,  nitrogen  and  phosphorus  in  surface  discharge  from 
Newnan's  Lake  prior  to  and  during  the  gravity  drawdown  and 
(2)  to  determine  the  effect  of  exposing  littoral  zone 
sediments  to  air  on  oxidative  removal  of  organic  matter  and 
long-term  consolidation  of  the  substrate. 

Materials  and  Methods 
Discharge.  Rainfall,  and  Lake  Stage 

Discharge  through  Prairie  Creek,  rainfall,  and  lake 
water-level  were  measured  starting  February  20,  1989,  9 
weeks  prior  to  drawdown,  and  continued  throughout  the 
drawdown  period  until  July  31,  1989,  2  weeks  after  the 
spillway  boards  were  re-installed. 

Surface  flow  from  the  lake  was  measured  weekly  using 
an  Ott  flow-meter  (Type  C2-10150)  approximately  30  m 
downstream  from  the  dam  in  Prairie  Creek  (Figure  3-1) .   The 
frequency  of  sampling  was  increased  during  storm  events  to 
include  a  minimum  of  1  prestorm  flow,  1  peak  flow,  and  1 
poststorm  flow  measurement.   The  creek  width  was  divided  in 
1  m  subsections,  and  flow  was  measured  in  each  section  at 
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Figure  3-1.  Map  of  Florida  showing  study  area  location  and 
map  of  Newnan's  Lake.   Open  circles  indicate  the  locations 
of  field  enclosures,  and  closed  circles  indicate  the 
sediment  collection  sites  for  the  oxidation/consolidation 
laboratory  experiments. 
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0.6  depth  (when  channel  depth  <  30  cm)  or  at  0.2  and  0.8 
depth  and  averaged  (when  channel  depth  >  30  cm) .   Discharge 
was  computed  by  multiplying  flow  in  each  subsection  by  its 
cross-sectional  area.   Total  discharge  equals  the  sum  of 
all  subsection  discharges.   If  the  measured  discharge 
differed  from  the  discharge  computed  using  USGS 
stage-discharge  ratings  by  more  than  10%,  then  discharge 
was  remeasured  and  averaged  with  the  first  measurement.   A 
calibration  flume  (Department  of  Civil  Engineering, 
University  of  Florida)  was  used  to  determine  the  accuracy 
of  the  Ott  flow-meter. 

The  Gainesville  Flight  Service  Station  (National 
Oceanic  and  Atmospheric  Administration)  kept  daily  rainfall 
records  during  the  study  period  using  a  standard  (20.32  cm 
diameter)  U.S.  Weather  Service  raingauge,  located  near  the 
center  of  the  watershed.   Lake  levels  were  recorded  daily 
by  the  U.S.  Geological  Survey  (station  #  02240900)  located 
on  the  west  side  of  the  lake. 
Water  Quality  Analysis 

Replicate  water  samples  were  collected  concurrently 
with  discharge  measurements  20  m  upstream  from  the  spillway 
in  the  center  of  the  creek.   Water  samples  were  collected 
in  acid  washed  Nalgene  containers,  immediately  stored  on 
ice,  and  frozen  within  1  hour  of  collection.   Samples  were 
analyzed  for  total  suspended  solids  (TSS) ,  particulate 
organic  matter  (POM) ,  total  Kjeldahl  nitrogen  (TKN) ,  and 
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total  phosphorus  (TP) .   The  analyses  were  completed  within 
90  days  of  collection. 

Analysis  for  TSS  was  according  to  Standard  Methods 
(method  209C,  A.P.H.A.  1985)  using  pre-weighed  and 
pre-muffled  glass  fiber  filters  (Whatman,  934-AH,  4.25  cm 
diameter)  with  an  effective  particle  retention  of  1.5  /xm. 
Replicate  analyses  were  done  on  0.5  liter  samples.   POM  was 
measured  as  weight  loss  on  ignition  (at  550°C  for  1  hour, 
followed  by  rehydration  with  distilled  water  and  drying  at 
95°C  for  24  hours)  (method  209D,  A.P.H.A.  1985)  of  oven-dry 
samples  (95°C  for  24  hours) .   TKN  analyses  were  according 
to  E.P.A.  method  351.2  (E.P.A.  1979)  using  a  40-cell  block 
digestor  and  a  Technicon  II  semi-automated  manifold.   The 
digested  sample  was  also  used  for  TP  determination.   The 
liberated  orthophosphate  in  the  digested  samples  was 
determined  using  the  ascorbic  acid  method  (method  424F, 
A.P.H.A.  1985) .   Absorbance  of  the  samples  was  read  at  880 
nm  on  a  Bausch  and  Lomb  spectrophotometer  (Model  21)  with  a 
light  path  of  2.5  cm.   Both  nitrogen  and  phosphorus 
analyses  were  done  on  replicate  samples  and  averaged. 
Sediment  Oxidation 

Oxidation  rates  of  organic  matter  in  exposed  littoral 
zone  sediment  were  measured  in  the  laboratory  and  the 
field.   Littoral  surface  sediment  was  collected  from  3 
stations  in  the  lake  (Figure  3-1)  and  homogenized  in  a 
blender  in  the  laboratory.   Subsamples  (n=25)  of  100  ml 
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each  were  incubated  in  crucibles  at  room  temperature  for  78 
days  under  3  treatments.   One  set  was  allowed  to  dry 
completely  and  remained  dry  during  the  test  period;  one  set 
was  kept  wet  without  standing  water  by  adding  up  to  8  ml  of 
distilled  water  per  week,  and  a  third  set  was  kept 
inundated  (with  2-3  cm  of  distilled  water) .   The  duration 
of  this  experiment  (78  days)  was  chosen  to  coincide  with 
the  expected  length  of  time  that  littoral  substrate  would 
be  exposed  in  the  field.   Crucibles  in  the  laboratory  were 
rotated  every  48  hours  to  promote  egual  exposure  to  light 
conditions  (generally  between  20-30  Lux,  measured  with  a 
Li-Cor  photometer,  Model  LI  188)  .   Organic  matter  content 
before  and  after  incubation  was  determined  using  weight 
loss  on  ignition  (at  550 °C  for  1  hour  followed  by 
rehydration  with  distilled  water  and  drying  at  95 °C  for  24 
hours)  of  oven-dry  samples  (95°C  for  24  hours) . 

PVC  enclosures  (surface  area  =  180  cm2)  were  installed 
in  pairs  at  3  locations  in  exposed  littoral  zone  sediments 
(Figure  3-1) .   The  enclosures  were  pushed  through  the  top 
layer  of  wet  organic  substrate  to  a  point  5-10  cm  into  the 
underlying  sand.   Enclosed  sediment  was  shielded  from 
precipitation  by  an  elevated  plexiglass  roof.   This  design 
prevented  removal  of  enclosed  substrate  by  water  and/or 
wind  erosion.   Hardware  cloth  (1.25  cm  mesh  size)  was  put 
on  top  of  the  enclosure  to  deter  animal  disturbance. 
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Organic  matter  content  at  the  beginning  of  exposure 
was  quantified  by  average  weight  loss  on  ignition  of  3 
"core"  samples  adjacent  to  the  enclosure.   Each  sample 
covered  an  area  equivalent  to  the  enclosure  and  was 
analyzed  in  its  entirety  down  to  sand.   After  28  days  of 
incubation,  the  organic  matter  content  of  the  enclosed 
substrate  down  to  sand  was  determined.   Net  oxidative 
removal  (g  m~2day_1)  was  computed  as  the  difference  between 
pre-  and  postincubation  measurements. 
Sediment  Compaction 

A  small-scale  sediment  compaction  study  was  also 
performed  to  evaluate  changes  in  substrate  bulk  density, 
water  content  and  percent  organic  matter  at  various  times 
after  exposure  and  after  reflooding.   Four  vessels  (glass, 
1  liter)  containing  500  ml  of  homogenized  littoral 
substrate  were  allowed  to  dry  in  the  laboratory  for  51 
days,  then  inundated  (to  a  depth  of  10  cm)  for  149  days, 
simulating  the  water  regime  in  the  field.   Incubation  of 
the  substrate  occurred  at  room  temperature,  indoor  light 
conditions,  and  approximately  40%  relative  humidity.   The 
locations  from  which  sediment  samples  were  collected  are 
shown  in  Figure  3-1.   Only  the  top  5  cm  of  substrate  was 
used.   Ten  water  content  (percent)  and  organic  matter 
content  (percent)  determinations  were  made  at  5  times 
during  exposure,  at  the  time  of  reflooding,  and  4  times 
during  reflooding. 
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Results  and  Discussion 
Discharge.  Rainfall,  and  Lake  Stage 

During  the  3 . 5  months  prior  to  drawdown  lake  stage 
dropped  from  20.33  m  (MSL)  to  20.07  m  (MSL)  due  to  lack  of 
rain.   The  boards  in  the  spillway  were  removed  on  April  24 
(1989)  after  which  the  water  level  in  the  lake  dropped  from 
20.07  m  (MSL)  to  19.80  m  (MSL)  in  1  month  (Figure  3-2). 
Lake  stage  remained  between  19.70  and  19.80  m  MSL  for  8 
weeks.   This  stage  was  the  lowest  recorded  during  the  last 
25  years,  with  the  exception  of  a  brief  period  during  the 
drought  of  1981  (when  the  lake  reached  19.65  m  MSL).   The 
maximum  drop  in  water  level  during  drawdown  was  36  cm. 
Severe  lack  of  rain  during  the  months  prior  to  drawdown 
(Table  3-1)  resulted  in  lake  water-levels  approximately  41 
cm  below  average  for  early  spring.   This  reduced  the 
amplitude  of  the  drop  in  lake  water-level.   Small  elevation 
gradients  in  the  basin  and  the  gradual  build-up  of 
obstructions  in  the  lake  outflow  upstream  and  downstream 
from  the  spillway  since  its  construction  (Crider  1972; 
personal  observation)  likely  prevented  a  more  dramatic 
drawdown  of  water  level.   Those  obstructions  to  the  flow  in 
Prairie  Creek  may  have  resulted  from  the  impact  of  the 
spillway  on  stage  level  and  discharge  volumes  in  the  creek. 
Low-flow  and  low-stage  conditions  downstream  from  the 
spillway  are  much  reduced  compared  to  the  natural 
hydroperiod  of  the  creek  prior  to  dam  construction  in  1967 
(Gottgens  1987) . 
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Figure  3-2.  Water  level  record  for  Newnan's  Lake  for  the 
period  of  January-September  1989  [Source:  U.S.  Geological 
Survey  station  #  02240900]. 
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Table  3-1.   Average  rainfall  (1897-1987)  and  rainfall 
during  first  7  months  of  1989.   Records  given  in  cm. 

Jan  Feb 

Mar 

Apr  May 

Jun   Jul   Total 

Average 

1989 

Deficit 

7.7   9.0 
2.6   3.0 
5.1   6.0 

9.1 
5.5 
3.6 

7.6   8.9 
7.4   5.2 
0.2   3.7 

16.7  18.2  77.2 
24.4  10.0  58.1 
-7.7   8.2  19.1 

[Source:  National  Oceanic  and  Atmospheric  Administration,  Gainesville 
Flight  Service  Station,  Florida] 


Assuming  a  reduction  of  the  mean  water  depth  from  1.5 
m  to  1.2  m  during  the  drawdown,  it  is  estimated  that  20%  of 
water  volume  of  the  lake  was  removed  during  the  drawdown 
period.   Assuming  unaltered  base-discharge  rates  from  the 
lake  through  Prairie  Creek  during  the  drawdown  period, 
integration  of  the  discharge  curve  (Figure  3-3)  and  then 
subtracting  base-discharge  yields  an  estimate  of  1.75  x  io6 
m3  of  water  flushed  from  the  lake  by  removal  of  the 
spillway.   Base-discharge  is  defined  as  the  discharge 
flowing  through  the  channel  with  the  spillway  in  place,  and 
is  estimated  by  measuring  average  discharge  prior  to  dam 
removal  and  after  re-installation  of  the  flashboards.   In 
both  cases  an  averaging  period  of  2  weeks  was  used.   Below 
average  rainfall  during  the  drawdown  period  (Table  3-1)  and 
high  evapo-transpiration  contributed  considerably  to  the 
decrease  in  lake  stage.   Discharge  through  Prairie  Creek 
dropped  to  approximately  twice  the  assumed  base-discharge 
within  6  weeks  after  spillway  removal  and  remained  at  that 
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Figure  3-3.   Discharge  (m3/sec)  through  Prairie  Creek  § 
Florida  S.R.  20  for  the  period  of  February-July  1989. 
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rate  until  closing  of  the  weir  (Figure  3-3) .   In  spite  of 
low  total  rainfall,  5  storm  events  of  low  to  moderate 
intensity  were  included  in  the  sampling  period (Table  3-2). 


Table  3-2.   Precipitation  and  wind  conditions  during 
sampled  storm  events,  Newnan's  Lake. 


Dates  (1989) 


Amount     Wind  conditions 
(cm) (km/hr) 


Pre- 

drawdown 

During 
drawdown 


March  22-23  1.09 

April  14-15  2.03 

May  29  2.26 

June  18-19  4.52 

July  16-17  1.42 


0-18  from  N 
0-12  from  N  or  W 

0-40  misc.  dir. 
0-16  misc.  dir. 
0-16  misc.  dir. 


Measured  at  spillway  with  standard  2.54  cm  glass  raingage  and  Dwyer 
handheld  windmeter 


Removal  of  Particulate  Matter  and  Nutrients 

Increased  flow  through  Prairie  Creek  after  spillway 
removal  increased  discharges  (kg/day)  of  particulate 
organic  matter  (POM) ,  total  Kjeldahl  nitrogen  (TKN) ,  and 
total  phosphorus  (TP)  (Figure  3-4a,b,  and  c) .   POM  was  a 
consistent  fraction  of  total  suspended  solids  during  the 
sampling  period  (POM  =  0.76  x  TSS;  N=64;  R2=0.96). 
Assuming  base-discharge  rates  (kg/day)  through  Prairie 
Creek  during  the  drawdown  period,  integration  of  total 
discharge  rates  and  subtraction  of  base-discharge  rates 
yields  an  estimate  of  the  amount  (kg)  of  POM,  TKN,  and  TP 
flushed  from  the  lake  owing  to  drawdown  (Table  3-3) . 
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Figure  3-4.   Characteristics  of  the  discharge  from  Newnan's 
Lake  through  Prairie  Creek.   Opening  and  closing  of 
spillway  indicated.   Discharge  rate  (kg/day)  of  POM,  TKN, 
and  TP  (Figure  3-4a,  b,  and  c) .   Concentration  (mg/1)  of 
POM,  TKN,  and  TP  in  discharge  (Figure  3-4d,  e,  and  f ) . 
Arrows  indicate  sampled  storm  events. 
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Table  3-3.   Amounts  of  total  suspended  solids,  particulate 
organic  matter,  total  Kjeldahl  nitrogen,  and  total 
phosphorus  flushed  in  excess  of  base-discharge  from 
Newnan's  Lake  during  partial  drawdown. 

kg  dry  weight   mg/m2  lake  area 
removed  by       removed  by 
drawdown         drawdown 

Total  suspended  solids  59,247  2,043 
Particulate  organic  matter  46,537  1,605 
Total  Kjeldahl  nitrogen  8,840  305 
Total  phosphorus 290 10 


Computed  on  a  m2-basis,  removal  during  drawdown  is 
small  compared  to  the  likely  stores  of  flocculent  sediment 
in  the  lake,  with  an  average  estimated  thickness  of  70  cm 
(Skoglund  1990) .   Two  factors  may  account  for  the 
relatively  low  removal  rates  of  particulate  material  from 
the  basin.   First,  the  lower  incidence  of  storm  events 
during  the  period  of  drawdown  may  have  precluded 
considerable  resuspension  of  flocculent  bottom  material 
(Pollman  1983).   Wind-induced  wave  action  will  stir  the 
water  column  and  may  promote  flushing  of  resuspended  matter 
(Maceina  and  Soballe  1990) .   Second,  the  low  water-level  at 
the  onset  of  drawdown  results  in  small  hydraulic  head  and 
discharge  rates.   Low  discharge  rates  depress  removal  of 
resuspended  material. 

Concentrations  of  POM,  TKN,  and  TP  in  the  lake 
discharge  are  shown  in  Figure  3-4d,  e,  and  f.   The  data 
points  represent  the  average  value  of  two  replicate 
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samples.   "Error  bars"  averaged  ±0.40  mg/1  for  POM  (range 
0-1.63),  ±0.15  mg/1  for  TKN  (range  0-0.42),  and  ±0.01  mg/1 
for  TP  (range  0-0.03).   The  concentrations  in  the  discharge 
were  higher  during  the  first  month  of  drawdown  compared 
with  the  pre-drawdown  period  (Figure  3-4d,  e,  and  f ) . 
These  differences  are  statistically  significant  (P<0.05, 
P<0.05,  P<0.01,  respectively)  using  a  pooled  analysis  of 
variance  (Byrkit  1975)  and  assuming  no  autocorrelation 
among  data  points  (Table  3-4) .   After  this  period,  when 
discharge  decreases  to  about  twice  the  assumed  base- 
discharge  (see  Figure  3-3)  concentrations  of  POM,  TKN,  and 
TP  drop  to  near  predrawdown  levels.   At  these  low  lake 
stages,  the  outflow  from  Newnan's  Lake  may  be  restricted  to 
the  less-turbid  surface  layer  of  water  released  through  the 
shallow  creek. 


Table  3-4.   Concentrations  of  particulate  organic  matter 
(POM  in  mg/1),  total  Kjeldahl  nitrogen  (TKN  in  mg/1),  and 
total  phosphorus  (TP  in  /xg/1)  in  surface  discharge  from 
Newnan's  Lake  prior  to  and  during  the  first  month  of 
drawdown.   Number  given  are  the  mean,  number  of 
measurements  (N) ,  and  standard  error  of  the  mean  (S.E.) 


Predrawdown 

During  1st  month 
of  drawdown 

Mean 

N 

S.E. 

Mean 

N 

S.E. 

POM 

9.48 

22 

0.25 

13.61 

8 

0.47 

TKN 

1.72 

22 

0.05 

3.02 

8 

0.03 

TP 

55.81 

22 

2.88 

83.13 

8 

11.42 
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Considerable  changes  in  water  quality  of  the  outflow 
occurred  during  the  sampled  storm  events.   Three  sampled 
storm  events  during  the  drawdown  period  produced  increases 
in  the  concentration  of  POM,  TKN,  and  TP  (Figure  3-4d,  e, 
and  f )  .   Other  peaks  in  these  time  patterns  may  have  been 
associated  with  wind  events  which  were  not  sampled. 
Although  a  more  detailed  analysis  (incorporating  daily  wind 
data  and  in-lake  water  quality)  would  provide  better 
evidence,  it  appears  that  sampled  storms  were  effective  in 
resuspending  surficial  flocculent  sediment  in  Newnan's 
Lake.   This  corresponds  with  findings  in  other  shallow  lake 
systems  (Sheng  and  Lick  1979;  Somlyody  1982).   Mixing  of 
these  deposits  in  the  water  column  enhances  their  removal 
through  flushing.   The  lack  of  high-intensity  wind  events 
(e.g.  in  excess  of  30  km/hr)  during  the  study  period 
contributed  to  the  low  particle  flushing  rates  encountered. 
Two  storms  of  low  intensity,  sampled  prior  to  removal  of 
the  spillway,  produced  small  or  no  increases  in  particulate 
matter  and  nutrient  concentrations  of  the  flow  through  the 
section  of  the  dam  with  the  missing  top-f lashboard. 
Evidently,  storms  are  effective  in  flushing  particulates 
and  nutrients  from  the  system.   The  flushing,  however,  is 
typically  significant  when  obstructions  such  as  the 
spillway  are  removed  from  the  lake  outlet. 
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Oxidation  of  Organic  Matter  in  Exposed  Littoral  Sediments 

A  second  objective  of  the  short-term  drawdown  of 
Newnan's  Lake  was  to  improve  littoral  zone  habitat. 
Exposure  of  littoral  zone  lake  bottom  to  air  may  allow 
oxidative  removal  of  organic  material  (Wegener  et  al.  1974) 
and  should  consolidate  flocculent  sediments  (Fox  et  al. 
1977) .   This  improves  the  area  for  future  sportfish 
spawning  and  provides  firm  rooting  for  aquatic  macrophytes 
(Holcomb  and  Wegener  1971) ,  which  may  contribute  to  a 
higher  standing  crop  of  aquatic  macro-invertebrates  and, 
eventually,  fish  (Wegener  and  Williams  1974).   A  vegetated 
littoral  zone  may  function  as  an  effective  nutrient  trap 
and  reduce  nutrient  input  from  runoff  into  the  lake  (Mickle 
and  Wetzel  1978)  .   Numerous  authors  have  suggested  that 
aquatic  macrophytes  can  inhibit  the  development  of  algae 
(Canfield  et  al.  1984;  Crisman  1986b).  Accumulation  of 
algae  are  perceived  as  a  persistent  problem  in  Newnan's 
Lake. 

No  significant  net  oxidation  of  littoral  substrate  was 
noted  after  28  days  of  in  situ  exposure  (Figure  3-5) .   The 
difference  in  organic  matter  content  between  pre-  and 
postexposure  enclosures  was  statistically  not  significant 
(paired  t-test,  P=0.05).   Field  observations  showed 
production  of  organic  matter  in  the  form  of  germinating 
seeds,  roots,  and  above-ground  plant  biomass  inside  the 
enclosures.   If  oxidation  of  the  sediments  occurred,  it  may 
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have  been  masked  by  this  organic  matter  production.   It  is 
possible  that  oxidation  rates  may  have  been  quite  low, 
particularly  if  the  organic  sediments  were  largely  humic 
compounds,  which  are  relatively  resistant  to  microbial 
degradation  (Sederholm  et  al.  1973) . 

The  laboratory  experiments  did  not  provide  evidence  of 
oxidation  of  these  sediments  either.   Oxidation  rates  were 
extremely  low  and  not  significantly  different  between  dry, 
moist,  and  permanently  inundated  substrate  (Table  3-5) . 


Table  3-5.   Mean  oxidation  rates  of  incubated  littoral 
sediments  in  the  laboratory  under  3  different  treatments. 


N    Mean  oxidation  rate  S.D.    Range 
(mg  g_1d-1) 


Inundated      9  0.23  0.05  0.17-0.30 

Moist           9  0.24  0.05  0.16-0.35 

DrY Z °-19  0.05  0.13-0.26 


Production  of  organic  matter  in  these  laboratory  chambers 
under  the  low-light  regime  was  less  likely  than  in  the 
field  enclosures.   Hence,  these  experiments  suggest  that 
oxidation  rates  were  indeed  low.   These  results  support 
findings  by  Fox  et  al.  (1977),  who  noted  no  significant 
decomposition  of  organic  matter  in  sediment  from  Lake 
Apopka,  Florida,  using  a  series  of  laboratory  experiments. 
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Consolidation  of  Littoral  Sediments 

In  the  laboratory  consolidation-experiment,  water 
content  of  exposed  sediments  decreased  from  90  to  50 
percent  (Figure  3-6)  and  remained  moderately  compact  149 
days  after  reflooding.   Organic  matter  content  remained 
constant  throughout  the  duration  of  the  experiment  (Figure 
3-6) .   Firmer  substrate  will  reduce  rates  of  resuspension 
during  periods  of  high  winds. 

Summary  and  Conclusions 

The  drought  of  1989  produced  a  low  lake-stage  (41  cm 
below  average)  at  the  start  of  the  drawdown  period  and  a 
low  incidence  of  storm  events  during  the  drawdown.   This 
reduced  the  discharge  rates  of  particulate  organic  matter 
and  nutrients  from  the  lake. 

Concentrations  of  POM,  TKN,  and  TP  in  the  lake 
discharge  were  significantly  higher  during  the  first  month 
of  drawdown  than  during  the  pre-drawdown  period.   These 
concentrations  dropped  to  near  predrawdown  levels  at  lower 
lake  stages,  when  the  sill  depth  at  the  mouth  of  Prairie 
Creek  may  have  restricted  the  outflow  from  the  lake  to  the 
less-turbid  surface  layer  of  water. 

Storm  events  produced  increases  in  concentration  of 
particulate  organic  matter  and  nutrients  in  the  discharge 
from  the  lake  during  drawdown.   Storms  sampled  prior  to 
opening  of  the  spillway  did  not  cause  such  increases. 
Storms  resuspended  fine  particulate  deposits  in  this 
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Figure  3-6.   Changes  in  water  and  organic  matter  content 
(weight/weight  ratio)  in  Newnan's  Lake  littoral  sediments 
upon  exposure  and  after  reflooding  in  a  laboratory  setting. 
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shallow,  exposed  lake  and  promoted  flushing  of  this 
material. 

The  sill  at  the  mouth  of  Prairie  Creek  reduces  the 
likelihood  of  a  gravity  drawdown  to  a  stage  much  lower  than 
that  accomplished  (19.70  m  MSL)  .   At  this  stage,  the  water 
level  is  barely  lakeward  of  the  cypress  tree  fringe  and 
drying/ consolidation  of  lake  bottom  is  limited  to  a  narrow 
littoral  zone  fringe. 

Field  and  laboratory  tests  did  not  show  oxidative 
removal  of  organic  matter  from  exposed  areas  of  the  lake 
bottom.   Consolidated  sediments  remained  moderately  firm 
after  reflooding  in  a  laboratory  experiment  and  may, 
therefore,  provide  improved  substrate  for  rooted  aguatic 
vegetation. 


CHAPTER  4 

REDISTRIBUTION  OF  ORGANIC  SEDIMENTS 

FOLLOWING  A  SHORT-TERM  DRAWDOWN. 


Water-level  drawdown  is  a  well  established  lake 
management  technique.   It  has  been  used  to  influence  the 
abundance  and  composition  of  aquatic  plant  communities 
(Holcomb  and  Wegener  1971;  Hestand  and  Carter  1975;  Cooke 
1980;  Tarver  1980),  increase  fish  standing  crop  (Lantz  et 
al.  1964;  Wegener  and  Williams  1974),  and  consolidate 
littoral  sediment  (Holcomb  et  al.  1975;  McKinney  and 
Coleman  1980) .   Drawdowns  have  also  produced  higher 
densities  of  littoral  macroinvertebrates  (Wegener  et  al. 
1974)  ,  increased  nutrient  concentrations  in  the  water 
column  (Serruya  and  Pollingher  1977) ,  and  reduced  dissolved 
oxygen  levels  by  disrupting  the  thermal  stratification  of 
the  water  column  (Richardson  1975) .   The  intent  of  the  1989 
short-term  drawdown  in  Newnan's  Lake  was  to  improve 
littoral  habitat  for  fish  growth  and  recruitment  by 
allowing  oxidation  and  consolidation  of  exposed  littoral 
zone  lake-bottom  (Chapter  3). 

In  spite  of  the  frequent  use  of  water-level  drawdown, 
little  is  known  about  its  impact  on  erosion  of  littoral 
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sediments.   Low  water  levels  during  drawdown  increase  wind- 
wave  action  on  littoral  substrate.   The  resulting  physical 
resuspension  of  bottom  material  may  then  be  followed  by 
transport  and  gravitational  settling  of  the  entrained 
particles  in  deeper  areas  (Bengtsson  et  al.  1990) .   As 
such,  material  can  be  focused  from  the  littoral  zone  to 
profundal  substrate  (Davis  1968) .   This  process  may  be 
particularly  pronounced  in  lakes  with  fine,  organic  bottom 
substrate,  i.e.  those  where  water-level  drawdown  is  most 
commonly  applied. 

The  effects  of  enhanced  resuspension  of  bottom 
material  may  include  increased  turbidity  and  reduced  light 
penetration  in  the  water  column.   It  may  also  exert 
considerable  oxygen  demand  (James  1974)  and  lead  to 
increased  availability  of  nutrients  in  the  overlying  water 
for  algal  utilization  (Holdren  and  Armstrong  1980;  Pollman 
1983).   "Sloughing"  of  eroded  material  to  deeper  parts  of 
the  lake  reduces  maximum  water  depth  and  may  lead  to  the 
development  of  substantial  shallow  areas.   This  increases 
the  potential  of  the  lake  ecosystem  to  support  extensive 
growth  of  rooted  macrophytes. 

The  objective  of  this  part  of  the  study  was  to 
guantify  the  removal  of  flocculent  sediments  from  the 
littoral  zone  and  determine  their  redistribution  to  deeper 
areas  of  the  lake  during  a  short-term  drawdown  of  Newnan's 
Lake. 
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Materials  and  Methods 
Sediment  Cores 

Removal  and  redistribution  of  lake  bottom  material  due 
to  drawdown  was  investigated  using  two  series  of  nine 
sediment  cores.   Marker  horizons  were  identified  in  the 
cores  either  from  profiles  of  bulk  density,  unsupported 

Pb,  ■LJ/Cs,  organic  matter  and  nutrient  content  or  from 
direct  field  evidence  of  distinct  stratigraphy.   Gain  or 
loss  of  sediment  at  each  site  was  then  quantified  by 
matching  marker  horizons  between  pre-  and  postdrawdown 
cores. 

Figure  4-1  shows  the  location  of  the  cores.   An 
electronic  long-range  navigation  system  (LORAN,  Si-Tex  797) 
was  used  to  ensure  agreement  in  location  of  pre-  and 
postdrawdown  sampling  sites.   Additionally,  triangular 
compass  measurements  with  permanent  landmarks  were  used  for 
the  littoral  cores.   Specifics  of  the  core  locations  are 
given  in  Table  4-1.   Latitude  and  longitude  records  were 
relative  to  the  calibration  site  located  near  the 
southwestern  boat  ramp  (Figure  4-1) . 

The  accuracy  of  this  LORAN  to  return  to  a  sampling 
site  is  limited  to  a  range  of  approximately  20  meters. 
This  inherently  results  in  error  when  comparing  pre-  and 
postdrawdown  sediment  stratigraphy.   This  error  is  reduced 
when  in-lake  variability  between  nearby  sites,  in  terms  of 
water  depth  and  bottom  stratigraphy,  is  low  as  in  Newnan's 
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Figure  4-1.    Map  of  Newnan's  Lake  with  core  locations 
indicated  (filled  circles) . 


63 
Lake  with   its  rather   flat  bottom  topography    (Holly   1976) 
and  homogeneous  deposits  of   soft  sediment    (Skoglund   1990) . 
By  matching  several  marker  horizons  between  cores,    rather 
than  one  single  horizon,    error  may  be  reduced  further. 


Table  4-1.      Newnan's  Lake  cores:      Water  depth 
and  location. 


Core 

Depth  (cm) 

Latitude 

Longitude 

1-11 

150 

29 

39 

03 

82 

13  82 

2-12 

200 

29 

39 

12 

82 

12  87 

5-15 

170 

29 

38 

02 

82 

12  71 

6-16 

150 

29 

38 

04 

82 

13  47 

7-17 

118 

29 

37 

69 

82 

14  74 

8-18 

115 

29 

36 

80 

82 

14  50 

NE 

65 

29 

39 

27 

82 

12  05 

SE 

62 

29 

36 

68 

82 

14  43 

W 

58 

29 

38 

16 

82 

14  22 

Note:      Measured  with  Si-Tex   797   LORAN   in  degrees,    minutes, 
and    (minutes/100).      Calibration  marker   at   southwest 
boat   ramp:    29   37   07/82   15   23. 


Water  depth  at  all  stations  was  estimated  by  carefully 
lowering  a  Secchi  disk  until  contact  between  the  disk  and 
soft  bottom  substrate  was  noticed. 

Cores  were  collected  and  preserved  as  described  in 
Chapter  2 .   Visual  observations  of  the  cores  were  made  in 
the  field  and  again  during  sectioning  in  the  laboratory  to 
detect  changes  in  sediment  color  or  texture  which  could 
serve  as  marker  horizons. 
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Bulk  Density  and  Nutrient  Analyses 

The  cores  were  sectioned  in  1  cm  intervals  to  a  depth 
of  30  cm,  into  2  cm  intervals  from  30  to  80  cm  depth,  and 
into  4  cm  intervals  below  that.   The  sections  were  sealed 
in  plastic  bags  and  stored  in  a  refrigerator.   Bulk 
density,  organic  matter  content,  total  Kjeldahl  nitrogen 
(TKN) ,  and  total  phosphorus  (TP)  were  determined  as 
described  in  Chapter  2.   Measurements  were  made  at 
intervals  selected  to  give  a  representation  of  the  entire 
core  profile. 
Radio-isotope  Analyses 

Marker  horizons  from  measurement  of  210Pb  and  137Cs 
levels  throughout  the  core  profiles  may  aid  in  a 
determination  of  sediment  redistribution  during  drawdown. 
210Pb  profiles  have  been  reliable  in  documenting  sediment 
removal  due  to  major  storm  surges  (Robbins  et  al.  1978). 
Pennington  (1981)  recorded  variations  in  137Cs  profiles  in 
a  shallow  lake  resulting  from  episodic  sediment 
redistribution  and  deposition.   Bengtsson  et  al.  (1990) 
successfully  used  settling  sediment  traps  to  investigate 
redistribution  of  fine  sediments  in  Swedish  lakes. 
Maintaining  suspended  sediment  traps  in  the  water  column, 
however,  was  not  an  option  in  a  public  use  lake,  such  as 
Newnan's.   Use  of  210pb  and  137Cs  in  the  determination  of 
depositional  markers  has  additional  significant  benefits, 
because  such  profiles  permit  the  calculation  of  the  age  of 
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deposited  material  in  the  cores  (Eakins  and  Morrison  1978; 
Appleby  and  Oldfield  1983) . 

Between  8  and  16  210Pb  and  137Cs  measurements  were  made 
depending  on  the  length  of  the  core.   210Pb  and  137Cs 
concentrations  were  measured  by  direct  y-assay  as  described 
in  Chapter  2. 

Results  and  Discussion 
Littoral  Cores 

Two  series  of  three  littoral  cores  were  taken; 
northeast,  southeast  and  west  cores  (Figure  4-1). 
Predrawdown  cores  were  taken  when  water  depth  averaged  62 
cm.   Following  the  drawdown  and  natural  rise  of  the  water 
level  to  predrawdown  stage  the  second  series  of  three  cores 
were  taken  at  the  same  locations.   Water  level  at  these 
sites  was  lowered  from  62  to  30  cm  for  8  weeks  during  the 
drawdown.   Sand  underlying  flocculent  brown  substrate 
served  as  a  suitable  marker  horizon. 

Sediment  depth  to  sand  decreased  by  an  average  of  42% 
after  drawdown  (Figure  4-2a) ,  and  bulk  density  of  the 
remaining  substrate  increased  by  an  average  of  250%  at  all 
three  sites  (Figure  4 -2b) .   Organic  matter,  TKN,  and  TP 
(expressed  as  weight /weight  ratios)  decreased  at  two  of  the 
three  littoral  zone  sample  sites  following  drawdown  (Figure 
4-2c,  d,  and  e) .   Analogously,  the  amount  of  organic  matter 
overlying  the  sandy  littoral  bottom  was  reduced  (Figure  4- 
2f ) .   Calculated  removal  rates  are  given  in  Table  4-2. 
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Figure  4-2.  Northeast,  southeast  and  west  littoral  cores, 
Newnan's  Lake.  Predrawdown  (open  bars)  and  postdrawdown 
(filled  bars) ;  a)  sediment  depth  to  sand  (cm) ,  b)  bulk 
density  (g/cm3) ,  c)  organic  matter  content  (g/g) ,  d)  total 
Kjeldahl  nitrogen  (mgN/gdry  ^j ,  e)  total  phosphorus 
(mgp/gdry  ^J ,  and  f )  organic  matter  to  sand  (g/cm2) . 
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Table  4-2.   Material  removal  rates  (g  m~2day-1) 
from  littoral  cores  during  drawdown,  Newnan's 
Lake. 


Removal  rates NE SE W 

Organic  matter      3.52      14.71       0.0 
TKN  0.14        1.58         0.01 

TP  0.01        0.11        -0.02 


Removal  rates  from  the  west  core  deviated  from  the 
observed  pattern.   This  may  have  been  due  to  the  location 
of  this  station  in  a  rather  quiet  littoral  cove,  less 
subject  to  wind -wave  action.   Sediment  depth  to  sand  (in 
cm)  decreased  at  this  site  following  drawdown,  but  no 
removal  of  organic  matter  (in  mg/cm2)  was  noted  (Figure  4-2f) 

These  observations  may  be  interpreted  in  two  ways. 
First,  sediment  thickness  may  have  decreased  without  actual 
removal  of  material  due  to  consolidation  of  the  substrate 
at  this  site.   This  appears  unlikely,  since  this  location 
remained  inundated  by  at  least  30  cm  of  water  at  all  times. 
Second,  erosion  of  material  did  in  fact  occur,  resulting  in 
the  observed  increased  bulk  density  of  the  remaining 
substrate  at  this  site.   However,  eroded  organic  material 
may  have  been  replaced  by  net  primary  production  during  the 
period  of  drawdown.   Increased  organic  matter  and  total 
phosphorus  content  at  this  site  (Figure  4-2c,  and  e) 
support  the  latter  interpretation.   Enhanced  primary 
productivity  of  benthic  algae  is  likely  at  low  water  depth 
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during  the  drawdown  when  more  light  reaches  the  littoral 
bottom.   This  is  particularly  plausible  at  the  protected 
west  core  site,  with  lower  light  inhibition  by  suspended 
solids. 

Removal  of  littoral  sediment  likely  resulted  from 
increased  shear  stress  caused  by  wind-wave  action.   Such 
surface  waves  generate  periodic  oscillations  in  the  water 
column,  which  attenuate  with  water  depth  (cf.  Wetzel  1983) 
but  may  reach  the  sediment -water  interface  in  shallow 
water.   When  this  stress  exceeds  the  bulk  shear  strength  of 
the  surficial  deposits,  sediment  resuspension  occurs  (Lick 
1982) .   Waves  are  largely  wind-induced,  although  wakes  from 
boat  traffic  may  also  be  considerable  in  this  public-use 
lake.   Resuspended  material  may  then  be  moved  from  shallow 
to  deeper  parts  of  the  basin  by  water  currents  (Davis  and 
Ford  1982;  Hakanson  1982;  Bengtsson  et  al.  1990)  and 
eventually  settle  where  water  depth  is  sufficient  to 
eliminate  stress  from  wind-wave  action. 

The  amount  of  organic  matter  was  not  significantly 
different  between  the  pre-  and  postdrawdown  cores  although 
the  northeast  and  southeast  cores  demonstrated  a 
considerable  reduction  (Figure  4-2f ) .   The  reduction  in 
thickness  of  the  sediment  layer  and  the  increase  in  bulk 
density  following  drawdown,  observed  at  all  three  sites, 
were  statistically  significant  (P=0.05,  paired-t) .   The 
assumption  in  the  study  was  that  these  three  sites  were 
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representative  of  the  littoral  zone.   A  different  approach 
may  be  to  make  many,  simple  to  perform  measurements  of 
sediment-thickness-to-sand  throughout  the  littoral  zone 
pre-  and  postdrawdown .   Such  sampling  will  represent  better 
the  entire  littoral  zone,  but  it  will  not  provide  data  on 
removal  rates  of  bulk  sediment,  organic  matter  and 
nutrients. 
Profundal  Cores 

The  locations  of  the  profundal  core  sites  are  shown  in 
Figure  4-1.   Material  in  the  cores  consisted  generally  of 
homogeneous,  black-brown  sediment  of  low  bulk  density  (<100 
mg/cm3  up  to  a  depth  of  30  cm) . 

Cores  l  and  11.   The  matching  profiles  for  cores  1 
(predrawdown)  and  11  (postdrawdown)  for  bulk  density, 
organic  matter,  unsupported  210Pb,  137Cs,  TKN,  and  TP  are 
shown  in  Figure  4-3  (left  panel).   During  sectioning  in  the 
laboratory,  the  first  occurrence  of  clay  was  noted  at  120 
cm  (core  1)  and  132  cm  (core  11)  (Table  4-3) .   This 
suggested  an  addition  of  12  cm  of  material  to  the  profile 
at  this  station  during  drawdown.   Profiles  for  bulk  density 
and  organic  matter  content  revealed  a  marker  at  a  depth  of 
82  cm  (core  1)  and  94  cm  (core  11) .   This  also  indicated  a 
gain  of  12  cm  during  drawdown. 

Unsupported  210Pb  and  137Cs  profiles,  however,  implied 
removal  of  material  postdrawdown.  Profiles  for  TKN  and  TP 
were  inconclusive  with  widely  fluctuating  nutrient  levels 
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Figure  4-3.  Depth  profiles  for  Newnan's  Lake  cores.  Left 
panel;  cores  l  (predrawdown,  dashed  line)  and  11  (post- 
drawdown,  solid  line) .  Right  panel;  cores  2  (predrawdown, 
dashed  line)  and  12  (postdrawdown,  solid  line) .  Arrows 
indicate  horizons  discussed  in  the  text.  A,  bulk  density  in 
mg/cm  ;  B,  organic  matter  content  in  mg/mg;  C,  unsupported 
Pb-210  in  pCi/g;  D,  Cs-137  in  pCi/g;  E,  total  Kjeldahl 
nitrogen  in  mg/gdry  wt;  and  F,  total  phosphorus  in  mg/gd   wt. 
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in  both  cores.   Average  TKN  concentrations  in  the  top  120 
cm  of  cores  1  and  11  were  33  and  27  mgTKN/gdry  weight 
respectively.   TP  levels  averaged  3.4  and  3.0  mgTP/gdry  weight, 
Examination  of  the  six  profiles  did  not  demonstrate  a 
distinct  removal  or  gain  of  bottom  material  due  to  drawdown 
for  this  station. 


Table  4-3.   Length  of  retrieved  cores  and 
depth  of  sand  and  clay  horizons  in 
Newnan's  Lake  cores. 


Core 

Core 

length 

Depth  (cm)  to 

site 

( 

cm) 

sand(s)/clay(c) 

1 

129 

120c 

11 

148 

132c 

2 

155 

12 

140 

5 

113 

85b, 113s, 95c 

15 

92 

76s, 86c 

6 

135 

135c 

16 

115 

7 

135 

100s 

17 

148 

100s 

8 

72 

38s, 64s 

18 

61 

36s 

Cores  2  and  12.   At  station  2,  bulk  density  and 
organic  matter  profiles  (Figure  4-3,  right  panel)  revealed 
clear  markers  at  68  cm  predrawdown  (core  2)  and  78  cm 
postdrawdown  (core  12) .   Peaks  and  dips  in  the  unsupported 
210Pb  and  137Cs  profiles  displayed  a  similar  gain  of 
approximately  10  cm  of  material  in  the  postdrawdown  core. 
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Predrawdown  peaks  in  the  210Pb  profile  at  depths  of  10,  18, 
and  44  cm  occurred  in  the  postdrawdown  core  at  16,  24,  and 
54  cm.   A  similar  gain  was  evident  from  the  137Cs  profiles. 
This  indicated  an  addition  of  approximately  9  cm  during 
drawdown  at  this  station. 

Based  on  the  bulk  density  of  the  top  10  cm  layer  of 
the  postdrawdown  core,  this  gain  was  equivalent  to  0.18 
g/cm2.  This  translated  into  an  additional  1.1  years  of 
sediment  deposition  during  the  245-day  drawdown  period 
compared  to  normal,  non-drawdown  sedimentation  rates 
(discussed  later) .   TKN  and  TP  profiles  were  inconclusive. 
TKN  levels  fluctuated  between  25  and  35  nig-n^/g^y  weight 
with  a  few  measurements  below  20  1^.^^^  weight.   Averages 
for  both  cores  approximated  32  mgTKN/gdry  weight  and 
corresponded  to  levels  in  cores  1  and  11.   Fluctuations  in 
TP  were  even  more  pronounced  with  an  average  concentration 
between  3.0  and  3.5  mgTp/gdry  weight. 

Cores  5  and  15.   Visual  analysis  of  cores  5 
(predrawdown)  and  15  (postdrawdown)  during  sectioning 
revealed  distinct  sand  layers  at  depths  of  85  and  113  cm 
for  core  5  and  at  76  cm  for  core  15.   First  occurrence  of 
clay  was  at  95  cm  (core  5)  and  86  cm  (core  15)  (Table  4-3) . 
Based  on  these  marker  horizons,  9  cm  of  sediments  were 
removed  during  drawdown.   Bulk  density  and  organic  matter 
profiles  appeared  to  display  a  phase  difference  of  a 
minimum  of  10  cm  (Figure  4-4,  left  panel)  also  showing 
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Figure  4-4.  Depth  profiles  for  Newnan's  Lake  cores.  Left 
panel;  cores  5  (pre-drawdown,  dashed  line)  and  15  (post- 
drawdown,  solid  line).  Right  panel;  cores  6  (pre-drawdown, 
dashed  line)  and  16  (post-drawdown,  solid  line) .  Arrows 
indicate  horizons  discussed  in  the  text.  A,  bulk  density  in 
mg/cm  ;  B,  organic  matter  content  in  mg/mg;  C,  unsupported 
Pb-210  in  pCi/g;  D,  Cs-137  in  pCi/g;  E,  total  Kjeldahl 
nitrogen  in  mg/gdry  wt;  and  F,  total  phosphorus  in  mg/gd 
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material  removal  during  drawdown.   Unsupported  210Pb  and 
137Cs  profiles  indicated  a  distinct  removal  of 
approximately  8  cm  during  drawdown.   In  addition,  TKN  and 
TP  profiles  suggested  material  removal  postdrawdown .   This 
implied  that  the  drawdown  removed  a  minimum  of  8  cm,  or 
0.26  g/m2,  from  the  sediment  profile.   This  equated  to 
approximately  4.7  years  of  sediment  removal  at  this  site 
during  drawdown  (discussed  later) . 

The  profiles  also  demonstrated  the  different  character 
of  the  substrate  at  station  5-15  compared  with  other 
sampling  sites.   Bulk  densities  were  significantly  higher 
throughout  the  profiles,  while  organic  matter  and  nutrient 
concentrations  were  depressed  compared  with  deposits  at 
similar  depths  in  other  locations.   Field  observations 
concurred  with  these  data  in  that  substrate  appeared 
firmer,  with  a  relatively  high  sand  content. 

Cores  6  and  16.   Sectioning  of  the  cores  in  the 
laboratory  revealed  a  clay  horizon  in  core  6  (predrawdown) 
at  135  cm  (Table  4-3),  well  below  the  bottom  section  of 
core  16  (postdrawdown) .   Bulk  density  and  organic  matter 
profiles  (Figure  4-4,  right  panel)  showed  a  rather 
homogeneous  top  100  cm  for  both  cores  containing  flocculent 
material  (bulk  density<100  mg/cm3)  and  an  average  organic 
matter  content  of  60%  (weight /weight  ratio) .   Firmer 
substrate  started  at  105  cm  depth  in  core  6  and  at 
approximately  110  cm  in  core  16.   A  horizon  of  inorganic 
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material  was  encountered  first  at  105  cm  depth  in  core  6. 
This  layer  may  have  occurred  immediately  below  the  last 
sampled  segment  in  core  16. 

Isotope  profiles  for  both  cores  also  displayed  a  shift 
in  phase  of  several  centimeters,  while  nutrient  profiles 
appeared  inconclusive.   Average  TKN  and  TP  levels  coincided 
with  cores  l-ll  and  2-12  (30  mgTKN/gdry  weight  and  3.0-3.5 
m9Tp/9dry  weight)  •   In  summary,  these  profiles  provided  weak 
support  for  the  conclusion  that  5  cm,  or  0.05  g/cm2,  of 
material  was  deposited  at  this  site  during  drawdown.   This 
equated  to  approximately  0.04  years  of  sedimentation  in 
excess  of  "normal"  sedimentation  occurring  during  the 
period  of  drawdown  (discussed  later) . 

Cores  7  and  17.   Core  sectioning  in  the  laboratory 
demonstrated  1  m  of  black,  soft,  organic  muck  overlying  a 
distinct  firm,  sandy  layer  in  both  cores  7  (predrawdown) 
and  17  (postdrawdown)  (Table  4-3) .  Water  and  organic 
matter  content  increased  again  below  this  layer.   This  was 
reflected  in  the  bulk  density  and  organic  matter  profiles 
at  this  core  site  (Figure  4-5,  left  panel). 

The  additional  similarity  of  unsupported  210Pb  and 
nutrient  profiles  for  pre-  and  postdrawdown  cores  suggested 
no  significant  redistribution  of  material  at  this  site. 
The  absence  of  a  137Cs  peak  below  40  cm  in  the  postdrawdown 
core  is  unexplained,  but  has  been  verified  by  repeated 
analysis  of  the  sample  material.   Average  TKN  and  TP 
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Figure  4-5.  Depth  profiles  for  Newnan's  Lake  cores.  Left 
panel;  cores  7  (pre-drawdown,  dashed  line)  and  17  (post- 
drawdown,  solid  line).  Right  panel;  cores  8  (pre-drawdown, 
dashed  line)  and  18  (post-drawdown,  solid  line) .  .Arrows 
indicate  horizons  discussed  in  the  text.  A,  bulk  density  in 
mg/cm  ;  B,  organic  matter  content  in  mg/mg;  C,  unsupported 
Pb-210  in  pCi/g;  D,  Cs-137  in  pCi/g;  E,  total  Kjeldahl 
nitrogen  in  mg/gdry  wt;  and  F,  total  phosphorus  in  mg/gd 
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concentrations  in  the  top  100  cm  corresponded  to  levels 
found  in  other  cores  (except  site  5-15) . 

Cores  8  and  18.   Minor  sand  was  encountered  at  38  cm 
in  core  8  (predrawdown) ,  with  a  clear,  pronounced  sand 
layer  at  64  cm.   Core  18  (postdrawdown)  showed  clear  sand 
at  36  cm  (Table  4-3) .   Both  the  bulk  density  and  organic 
matter  profiles  displayed  these  markers  (Figure  4-5,  right 
panel)  with  bulk  density  values  in  excess  of  1200  mg/cm3 
and  organic  matter  content  less  than  5%.   Likewise,  TKN  and 
TP  concentrations  approached  zero  at  65  (core  8)  and  36  cm 
(core  18) .   Evidence  from  these  four  profiles  suggested 
that  approximately  28  cm  of  material  eroded  from  the 
profile  at  this  station  during  drawdown. 

The  absence  of  the  210Pb  peak  in  the  top  10  cm  of  the 
postdrawdown  core  also  suggested  such  removal,  but  to  a 
much  lesser  extent.   Analogously,  elimination  of  the  137Cs 
peak  in  the  top  20  cm  of  deposits  in  core  18  and  the  onset 
of  detectable  levels  of  this  radionuclide  about  10  cm 
deeper  in  the  predrawdown  profile  implied  removal  of 
material  during  drawdown.   Such  massive  removal  rates, 
however,  seemed  implausible  in  light  of  the  small  flushing 
rates  of  particulate  matter  through  the  nearby  outflow 
recorded  during  drawdown  (Gottgens  and  Crisman  1991) .   It 
is  more  likely  that  either  the  exact  location  of  the  pre- 
and  postdrawdown  cores  did  not  match  or  that  a  significant 
disturbance  of  the  sediment  profile  occurred  at  this 
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station  during  drawdown.   Hence,  no  inference  on  gain  or 
removal  of  bottom  material  at  this  site  was  made. 
All  stations 

Rates  of  removal  or  gain  of  material  for  the  profundal 
sampling  stations  are  summarized  in  Table  4-4.  The  245-day 
time  period  covered  started  immediately 


Table  4-4.   Profundal  cores:  Gain/removal  (-)  rates  of  bulk 
sediment,  organic  matter,  TKN,  and  TP  (g  m"2day-1)  . 


Gain/removal(-)  Core  station 

rate         1-11   2-12   5-15    6-16   7-17   8-18 


Bulk  sediment 

* 

7.2 

-10.6 

2.2 

0.0 

* 

Organic  matter 

* 

4.3 

-5.9 

1.5 

0.0 

* 

TKN 

* 

0.2 

-0.3 

0.1 

0.0 

* 

TP 

* 

0.02 

-0.03 

0.01 

0.0 

* 

*)   No  clear  record  or  conflicting  evidence. 

prior  to  removal  of  the  spillway  and  lasted  until  the 
natural  return  of  the  water  level  to  predrawdown  levels 
following  spillway  re-installation. 

Consequently,  out  of  six  profundal  cores  analyzed  for 
this  study,  two  showed  evidence  of  added  sediment 
deposition  during  drawdown,  one  showed  no  gain  or  loss,  one 
showed  removal,  and  two  were  inconclusive.   No  conclusion, 
therefore,  can  be  reached  by  this  study's  data  regarding 
flocculent  sediment  accumulation  in  the  profundal  zone  of 
Newnan's  Lake  during  the  drawdown. 
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While  sediment  removal  from  the  littoral  zone  was 
indicated  by  the  results  of  this  study,  a  quantitative 
record  of  this  transfer  to  profundal  sites  was  difficult  to 
obtain.   This  is  particularly  true  when  profundal  substrate 
consisted  of  a  thick  pack  of  near  homogeneous  material 
without  clear  marker  horizons.   Furthermore,  the  small  drop 
in  lake  level  during  this  drawdown  did  not  create  large 
areas  of  erosion  in  the  littoral  zone.   This  reduced  the 
magnitude  of  potential  redistribution  of  sediments  in  the 
lake  which  reduced  the  signal  in  profundal  core  profiles. 

Using  the  210Pb  profiles,  sedimentation  rates  were 
computed  for  each  profundal  station.   Calculations  followed 
the  constant  rate  of  supply  model  (Goldberg  1963;  Appleby 
and  Oldfield  1983).  These  accumulation  rates  were  then 
compared  to  the  gain  or  loss  of  material  at  those  sites  due 
to  drawdown.   As  such,  this  gain  or  loss  was  equated  to  a 
time  period  of  "normal"  sedimentation.   For  instance,  the 
gain  of  material  due  to  drawdown  at  station  2-12  was 
equivalent  to  1.1  years  of  sedimentation  (Table  4-5). 
Uncertainty  analysis 

Because  inferences  were  made  from  210Pb  profiles  in 
different  aspects  of  this  work,  an  assessment  of  the  level 
of  confidence  in  these  profiles  is  appropriate.   This  was 
particularly  pertinent  since  profiles  were  established  in 
soft  lake  sediment  with  the  potential  of  disturbance  in  the 
chronology  of  deposits.   Three  independent  observations 
aided  in  such  an  assessment. 
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Table  4-5.      Newnan's  Lake  profundal  cores:      Comparisons  of 
recent    (5  yrs.    B.P.)    dry-sediment  accumulation  rates  with 
gain/ loss    (-)    of  material  during  drawdown,    and  with 
calculated  dry-sediment  accumulation  rates. 


Core 

Recent 

Gain/loss 

Gain/loss 

Calc. recent 

Cumulative 

site 

sed . rt . 

during 

during 

sed . rt . 

residual 
un..2™Pb 
(pCi/cm2) 

(g  cm'^r"1) 

drawdown 
(g/cni2) 

drawdown 
(yrs)  ') 

-2  \ 
(g  cm  yr  ') 

1-11 

0.06 

* 

* 

0.05 

30.2 

2-12 

0.10 

0.18 

1.10 

0.10 

26.4 

5-15 

0.07 

-0.26 

-4.70 

0.07 

17.4 

6-16 

0.08 

0.05 

0.04 

0.07 

28.0 

7-17 

0.06 

0.00 

-0.70 

0.05 

27.0 

8-18 

0.07 

* 

* 

0.07 

18.9 

*)      no  clear  record  or  conflicting  evidence. 

1)  sedimentation  normally  occurring  during  the  245  days  between  pre- 
and  post-measurements   is  subtracted. 

2)  Using  a  model  developed  by  Binford  and  Brenner  (1986).  Measured 
average  cumulative  residual  unsupported  210Pb=24.65  pCi/cm2;  Flux  for 
210Pb-fallout=0.77  pCi  cm'Zyr'1. 


First,  recent  210Pb  based  deposition  rates  were  not 
statistically  different  (P=0.05;  two-tailed  correlated 
test)  from  calculated  values  using  an  earlier, 
independently  developed  model  (Binford  and  Brenner  1986) 
(Table  4-5) .   In  this  model  fallout  210Pb  is  used  as  a 
dilution  tracer  to  compute  net  accumulation  rates  of  any 
material  in  surface  mud  according  to: 

r=  F210Pb  x  A_1  (4-1) 

where  F210Pb  equals  the  flux  of  fallout  210Pb  (pCi  cm"2y_1) 
and  A  is  the  activity  of  210Pb  in  the  sediment  sample 
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(pCi/gdr  weight)  •   Furthermore,  the  different  cores  from 
Newnan's  Lake  had  comparable  210Pb  residuals  (i.e.  total 
residual  unsupported  210Pb  contents)  despite  differences  in 
accumulation  rates  (Table  4-5) .   The  210Pb  residuals  of  the 
cores  reflected  the  210Pb  fallout  from  the  atmosphere. 
Since  this  fallout  lies  in  the  range  0.5-0.9  pCi  cm"2yr-1 
(Nozaki  et  al.  1978),  depending  on  locality,  the  210Pb 
residuals  should  lie  in  the  range  16-30  pCi/cm2  (the  210Pb 
radioactive  decay  constant=0. 03114  yr-1)  .   This 
corresponded  to  the  measurements  in  Newnan's  Lake  cores 
(Table  4-5) . 

Second,  recent  210Pb  based  dry-sedimentation  rates 
correlated  well  with  water  depth  (R2=0.81;  N=5) ,  when 
station  5-15  is  excluded.   The  proximity  of  this  site  to  a 
fish  attractor  (e.g.  submersed  brush  attached  to  an 
anchored  buoy)  and,  hence,  higher  boat  traffic  and  boat 
wake  may  have  disturbed  the  sediments  and  depressed 
sedimentation  rates.   Conseguently,  direct  comparisons 
between  this  sampling  site  and  others  in  the  lake  may  be 
misleading.   The  close  relationship  between  water  depth  and 
material  accumulation  rates  is  well-established  (Evans  and 
Rigler  1980) .   Including  the  5-15  site,  57%  of  the 
variability  in  dry-sedimentation  rates  was  explained  by 
depth  of  the  water  column.   While  this  relationship  does 
not  necessarily  underwrite  the  accuracy  of  recent  210Pb 
levels,  it  does  demonstrate  that  these  levels  correlate 
well  with  each  other  (i.e.  their  precision  is  supported) . 
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Finally,  counting  statistics  and  error  prediction  were 
applied  to  the  recorded  unsupported  210Pb  data  to  compute 
statistical  precision.   This  error  analysis  only  addressed 
internal  uncertainty  associated  with  the  accuracy  of  y-ray 
detection.   It  did  not  consider  uncertainty  controlled  by 
external  factors  such  as  error  associated  with  the  (sub) 
sampling  design,  smearing  of  the  core  (Chant  and  Cornett 
1991) ,  bulk  density  determinations,  and  post  depositional 
mobility  of  constituents  in  the  core  (Anderson  et  al. 
1987) .   "Error  bars"  associated  with  the  experimental  data 
are  illustrated  for  cores  1  and  2  in  Figure  4-6.   Other 
cores  showed  a  similar  magnitude  of  error.   The  length  of 
the  error  bar  eguals  one  standard  deviation  (a)  on  either 
side  of  the  point  (i.e.  68.3%  confidence  limits),  which  is 
standard  practice  in  expressing  uncertainty  in  nuclear 
measurements  (Wang  et  al.  1975).   Because  the  recorded 
counts  in  nuclear  counting  experiments  follow  a  Poisson 
distribution  (Knoll  1979) ,  the  predicted  standard  deviation 
is  the  sguare  root  of  the  mean  number  of  counts.   To  arrive 
at  this  mean  (counts/hr) ,  samples  were  counted  from  14  to 
45  hours  depending  on  sample  weight.   Since 

Net  counts  =  Total  counts  -  Background  counts   (4~2) 
uncertainty  in  the  net  counts  is  propagated  according  to 
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Figure  4-6.  Unsupported  Pb-210  profiles  for  core  1  (top)  and 
core  2  (bottom).  The  length  of  the  error  bar  equals  la  on 
either  side  of  the  data  point. 
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°n  -  V<°t2  +  5?5  (4-3) 

where  an,  at,  and  ab  are,  respectively,  the  standard 
deviations  of  the  net  count,  total  count,  and  background 
count . 

The  effect  of  the  computed  internal  uncertainties  on 
the  interpretation  of  the  profile  in  terms  of  marker 
horizons  or  sediment  accumulation  rates  was  small.   For 
instance,  recent  sediment  accumulation  rates  for  core  1 
varied  between  0.049  and  0.055  g  cm~2yr-1  when  mean 
unsupported  210Pb  concentrations  (5  years  B.P.)  plus  and 
minus  1  a,   respectively,  were  substituted  into  equation  (4- 
1) .   Analogously,  these  rates  for  core  2  ranged  from  0.085- 
0.105  g  cm~2yr-1.   Error  analysis,  such  as  Monte  Carlo 
simulation,  may  be  used  to  estimate  the  effect  of  internal 
uncertainty  on  the  assigned  dates  for  an  entire  core 
profile.   Binford  (1990)  applied  this  technique  to  12  cores 
from  north  Florida  lakes  and  found  95%  confidence  intervals 
ranging  from  about  1-2  years  at  sediments  10  years  of  age, 
10-20  at  100  years,  and  8-90  at  150  year  old  deposits.   No 
estimates  of  the  effect  of  external  uncertainty  on  210Pb 
derived  sediment  chronology  have  been  documented. 

Summary  and  Conclusions 
An  average  of  6.08  g  m"2day-1  of  organic  matter  (0.58 
gTKN  m"2day-1  and  0.03  gTP  m~2day-1)  eroded  from  the  littoral 
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zone  by  lowering  the  water  depth  from  62  to  3  0  cm  for  8 
weeks.   Littoral  sediments  with  low  bulk  density  eroded 
fastest,  and  bulk  density  of  remaining  substrate  increased 
by  an  average  of  250%.   In  the  absence  of  direct  field 
evidence  of  distinct  stratigraphy,  transfer  of  this 
material  to  deeper  areas  of  the  lake  was  measured  using 
sedimentary  cores  with  marker  horizons.   The  latter  were 
derived  from  measurements  of  bulk  density,  organic  matter, 
radionuclides,  and  nutrients.   The  quantitative  record  of 
this  transfer,  however,  was  unclear.   The  small  drop  in 
lake  level  during  drawdown  may  have  contributed  to  this. 
Since  the  drawdown  did  not  create  large  areas  of  erosion  in 
the  littoral  zone,  the  magnitude  of  potential 
redistribution  of  sediments  in  the  lake  and,  thereby,  the 
signal  in  profundal  cores,  is  reduced. 

Two  options  may  now  be  identified  for  further  work. 
First,  measurements  limited  to  the  littoral  zone  can  be 
followed  by  calculations  of  deposition  rates  in  the 
profundal  if  the  extent  of  the  zones  of  erosion,  transport, 
and  accumulation  of  sediment  in  the  lake  are  known. 
Second,  use  of  settling  sediment  traps  to  collect 
resuspended  and  transported  material  may  give  additional 
insight.   In  Newnan's  Lake,  however,  installation  of  these 
traps  on  the  profundal  bottom  is  difficult  because  of  the 
soft  nature  of  the  substrate.   Suspension  of  the  traps  in 
the  water  column  may  jeopardize  public  use  of  the  lake. 
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Instead,  this  technique  may  be  used  in  other  systems. 
Since  resuspension  and  erosion  of  fine  sediments  may 
influence  the  ecology  of  a  lake  through  habitat  alteration, 
release  of  nutrients,  high  turbidity  and  enhanced  sediment 
oxygen  demand,  erosion  processes  must  be  considered  when 
drawdowns  are  attempted. 


CHAPTER  5 
SEDIMENT  AND  PHOSPHORUS  DYNAMICS  MODEL 


Analysis  of  field  data  demonstrated  that  the  spillway 
caused  accelerated  accumulation  of  flocculent  and  nutrient- 
rich  sediment  in  Newnan's  Lake  (Chapter  2) .   Plausible 
consequences  of  such  increased  material  transfer  between 
water  and  sediment  for  shallow,  productive  lake  ecosystems 
were  discussed  in  the  introduction  (Chapter  1) .   They 
included:  (1)  increased  resuspension  of  flocculent 
sediments  by  wind-wave  action  leading  to  high  turbidity  and 
diminished  light  penetration  in  the  water  column,  (2) 
enhanced  biological  and  chemical  oxygen  demand  from  organic 
and  inorganic  suspended  matter  (Hargrave  1969;  James  1974), 
and  (3)  expansion  of  flocculent  substrate  into  the  littoral 
zone. 

Each  of  these  consequences,  in  turn,  has  implications 
for  lake  management  criteria  such  as  water  clarity,  level 
of  primary  production,  extent  of  the  littoral  zone,  lake 
access,  and  standing  stock  of  sportfish.   When  acting 
simultaneously,  however,  their  net  effect  on  these  lake 
management  criteria  is  difficult  to  determine.   A  few 
examples  may  illustrate  this.   Increased  suspended  matter 
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in  the  water  column,  for  instance,  decreases  light 
penetration  and  habitat  for  submerged  vegetation,  but 
accelerated  deposition  of  sediments  and  the  resulting 
increased  rate  of  filling  of  the  lake  may  actually  lead  to 
a  larger  area  that  can  be  invaded  by  bottom  vegetation. 

Increased  suspended  matter  may  also  reduce  water 
clarity  and  decrease  the  aesthetic  value  of  the  lake.   This 
opposes  management  objectives  designed  to  secure  the  lake 
for  water  recreation.   Enhanced  primary  production  may 
follow  when  higher  rates  of  wind-wave  induced  resuspension 
of  bottom  material  increase  internal  loading  of  nutrients 
in  the  lake  through  entrainment  of  nutrient-rich 
interstitial  (pore)  water  from  the  surface  sediments  (Lam 
and  Jaquet  1976)  or  desorption  from  suspended  particles 
(Pollman  1983).   Oxygen  demand  from  increased  suspended 
solids  lowers  dissolved  oxygen  level  in  the  water  column, 
which  may  promote  nutrient  release  from  the  sediments 
(Mortimer  1971;  Theis  and  McCabe  1978). 

Low  concentrations  of  dissolved  oxygen  may  stress  fish 
and  other  heterotrophs  and  selectively  favor  species  with 
physiological  and  behavioral  mechanisms  that  aid  them  in 
surviving  under  such  conditions.   The  Florida  gar 
(McCormack  1967) ,  bowfin  (Johansen  1970) ,  and  bullheads 
(Loftus  and  Kushlan  1987)  are  able  to  utilize  atmospheric 
oxygen,  while  small  fishes  with  upturned  mouths  can  extract 
oxygen  from  the  we 11 -oxygenated  surface  film  (Lewis  1970) . 
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Centrarchids,  a  significant  component  of  the  standing  stock 
of  sportfish  in  Florida  lakes,  are  usually  the  first  fishes 
to  die  under  low-oxygen  stress  (Kushlan  1974)  .   This  would 
conflict  with  lake  management  strategies  aimed  at  promoting 
sportfish.   Expansion  of  flocculent  sediments  into  the 
littoral  zone  may  also  reduce  nesting  habitat  for 
centrarchids  (Bruno  1984),  which  rely  on  firm  substrate  for 
deposition  of  eggs  (Eddy  and  Underhill  1978) . 

Changes  in  components  of  the  Newnan's  Lake  ecosystem 
have  been  documented,  such  as  the  consistent  decline  in  the 
number  of  harvestable  sportfish  per  hectare  from  1850  to 
227  between  1982  and  1988  (Florida  Game  and  Fresh  Water 
Fish  Commission  1982-1989)  .   However,  no  work  has  been  done 
to  analyze  the  response  of  the  lake  to  simultaneously 
acting  factors  and  assumptions  such  as  those  described 
above.   In  other  systems,  dynamic  simulation  models  have 
been  used  successfully  for  such  analysis  by  synthesizing 
the  best  current  understanding  of  succession  in  a  prairie 
ecosystem  (Gutierrez  and  Fey  1975) ,  eutrophication  in  a 
lake  (Anderson  1973),  and  predator-prey  dynamics  in  a 
natural  two-species  system  (Montague  et  al.  1982)  .   The 
computer  model  keeps  track  of  all  interrelationships  in  the 
system  once  they  have  been  identified.   It  then  builds  on 
the  reliable  part  of  our  understanding  of  the  system,  while 
compensating  for  the  "unreliable"  part  by  requiring  an 
explicit  statement  of  assumptions  and  showing  the 
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consequences  of  these  in  a  model  analysis.   Hence,  the 
first  objective  of  this  model  is  to  integrate  existing 
information  concerning  the  Newnan's  Lake  ecosystem  and  to 
program  this  information  on  the  computer  in  the  form  of  a 
testable  hypothesis. 

A  second  purpose  of  the  model  is  to  identify  knowledge 
gaps  that  are  most  critical  in  improving  the  understanding 
of  the  behavior  of  this  ecosystem.   This  is  accomplished  by 
systematically  changing  the  value  of  the  parameters  used  to 
develop  the  model  and  determining  the  effect  of  this 
systematic  change  on  model  output.   Those  parameters 
generating  considerable  changes  in  one  or  more  model  output 
characteristics  require  further  field  and/or  laboratory 
study  to  increase  confidence  in  the  model.   As  such,  the 
model  helps  to  guide  limnological  research,  while,  in  turn, 
research  data  can  be  used  to  further  develop  the  model. 

Finally,  the  model  aims  to  evaluate  the  consequences 
of  alternative  management  actions  in  Newnan's  Lake.   Models 
of  different  hypotheses  of  the  causes  of  past  lake 
responses  (such  as  an  increase  in  the  rate  of  sediment 
deposition)  to  perturbation  (such  as  water-level 
stabilization)  may  be  used  to  make  inferences  about  changes 
to  be  expected  following  future  lake  management  actions. 
Different  hypotheses  by  which  the  dam  may  influence 
sediment  and  detritus  accumulation  are  illustrated  in 
Figure  5-1.   They  include  spillway-induced  changes  in  lake 
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hydrology  (e.g.  water  residence  time),  in  sediment  physics 
(e.g.  sediment  resuspension  and  flushing),  and  in  primary 
production  (e.g.  light  and  nutrient  conditions).   By 
testing  the  effect  of  each  of  these  mechanisms  on  model 
output  and  ranking  them  in  terms  of  risks  and  benefits  to 
intended  uses  of  the  lake,  the  model  can  provide  options 
for  lake  management.   As  an  example,  the  possible  effects 
of  installation  of  a  dam  designed  to  release  bottom-water, 
i.e.  one  which  promotes  sediment  flushing,  can  be  tested  in 
the  model.   Similarly,  the  effects  of  hand-planting  of 
native  aguatic  macrophytes,  a  recently  started  program 
(Krummrich,  Florida  Game  and  Fresh  Water  Fish  Commission, 
pers.  comm.),  or  a  proposed  pumpdown  of  the  lake  (KBN  1991) 
may  be  predicted. 

Methods 
Feedback  Dynamics 

Different  dynamic  modeling  methods  (as  compared  with 
static  modeling  of  flow  analysis,  energy  analysis,  or 
statistical  modeling)  may  be  used  to  develop  and  test 
causal  hypotheses  for  ecosystem  behavior.   Two  of  these 
methods  are  the  energy  modeling  technigues  of  Odum  (1983 
1988)  and  feedback  dynamics,  originally  described  by 
Forrester  (1961,  1968).   Odum's  modeling  approach,  based  on 
energy  flow  in  open  systems,  has  been  used  to  test 
hypothetical  scenarios  relative  to  lake  management  (Gayle 
1975;  Fontaine  1978).   Feedback  dynamics  emphasizes  the 
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identification  of  a  closed  system  of  influences  (i.e.  a 
system  of  feedback  loops)  that  produces  time  patterns  in  an 
ecosystem.   By  limiting  the  analysis  to  an  endogenous, 
closed  system  of  feedback  loops,  attention  is  focused  on 
those  influences  that  are  most  important  in  generating  and 
controlling  system  behavior.   External  features,  such  as 
variables  that  only  vary  as  a  function  of  time  (e.g. 
seasons)  and  random  variation  are  excluded.   They  do  not 
require  an  explanation  and  can  not  be  managed  to  generate 
different  system  behavior.   Feedback  dynamics  methodology 
was  used  in  the  development  of  this  simulation  model. 

The  immediate  goal  is  to  develop  a  "dynamic 
hypothesis",  which  explains  the  observed  pattern  of  changes 
in  the  ecosystem.   The  feedback  relationships  involved  in 
the  hypothesis  are  quantified  using  field  and  laboratory 
data,  literature  values,  theory  and  opinion.   Next,  this 
hypothesis  is  represented  in  equation  form  and  the  solution 
is  simulated  on  the  computer.   This  solution  is  then 
analyzed  to  determine  whether  the  model  is  consistent  with 
theory  and  a  plausible  representation  of  what  is  known  to 
be  true  about  the  ecosystem  of  interest.   If  the  model 
fails  this  test,  the  hypothesis  is  rejected  and  a  modified 
hypothesis  may  be  formulated  and  tested.   In  addition,  a 
sensitivity  analysis  is  done  in  which  the  consequences  of 
inaccurate  parameter  estimates  are  determined.   Parameters 
that  produce  substantial  changes  in  one  or  more  output 
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characteristics  when  altered  from  their  initial  values  are 
key  areas  for  further  field  and  laboratory  research  to 
enhance  confidence  in  the  model.   Finally,  model  validation 
(i.e.  a  rigorous  attempt  to  falsify  the  dynamic  hypothesis) 
is  carried  out.   This  may  involve  field  testing  of  the 
model,  such  as  comparing  the  results  of  a  change  in  the 
model  with  the  outcome  of  the  same  change  made  in  the 
field.  A  flow  chart  of  the  feedback  dynamics  procedure  is 
given  in  Figure  5-2.   The  coupling  between  limnology  and 
modeling  is  emphasized  in  that  field  and  laboratory 
research  assist  in  model  development  and  analysis 
(particularly  in  recording  initial  observations, 
quantifying  the  feedback  relationships,  and  field  testing 
the  model) ,  while  the  model  analysis  identifies  the  most 
sensitive  areas  for  further  research.   As  such,  the  model 
is  a  useful  tool  in  the  analysis  of  complex  ecosystems. 
Dynamic  Hypothesis 

The  overall  hypothesis  is  that  the  spillway  in  the 
outlet  of  Newnan's  Lake  influences  a  system  of  feedback 
loops  in  the  lake  due  to  increased  water  residence  time, 
and  reductions  in  both  bottom-water  drainage  and  amplitude 
of  seasonal  water-level  fluctuations.   These  influences 
contribute  to  the  lake  management  problems  described  in  the 
introduction  of  this  chapter.   Feedback  relationships 
involved  in  the  modeled  hypothesis  for  the  profundal  zone 
are  illustrated  in  Figure  5-3.   Causality  in  the 
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hypothesized  relationships  is  implied  by  showing  the 
influencing  parameter  at  the  arrow  tails  and  the  influenced 
parameter  at  the  arrow  heads.   Positive  and  negative  signs 
at  the  arrow  heads  indicate  whether  the  parameter  at  the 
tail  of  the  arrow  and  the  parameter  at  the  head  of  the 
arrow  change  in  the  same  or  in  the  opposite  direction, 
respectively.   The  overall  sign  of  a  closed  loop  is 
indicated  by  a  sign  within  a  small  circular  arrow  within 
the  loop.   In  a  positive  feedback  loop,  the  original 
direction  of  any  change  in  any  loop  component  is 
reinforced.   Conversely,  a  negative  loop  opposes  such 
change.   The  parameters  used  in  the  simulation  model  are 
summarized  in  Table  5-1.   The  simulation  model  listing  is 
given  in  Appendix  B.   Quick  Basic  4.5  (Microsoft  Corp.)  was 
the  programming  language  used  to  develop  the  computer  code. 
Profunda 1.   Causal  feedback  relationships  hypothesized 
for  the  profundal  zone  of  the  lake  (Figure  5-3)  consist  of 
water-column  and  surface-sediment  components  with  processes 
of  resuspension  and  settling  occurring  at  the  sediment- 
water  interface.   The  importance  of  phosphorus  (P)  in 
controlling  primary  production  in  this  lake  was  documented 
in  Chapter  1.   Therefore,  P  is  used  as  a  tracer  of  material 
in  this  model,  and  components  are  represented  in  grams  of  P 
per  m2  (in  the  surface  sediment)  or  grams  of  P  per  m3  (in 
the  water  column) .   Similarly,  resuspension  and  settling 
are  computed  as  grams  of  P  per  m2  per  year  or  grams  of  P 
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per  m3  per  year,  respectively.   Areal  rates  are  converted 
to  volumetric  rates  by  dividing  by  water-column  depth. 

The  dam  influences  the  modeled  relationships  in  three 
ways  (Figure  5-3)  . 

(1)   It  reduces  the  water  flow  rate  through  the  lake 
by  increasing  the  water  residence  time  (tw) . 


Tw  - 


Lake  Volume    Lake  Volume 


£  Outflows    Et  +  Oaw  +  0™ 


'gw 


Et  (evapo-transpiration)  averaged  106.849  m3/day  for 
Newnan's  Lake,  using  coefficients  from  Opper  (1982)  from 
nearby  Lake  Wauburg.   Ogw  (surface  water  outflow)  was 
computed  from  stage-discharge  relationships  developed 
during  1989  when  the  spillway  was  removed  for  a  period  of  3 
months.   The  spillway  reduced  mean  Osw  from  103.680  to 
51.840  m3/day  (Gottgens  and  Crisman  1991).   0 
(groundwater  outflow)  was  calculated  using  Walton's  (1970) 
equation  for  downward  leakage  (Qc) ; 

Qc  =  -E  x  Ah  (5-2) 

m 

where  p  is  the  coefficient  of  permeability  of  deposits,  m 
is  the  saturated  thickness  of  the  lake  bottom  sediments  and 
Ah  represents  the  average  difference  between  the  lake  level 
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and  the  potentiometric  surface  of  the  Floridan  aquifer. 
Following  work  by  Opper  (1982),  estimates  of  1280  m3/day,  3 
m,  and  4  m  were  used,  respectively,  for  p,  m,  and  Ah.   This 
produces  a  Qc  estimate  of  1706  m3/day.   Substituting  values 
for  Et,  Osw,  and  Qc  into  equation  (5-1)  yields  a  tw  of  0.6 
years  for  the  lake  without  spillway  and  0.8  years  with  the 
spillway  in  place.   This  increase  in  water  residence  time 
enhances  settling  of  suspended  matter.   It  also  increases 
lake  volume,  which  results  in  a  larger  surface  area  and 
hydraulic  head  pressure.   Therefore,  rates  of  evapo- 
transpiration  and  losses  to  groundwater  increase.   Both 
involve  loss  of  water  without  loss  of  particulate  material. 
The  result  is  increased  accumulation  of  flocculent  bottom 
material . 

(2)   The  spillway  reduces  flushing  of  suspended  matter 
by  eliminating  discharge  of  bottom-water.   This  layer  of 
water  is  rich  in  solids  due  to  sediment  resuspended  by  wind 
and  wave-induced  currents  (Sheng  and  Lick  1979) .   Outflow 
over  the  dam  is  restricted  to  less  turbid  surface  waters. 
The  concentration  of  suspended  particulate  matter  in  the 
surface  outflow  of  Newnan's  Lake  averaged  12.5  mg/1  (n=22) 
during  45  days  immediately  prior  to  removal  of  the  spillway 
in  1989  and  18.0  mg/1  (n=8)  for  the  same  length  of  time 
after  the  stoplogs  were  removed  (Gottgens  and  Crisman 
1991) .   Flushing  of  dissolved  material  (such  as  P04)  is  not 
affected  by  bottom-water  drainage  and  is  simply  a  function 
of  water  residence  time  and  P  concentration. 
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(3)   The  spillway  increases  the  lake  stage  and  reduces 
the  amplitude  of  the  seasonal  water-level  fluctuations. 
The  Newnan's  Lake  stage  hydrograph  from  1945  to  present  has 
stages  recorded  between  19.51  and  21.77  m  (msl)  (Adkins 
1991) .   This  record  is  complete  with  the  exception  of  the 
period  1954-57  and  1962-64.   During  these  times  severe 
droughts  reduced  the  lake  levels  in  all  the  lakes  in  the 
Orange  Creek  basin.   Average  stage  for  the  prespillway 
period  (1945-1967)  was  20.24  m  (msl)  with  an  average  annual 
amplitude  of  oscillation  of  0.55  m  (Skoglund  1990;  Adkins 
1991) .   For  the  postspillway  period  (1967-1990)  these 
numbers  were  respectively  20.37  and  0.39.   Higher  lake 
stage  and  smaller  amplitude  of  water-level  fluctuation  both 
increase  water  depth  which,  in  turn,  reduces  the  shear 
stress  from  wind-wave  action  on  surface  sediments. 
Consequently,  the  magnitude  of  resuspension  of  this 
material  decreases. 

Surface  sediments  are  operationally  defined  in  the 
model  as  the  layer  of  material  immediately  below  the 
sediment-water  interface  involved  in  the  settling  and 
resuspension  cycle.   Based  on  examination  of  bottom  cores 
from  Newnan's  Lake,  the  thickness  of  this  "fluid"  zone  was 
estimated  between  3  and  7  cm.   The  actual  thickness  of  this 
layer  may  depend  on  the  intensity  and  frequency  of  wave 
action,  water  depth,  and  the  bulk  density  of  the  bottom 
material.   The  initial  value  of  4.5  cm,  used  in  this  study, 
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corresponds  with  values  reported  by  Hwang  and  Mehta  (1989) 
for  similar  sediment  in  a  non-stratified  Florida  Lake.   An 
increase  in  the  amount  of  surficial,  flocculent  sediment 
will  increase  resuspension.   More  suspended  particles 
(seston)  lead  to  increased  settling  rates  and,  in  turn, 
more  surface  sediments.   This  positive  feedback  loop  is 
illustrated  for  the  inorganic,  live  organic  (benthos  and 
plankton) ,  and  dead  organic  (detritus)  fractions  in  Figure 
5-3  (loops  l,  2, and  3).   The  effect  of  each  of  these 
positive  loops  may  be  bounded  by  negative  loops  between 
seston  and  settling  (loops  4,  5,  and  6)  and  between  surface 
sediments  and  resuspension  (loops  7,  8,  and  9).   Core 
profile  data  were  used  to  partition  surface  sediments  in  an 
organic  and  inorganic  fraction  (0.66  and  0.34, 
respectively)  to  include  the  different  P  dynamics  (e.g. 
uptake,  sorption,  decomposition,  regeneration)  and  settling 
velocities  of  each  fraction  in  the  hypothesis. 

The  P  content  of  each  fraction  was  initially  estimated 
using  Redfield  et  al.  (1963),  whereby  1  gram  of  plankton 
contains  0.012  grams  of  P,  assuming  the  carbon  content  of 
the  plankton  is  50%.   When  back-calculated  from  a  water- 
column  P  concentration  of  0.08  g/m3  (averages  from  Shannon 
1970  and  Canfield  1981)  and  an  average  seston  concentration 
of  20  g/m  (field  data),  a  much  lower  P-content  of  0.004 
gp/g  was  found.   If  80%  of  the  total  seston  is  organic 
(Chapter  3)  and  half  of  this  fraction  is  assumed  to  be 
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planktonic  (i.e.  8  grains  of  plankton  per  m3  and  8  grams  of 
detrital  seston  per  m3)  ,  a  P  content  of  0.006  gp/g  for 
plankton,  0.002  gp/g  for  detritus,  and  0.004  gp/g  for  the 
remaining  4  grams  of  inorganic  seston  per  m3  yield  a  water- 
column  P  concentration  of  0.08  g/m3.   These  P-content 
values  correspond  with  field  data  discussed  by  Fuhs  et  al. 
(1972)  and  Wetzel  (1983)  in  which  C:P  atom-ratios  in  excess 
of  200:1,  i.e.  approximately  double  the  ratio  proposed  by 
Redfield,  are  encountered  commonly. 

Settling  velocities  for  the  inorganic  and  detrital 
fractions  were  taken  from  Sheng  and  Lick  (1979) .   These 
velocities  may  range  over  more  than  four  orders  of 
magnitude  depending  on  particle  size,  concentration  and 
mineral  composition  of  the  suspended  sediments,  and  ionic 
strength  of  the  water  (Fukuda  and  Lick  1980;  Lee  et  al. 
1981) .   In  this  model,  settling  rate  was  assumed  to  be 
proportional  to  the  concentration  of  seston  (Lick  1982; 
Hwang  and  Mehta  1989) .   Plankton  settling  was  estimated  one 
order  of  magnitude  smaller,  due  to  characteristics  that 
improve  flotation  or  reduce  sinking  rates.   Examples  of 
these  characteristics  are  density  reduction  by  the 
production  of  mucilage  sheaths  (most  blue-green 
phytoplankton)  or  gas  vacuoles  (bacteria  and  algae) ,  as 
well  as  active  swimming  (eukaryotic  heterotrophs) . 

Resuspension  of  bottom  material  due  to  wind-driven 
currents  was  assumed  to  be  negligible  in  the  profundal  zone 
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of  the  lake  (Lam  and  Jaquet  1976;  Luettich  et  al.  1990). 
All  entrainment  occurs  as  a  result  of  oscillations  produced 
in  the  water  column  by  wind-wave  action  at  the  water 
surface.   These  surface  waves  produce  a  series  of 
oscillating  motions  downward  in  the  water  column  which  may 
extend  to  the  sediment-water  interface  in  shallow  lakes 
(Pollman  1983).   The  height  of  vertical  oscillations 
attenuates  rapidly  with  depth  with  an  approximate  halving 
of  the  cycloid  diameter  for  every  depth  increase  of  A/9, 
where  A.  is  the  wavelength  of  the  surface  waves  (Wetzel 
1983) .   Near  the  sediment  boundary  layer,  the  oscillatory 
motion  is  reduced  to  a  simple  reciprocating  horizontal 
motion  (U.S.  Army  Coastal  Engineering  Research  Center 
1977) .   When  the  shear  force  of  this  displacement  exceeds 
the  bulk  shear  strength  of  the  surficial  sediment  deposits, 
sediment  resuspension  occurs  (Terwindt  1977;  Hwang  and 
Mehta  1989) .   When  A  and  water  depth  are  known,  this 
displacement  can  be  calculated.   While  not  directly 
proportional  to  wave  height  (H) ,  a  ratio  of  1:20  for  H:A  is 
a  common  average  (Wetzel  1983).   H  is  generally 
proportional  to  the  square  root  of  the  fetch  (x) ,  or 
distance  over  which  the  wind  blows  uninterrupted  by  land. 
The  height  of  the  highest  waves  observed  is  approximated  by 

H  =  0.105  v'x      (Wetzel  1983)  (5-3) 


Ill 

Average  fetch  for  Newnan's  Lake  is  3.5  km.   Hence,  maximum 
wave  height  is  0.62  m  and  A.  is  12.4  m.   Halving  the  cycloid 
diameter  for  every  depth  increase  of  A./9  yields  a  vertical 
displacement  of  0.28  m  at  the  sediment -water  interface  in 
Newnan's  Lake,  with  an  average  water  depth  of  1.5  m.   Field 
data  from  a  similar  shallow,  exposed  lake  (Somly6dy  1982) 
with  a  mean  displacement  of  0.24  m  at  the  boundary  layer 
demonstrated  a  resuspension  flux  of  12,410  x  w  g  «~ay~lf 
where  W  is  the  wind  speed  (m/sec) .   When  applied  to 
Newnan's  Lake  a  flux  of  44,200  g  m^y"1  is  computed, 
assuming  an  average  wind  speed  of  3  m/sec.   Resuspension  is 
assumed  to  be  proportional  to  the  amount  of  surficial 
sediment,  although  it  may  be  influenced  by  composition 
(particle  size,  organic  matter  content) ,  water  content,  and 
activity  of  benthos  (Lee  et  al.  1981).   Surface  sediment 
was  estimated  at  1,300  g/m2  (using  field  data  on  bulk 
density  and  the  thickness  of  this  layer) .   Therefore,  if 
the  total  resuspension  flux  eguals  44,200  g  m"2y_1,  the 
specific  rate  of  resuspension  (SRR)  eguals  34  g  g_1y-1' 
This  value  for  SRR  applies  to  conditions  of  maximum  wave 
height  and  a  water  depth  of  1.5  m.   Since  the  model 
simulates  resuspension  under  average  wave  height,  I  applied 
this  value  of  SRR  to  estimate  resuspension  under  conditions 
of  the  lowest  water  level  (i.e  highest  shear  stress),   in 
the  model,  SRR  varies  with  water  depth  (m)  according  to 
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V  x  Depth  /  r-  .  \ 

SRR  =  a  x  e  (5-4) 

Values  of  36.85  and  -0.084  for  the  model  coefficients  a  and 
Y,  respectively,  result  in  an  SRR  value  of  34  g  g^y"1  at 
minimum  water  depth  and  plausible  resuspension  rates  with 
increasing  water  depth.   The  relationship  between  SRR  and 
water  depth  is  graphed  in  Figure  5-4a.   Entrainment  of 
pore-fluid  soluble-reactive  phosphorus  (SRP)  is  assumed  to 
occur  at  the  same  specific  rate  (Figure  5-3,  loop  10). 
Water-column  SRP  only  settles  after  incorporation  into 
organic  or  inorganic  seston. 

Surface  sediment  is  lost  downward  to  permanent 
sedimentation.   For  each  fraction,  this  loss  rate  is 
proportional  to  the  amount  of  surface  material  present 
(Figure  5-3,  loops  11,  12,  and  13).   The  model  eguations 
for  sedimentation  were  calibrated  to  correspond  with  an 
average  prespillway  sedimentation  rate  of  650  g/m2 
calculated  from  210Pb  dated  profundal  cores.   No  downward 
loss  of  pore-fluid  SRP  is  assumed.   Instead,  this 
phosphorus  may  be  first  absorbed  to  inorganic  particles  or 
taken  up  by  benthos  after  which  sedimentation  as 
particulate  matter  may  occur. 

Direct  absorption  (and  desorption)  of  SRP  from 
sedimentary  detritus  is  assumed  to  be  negligible.   The 
bacterial  coating  of  detrital  particles,  a  component  of  the 
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Figure  5-4.      Plots  of  seven  hypothesized  relationships 
used   in  the  model. 
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benthos,  likely  traps  the  large  majority  of  SRP.   P  derived 
by  the  bacteria  from  detritus  may  flow  back  to  pore-fluid 
SRP  via  P  regeneration.   The  same  pathways  are  assumed  to 
exist  in  the  water  column. 

Absorption  and  desorption  processes  of  P  between 
inorganic  sediment  and  interstitial  pore-fluid  (Figure  5-3, 
loops  14  and  15)  and  between  inorganic  seston  and  water 
(loops  16  and  17)  are  modeled  as  a  function  of  a  specific 
reaction  rate  and  the  extent  to  which  the  respective  P 
concentrations  deviate  from  eguilibrium  conditions.   For 
example,  the  rate  of  P  desorption  from  inorganic  seston  to 
the  surrounding  water  (in  gSrpwater  m_2Y_1)  can  be  expressed 


as 


Sorp  =  v  x  (Pinorg  -  KR  x  Srpwater)  (5"5) 


where  v  is  the  specific  reaction  rate  (gSrpwater  ^Pinorg'V1) 
and  KR  is  the  eguilibrium  constant.   At  eguilibrium, 


p 

v-      _        c  Inorg  . 

Kr  "   c7^ (5~6) 

arPwater 


Using  initial  model  values  for  PINORG  and  SRPWATER,  the 
eguilibrium  constant  for  this  reaction  equals  2.67. 
Analogously,  the  KR  value  for  the  sorption  reaction  between 
P  in  inorganic  sediment  and  pore-fluid  equals  3.06. 


115 
Feedback  loops  involved  in  interactions  between 
plankton,  water-column  SRP,  and  detrital  seston  are  made 
more  apparent  in  Figure  5-5.   The  same  influences  are 
hypothesized  in  the  model  for  interactions  between  benthos, 
pore-fluid  SRP,  and  sedimentary  detritus.   The  difference 
is  the  absence  of  light.   This  reduces  the  rate  of  P  uptake 
in  benthic  algae,  a  process  which  reguires  energy  normally 
supplied  by  photosynthesis  (Kuenzler  and  Ketchum  1962).   To 
prevent  redundancy,  the  following  description  is  limited  to 
the  feedback  relationships  between  plankton,  SRP,  and 
detrital  seston  in  the  water  column. 

In  this  model,  plankton  consist  of  bacteria,  algae, 
and  eukaryotic  heterotrophs .   Since  the  model  emphasizes  P 
dynamics,  these  fractions  were  treated  separately.   Each  of 
the  fractions  may  affect  the  flow  of  P  through  the  system 
in  a  different  way.   P  is  taken  up  by  both  bacteria  and 
algae,  while  being  regenerated  by  eukaryotic  heterotrophs. 
P  is  recycled  from  the  detritus  by  bacteria  and  eukaryotic 
heterotrophs.   Bacteria,  however,  fix  rather  than 
mineralize  P  (Day  et  al.  1989).   Therefore,  no  direct  P 
regeneration  to  the  water  column  from  bacteria  is  assumed 
to  occur.   Direct  P  regeneration  to  the  water  does  occur 
through  ingestion  and  subseguent  excretion  by  eukaryotic 
heterotrophs.   Literature  values  for  similar  lakes  and 
field  data  from  this  study  were  used  to  obtain  the  relative 
size  of  each  fraction.   Again,  a  plankton  concentration  of 
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8  g/m  was  used.   Bacterial  biomass  was  estimated  4  g/m3 
(Saunders  et  al.  1980).   An  algal  biomass  of  3.35  g/m3  was 
derived  using  field  data  on  average  chlorophyll  a  in 
Newnan's  Lake  and  assuming  that  this  pigment  constitutes 
1.5%  of  algal  biomass  (A.P.H.A.  1985).   The  remaining  0.65 
g/m  is  then  assumed  to  be  eukaryotic  heterotrophs .   This 
corresponds  with  data  reported  by  Bays  and  Crisman  (1983) 
for  the  productive  Florida  lakes  from  their  zooplankton- 
biomass  study.   Bacterial,  algal,  and  eukaryotic 
heterotroph  fractions  were,  therefore,  respectively  0.5, 
0.42,  and  0.08. 

Specific  P-uptake  by  phytoplankton  is  affected  by 
light  (e.g.  photosynthetically  active  radiation,  PAR)  and 
concentration  of  SRP  in  the  water  column  (Fuhs  et  al. 
1972) .   More  suspended  matter  in  the  water  column  reduces 
PAR  available  for  photosynthesis  which  lowers  the  specific 
P-uptake  rate  by  phytoplankton  (Figure  5-5,  loop  18) . 
Light  reduction  by  seston  is  assumed  to  be  independent  of 
the  type  of  seston  (i.e.  inorganic,  detrital  or 
planktonic) .   A  reduced  P-uptake  rate  leaves  more  dissolved 
P  (loop  19).   It  also  yields  less  phytoplankton  which,  in 
turn,  depresses  P-uptake  (loop  20) .   As  seston 
concentration  (g/m3)  increases,  available  PAR 
(cal  cm"2hr"1,  at  mean  water  depth)  declines  according  to 
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___  m-  -0.069  x  Seston  /c_-7\ 

PAR  =  16   x  e  (5-7) 


where  16  is  the  value  of  PAR,^  (Gates  1963)  and  -0.069  is 
the  extinction  coefficient  (%  per  g/m3  of  seston)  , 
estimated  using  data  from  Beeton  (1958).   This  relationship 
is  shown  in  Figure  5-4b.  A  hypothesized  effect  of  PAR  on 
the  specific  P-uptake  rate  of  phytoplankton  (PAREFF)  is 
illustrated  in  Figure  5-4c.   Light  of  high  intensity  is 
detrimental  to  many  algae  and  is  particularly  associated 
with  photo-oxidative  destruction  of  enzymes  involved  in 
photosynthesis  (Steemann  Nielsen  1962) .   An  intensity  of  8 
cal  cm~2hr_1  was  used  in  the  model  as  optimum  PAR. 

The  hypothesized  effect  of  SRP  concentration  on 
specific  P-uptake  by  phytoplankton  (PEFF)  is  shown  in 
Figure  5-4d.   It  follows  Michaelis-Menten  kinetics  with  a 
half-saturation  constant  of  0.036  gSRP/m3  (Fuhs  et  al. 
1972)  . 

The  specific  P-uptake  rate,  V  (gSRp  q^^^T1)  , 
may  then  be  given  by 

V  =  v^  x  V  (PAREFF  x  PEFF)  (5-8) 

Thus,  the  specific  P-uptake  rate  is  zero  in  the  absence  of 
light  or  nutrients  and  the  maximum  rate  only  occurs  when 
both  the  light  and  nutrient  regime  are  optimal.   Fuhs  et 
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al.  (1972)  found  v'max  in  excess  of  1000  gSRP  g^a^ton^y"1 
in  a  detailed  study  of  P  kinetics  in  P-starved  planktonic 
diatoms.   They  cautioned,  however,  that  these  laboratory 
results  from  a  few  diatom  species  should  not  be  generalized 
to  algal  assemblages  in  the  field.   Turnover  rates  (i.e. 
specific  uptake  rates)  of  phosphorus  have  been  studied  in 
lakes  of  different  productivity  and  at  different  times  of 
the  year.   The  reported  rates  vary  by  several  orders  of 
magnitude  and  are  faster  under  oligotrophic  conditions. 
Peters  (1975)  reported  P  turnover  rates  ranging  from  7  to 
131,000  g  g'V1  for  epilimnia  of  15  European  lakes.   With 
no  literature  values  specific  for  the  conditions  in 
Newnan's  Lake,  v'max  was  estimated  at  10  gSRp  gp^^^o^Y'1 . 
This  guess  was  based  on  the  low  turnover  rates  found  by 
Peters  (1975)  for  systems  with  high  primary  productivity, 
similar  to  Newnan's  Lake.   The  specific  P-uptake  rate 
multiplied  by  the  algal-P  concentration  (gp/m3)  yields  the 
P-uptake  rate  (gSRP  m"3y-1)  . 

Estimates  for  the  v'max  for  planktonic  bacteria, 
benthic  bacteria,  and  benthic  algae  were  based  on  the  value 
used  for  v'max  for  phytoplankton.   Because  bacteria  are 
capable  of  storing  considerable  P  (represented  in  the  model 
by  an  elevated  P-content)  their  rates  will  be  higher 
relative  to  phytoplankton  specific  P-uptake  (Fuhs  et  al. 
1972) .   The  specific  P-uptake  rate  for  benthic  algae  is  low 
due  to  the  absence  of  light  (Kuhl  1974) .   Michael is-Menten 
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kinetics  were  used  to  describe  the  increase  in  specific  P- 
uptake  rate  by  planktonic  bacteria  with  increasing  SRP  in 
the  water  column  (Figure  5-4e) .   The  resulting  increase  in 
P-uptake  rate,  in  turn,  reduces  available  SRP  (Figure  5-5, 
loop  21) .   It  also  yields  more  planktonic  bacteria,  which 
stimulates  additional  P-uptake  (loop  22) . 

Components  of  the  plankton  and  detrital  seston  are 
linked  in  a  positive  feedback  loop.   As  plankton  increases, 
death  rates  increase  proportionally  and  produce  more 
detritus.   This  leads  to  higher  rates  of  decomposition  and 
an  increase  in  decomposer  biomass  (loop  23) .   Decomposers 
include  bacteria  and  eukaryotic  heterotrophs .   The  effect 
of  this  positive  loop  is  moderated  by  negative  loops 
between  components  of  the  phytoplankton  and  death  rate 
(loops  24,  25,  and  26)  and  between  detritus  and 
decomposition  rate  (loop  27) .  Average  lifetime  (for  all 
fractions  of  plankton  combined)  was  estimated  to  be  0.05 
years  (using  Wetzel  1983) .   This  lifetime  was  halved  for 
the  benthos  because  of  the  larger  fraction  of  short-lived 
microplankton  (i.e.  bacteria).   A  specific  rate  of 
decomposition  of  20  (g  g^y-1)  was  used,  indicating  a 
turnover  time  of  approximately  18  days.   Such  high  rates 
are  plausible  under  conditions  of  year-round  warm 
temperatures  and  correspond  with  data  reported  in  Wetzel 
(1983) . 
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Mineralization  of  P  occurs  during  egestion  and 
excretion  by  eukaryotic  heterotroph  consumers  (Figure  5-5, 
loop  28)  .   A  regeneration  rate  of  50  gSRp  gpiankton'V1  was 
reported  by  Ferrante  (1976)  for  a  mixed  epilimnetic 
zooplankton  assemblage  in  Lake  Michigan  with  water 
temperatures  between  3  and  24 °C.   As  an  approximation,  a 
value  of  500  gSRP  gPiankton~1y~1  was  used  to  account  for 
higher  year-round  temperatures  and  smaller-sized 
zooplankton  in  Newnan's  Lake  (Crisman  1986a). 

Littoral.   Causal  feedback  relationships  hypothesized 
for  the  littoral  zone  of  the  lake  are  shown  in  Figure  5-6. 
The  outflows  from  this  zone  connect  with  the  corresponding 
components  in  the  profundal  zone.   The  emphasis  in  the 
littoral  influence  diagram  is  on  the  effect  of  the  dam  on 
the  flow  of  particles  and  phosphorus  from  the  littoral  to 
the  profundal  area  of  the  lake.   As  discussed  earlier,  the 
reduction  in  amplitude  of  annual  water-level  fluctuations, 
resulting  from  the  dam,  may  gradually  reduce  the  extent  of 
the  vegetated  littoral  zone  along  the  periphery  of  the 
lake.   This  may  reduce  the  ability  of  this  zone  of  plants 
and  associated  epiphytes  to  filter  and  trap  phosphorus  and 
sediments  carried  into  the  lake  from  the  surrounding 
watershed  (Mickle  and  Wetzel  1978).   This,  in  turn,  may 
contribute  to  enhanced  primary  production  in  the  lake. 

As  the  amplitude  of  annual  water-level  fluctuations 
decreases  from  0.55  to  0.39  m  with  the  installation  of  the 
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dam,  the  wetland  fringe  in  the  lake  should  be  reduced  from 
20%  to  perhaps  10%  of  lake  surface  area.   This  reduction  is 
modeled  with  a  half life  of  5  years,  i.e.  5  years  after  the 
dam  is  installed  15%  of  the  lake  surface  area  remains  as 
vegetated  littoral  zone.   The  extent  of  this  vegetated  zone 
(y,  in  %  of  lake  surface  area)  as  a  function  of  amplitude 
of  water-level  fluctuation  (x,  in  m)  is  modeled  as  a 
sigmoid  relationship 


y  = 


-5  07  x  (  5~ 9  ) 

1  +  64.42  x  e 


This  eguation  was  selected  because  it  generates  plausible 
values  for  the  wetland  fringe  of  2%  with  no  water-level 
fluctuations  and  50%  when  the  amplitude  of  fluctuation 
approximates  half  of  the  mean  depth  of  the  lake.   In 
addition,  it  generates  the  proposed  extent  of  wetland 
fringes  under  pre-  and  postspillway  conditions  discussed 
above . 

Littoral  particles  brought  into  the  lake  from  the 
surrounding  watershed  may  either  be  temporarily  trapped  in 
littoral  vegetation  (Figure  5-6,  loop  30)  or  may  pass 
through  this  zone  (loop  29) .   Trapped  particles  erode  from 
littoral  substrate  after  a  delay  time  (0.15  years)  and  move 
to  deeper  areas  of  the  lake  (loop  31),  while  a  small 
fraction  is  lost  to  sedimentation  (loop  32)  .   Both  the  rate 
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of  trapping  and  the  rate  at  which  P  is  filtered  from  the 
overlying  water  by  submersed  littoral  vegetation  and 
epiphytes  (loop  33)  are  assumed  to  be  proportional  to  the 
extent  of  vegetated  littoral  area  (Figure  5-4f ) . 

In  summary,  the  dam  in  the  outlet  of  Newnan's  Lake 
influences  a  system  of  feedback  relationships  in  the  lake 
through  increased  water  residence  time,  and  reductions  in 
bottom-water  drainage  and  amplitude  of  annual  water-level 
fluctuations.   These  influences  contribute  to  the  lake 
management  problems  described  in  the  introduction  of  this 
chapter.   The  relationships  in  this  hypothesis  are  based  on 
data  from  field  and  laboratory  work  or  are  independently 
derived  from  the  literature.   In  the  absence  of  data  or 
theory,  relationships  are  based  on  estimate  or  opinion. 
Many  assumptions  have  been  stated  about  unknown  parameter 
values  and  eguation  forms.   The  following  simulation  and 
model  analysis  subjects  the  assumptions  to  tests  that 
attempt  to  falsify  the  feedback  hypothesis,  identify 
parameters  to  which  the  model  is  sensitive,  and  compare  the 
effect  of  changes  in  the  model  with  similar  changes  made  in 
the  field. 

Results  and  Discussion 
Standard  Run  Output 

The  output  from  the  simulation  is  shown  for  five  model 
variables  (Figure  5-7) .   These  variables  were  chosen 
because  of  their  potential  significance  for  lake  management 
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STANDARD  RUN  OUTPUT  FROM  SIMULATION  MODEL 
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Figure  5-7.   Standard  run  output  from  simulation  model 
with  spillway  installation  at  15  years. 
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(e.g  seston,  water-column  SRP)  and  because  they  may  be 
compared  to  collected  field  data  (e.g.  sedimentary  P  and 
organic  matter  content,  sedimentation  rate) .   The 
simulation  represented  40  years  with  spillway  installation 
initiated  after  15  years.   Dependent  variables  were  scaled 
individually  as  indicated  on  the  vertical  axis.   An 
integration  interval  of  0.0005  years  (of  simulated  time) 
was  selected  to  avoid  model  instability  and  serious 
integration  error,  while  minimizing  computer  run-time.   The 
selected  interval  was  1/4  of  the  inverse  of  the  highest 
specific  rate  used  in  the  model.   The  sinusoid  pattern  of 
the  output  was  produced  by  the  simulated  annual  water-level 
fluctuations. 

Prespillway  seston  concentrations  averaged  28.8  g/m3 
and  increased  by  19%  after  spillway  installation.   The 
model  output  reflected  the  hypothesis  that  the  dam 
increases  concentration  of  suspended  solids.   Separate 
analysis  of  the  pre-  and  postspillway  dynamics  of  the 
inorganic,  planktonic,  and  detrital  fractions  did  not  show 
implausible  shifts  in  the  composition  of  the  seston  (Table 
5-2) .   The  ratios  of  the  average  concentrations  of 
inorganic,  planktonic,  and  detrital  seston  were  1  :  2.4  :  2.1 
and  1  :  2.0  :  1.8  pre-  and  postspillway,  respectively. 

Soluble  P  concentrations  in  the  water  column  increased 
from  an  average  of  0.016  g/m3  prespillway  to  0.028  g/m3 
approximately  10  years  postspillway  installation  after 
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which  the  average  remained  constant  (Figure  5-7) .   This 
corresponded  to  the  pattern  for  the  other  model  output 
variables  shown  for  the  standard  run,  with  the  exception  of 
sedimentary  organic  matter  content.   The  concentration 
range  for  water-column  SRP  compared  well  with  reported 
ranges  of  0.003-0.028  (Shannon  1970)  and  0.01-0.06  g/m3 
(Crisman  1986a) . 


Table  5-2.   Pre-  and  postspillway  concentrations  of 
inorganic,  planktonic,  detrital,  and  total  seston  from 
standard  model  run.   Values  given  are  the  average,  minimum, 
and  maximum  (in  g/m3)  ;  time-of -minimum  and  time-of -maximum 
(in  years  after  model  start) . 


Prespillway  Postspillway 
Avg   Min   Tmin  Max   Tmax  Avg   Mln   Tmin  Max   Tmax 

Inorganic    5.3   4.0   0.0   6.5   0.6   7.3   4.4   15.0   8.2   30.9 
seston 

Planktonic   12.6   8.0   0.0   14.4   6.8   14.6   11.6   16.3   17.1   32.1 
seston 

Detrital    10.9   8.0   0.0   11.7   5.9   12.9   10.7   16.3   14.0   32.1 
seston 

Total       28.8   20.0    0.0   31.3    4.8   34.8   27.4   15.3   38.3   31.8 
seston 


Organic  matter  content  of  the  surface  sediment, 
computed  in  the  model  as  (gBenthos+<?Detritus)  /^Surf  .sed.  , 
decreased  after  spillway  installation.   This  was  contrary 
to  in  situ  measurements,  which  showed  higher  organic  matter 
content  in  recent  sediments  (Chapter  2) ,  assuming  no 
decomposition  of  organic  matter  after  permanent 
sedimentation  (Gale  et  al.  1992).   Regardless  of  the 
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validity  of  this  assumption,  a  reduction  in  the  organic 
matter  content  of  the  surface  sediments  is  unlikely  when 
increased  water-column  SRP  concentration  probably  elevates 
primary  production.   The  reduction  was  produced  in  the 
model  by  a  44%  increase  in  the  inorganic  fraction  of 
surface  sediment  compared  with  15%  increases  in  benthos  and 
detritus  following  spillway  installation.   This  inorganic 
fraction  was  computed  in  the  model  by  dividing  the  ratio  of 
P  in  inorganic  sediment  (gp/m2)  and  the  P  content  of 
inorganic  sediment  (gP/ginorg.sed.)  bY  the  total  amount  of 
surface  sediments  (g/m2) .   Decreasing  the  value  of  the 
equilibrium  constant  (KR2)  in  the  sorption  reaction  between 
P  in  inorganic  sediment  and  P  in  sedimentary  pore-fluid 
from  3.06  to  2.67  (i.e.  equal  to  the  equilibrium  constant 
for  the  reaction  between  P  in  inorganic  seston  and  water- 
column  SRP)  eliminated  the  decrease  in  sedimentary  organic 
matter  postspillway  without  significantly  affecting  other 
model  output.   KR2  represented  the  ratio  at  equilibrium  of 
P  in  inorganic  sediments  and  SRP  in  interstitial  pore- 
fluid.   In  isotopic  exchange  studies,  using  32P  as  a 
tracer,  Li  et  al.  (1972)  found  values  for  this  equilibrium 
constant  ranging  from  1.33  to  4.26  for  four  Wisconsin 
lakes.   Thus,  decreasing  KR2  from  3.06  to  2.67  kept  this 
equilibrium  constant  within  a  plausible  value  range. 

Soluble  P  concentration  in  the  interstitial  pore-fluid 
doubled  after  spillway  installation  (Figure  5-7) .   No 
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implausible  shifts  in  the  composition  of  total  sedimentary 
phosphorus  occurred  (Table  5-3) .   The  ratio  of  average  P 
concentration  in  inorganic  sediment,  benthos,  detrital 
sediment,  and  pore-fluid  was  3.2  :  2.2  :  4.5  :  1 
prespillway  and  3.1  :  1.7  :  3.5  :  1  after  spillway 
installation. 

In  order  to  compare  simulated  total  sedimentary  P 
concentrations  with  those  measured  in  the  cores,  the  former 
needed  to  be  converted  to  gP/gSurf  .sed.  •   Using  the  average 
simulation  values  for  pre-  and  postspillway  surface 
sediment  (1293  and  1602  gSurf .sed./m2,  respectively),  both 
pre-  and  postspillway  sedimentary  P  concentration  equaled 
3.1  mgp/gSurf  8ed  .   These  values  were  identical  to  those 
measured  in  the  center  of  the  lake  (Chapter  2). 


Table  5-3.   Pre-  and  postspillway  concentrations  of 
sedimentary  P  in  the  inorganic,  benthos,  detrital,  and 
soluble  form.   Data  from  the  standard  model  run.   Values 
reported  are  the  average,  minimum,  and  maximum  (in  g/m2) ; 
time-of -minimum  and  time-of -maximum  (in  years  after  model 
start) . 


Prespillway  Postspillway 

Avg   Min  Tmin   Max  Tmax   Avg   Min  Tmin   Max  Tmax 


Inorganic    1.17  0.96  15.0   1.46   0.5   1.68  0.97  15.0   2.06  29.5 
sed.  P 

Benthos      0.80   0.40   0.0    0.85   5.9    0.92   0.80  16.4    1.08  30.8 
sed.  P 

Detrital      1.65   1.37   0.1    1.73   6.0    1.90   1.64  15.6    2.20  30.9 
sed.  P 

Pore-fluid    0.37   0.31  15.0    0.47   0.5    0.54   0.31  15.0    0.66  29.5 
sed.  P 

Total        4.00   3.76   0.0    4.17   3.5    5.04   3.81  15.0    5.61  29.5 
sed.  P 
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Postspillway  in  situ  sedimentation  rates  were  45,  19, 
and  10%  higher  at  increasing  distances  from  the  dam 
(Chapter  2).   In  the  model,  sedimentation  rate  increased 
27%  from  647  to  823  g  m"2y_1.   Hence,  in  situ  and  simulated 
rates  compared  reasonably  well. 

Both  the  similarity  between  measured  and  modeled 
changes  for  sedimentary  P  and  for  sedimentation  rate,  as 
well  as  the  plausible  behavior  of  the  model  for  water- 
column  P,  pore-fluid  P,  and  seston  add  support  to  the 
feedback  loop  hypothesis.   The  decrease  in  organic  matter 
content  of  surface  sediments  in  the  simulation  does  not 
correspond  to  field  data.   This  may,  in  part,  be  corrected 
by  adjusting  the  estimate  for  the  equilibrium  constant  in 
the  sorption  reaction  between  P  in  pore-fluid  and  P  in 
inorganic  sediment. 
Sensitivity  Analysis 

The  objectives  of  the  sensitivity  analysis  were:  (1) 
to  enhance  understanding  of  model  behavior,  (2)  to  test  the 
model  for  sensitivity  to  errors  in  parameter  estimates,  and 
(3)  to  identify  sensitive  parameters  which  require  further 
field  and/or  laboratory  study  to  increase  confidence  in  the 
model . 

All  sixty-eight  model  parameters  were  individually 
halved  and  doubled.   Parameters  involved  with  the 
fractionation  of  plankton,  benthos,  and  littoral  material 
were  changed  by  increasing  and  decreasing  the  largest 
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fraction  by  20%,  while  changing  the  remaining  two  fractions 
proportionally  to  maintain  a  total  of  100%.   The  effect  of 
all  parameter  changes  on  water  depth,  seston,  water-column 
SRP,  organic  matter  content  of  surface  sediments,  and 
sedimentation  rate  is  given  in  Table  5-4.   Output 
characteristics  of  interest  included  for  both  pre-  and 
postspillway  conditions:  (1)  changed  averages,  minima,  and 
maxima  as  percent  of  standard  run  output,  and  (2)  altered 
time-of-minimum  and  time-of -maximum  as  percent  of  total 
simulation  time  (40  years) .   The  parameters  in  Table  5-4 
were  grouped  per  output  variable  and  were  ranked  in 
descending  order  of  activity  relative  to  the  percent  change 
in  the  average  value  of  the  standard  run  output.   Parameter 
activity  is  defined  as  the  absolute  value  of  the  percent 
change  in  an  output  variable  when  the  parameter  is  halved, 
plus  that  when  the  parameter  is  doubled.   A  "+•«  sign 
indicates  that  doubling  the  parameter  increases  and  halving 
decreases  the  value  of  the  output  variable.   A  "-"  sign 
indicates  the  opposite.   Values  in  italics  indicate  that 
both  changes  of  the  parameter  cause  the  output  to  move  in 
the  same  direction.   This  occurs  when  the  initial  parameter 
value  was  optimal.   Only  changes  greater  than  25%  are 
reported . 

Water  depth  was  relatively  insensitive  to  halving  and 
doubling  of  individual  parameter  values  (Table  5-4) .  Only 
a  few  parameters  produced  changes  greater  than  25%  in  water 
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Table  5-4.   Results  of  sensitivity  analysis.   Values 
reported  under  average,  minimum,  and  maximum  are  the  sum 
of  the  absolute  values  of  the  percentage  changes  in  the 
output-variable  (compared  to  standard  run  output)  when 
the  parameter  is  doubled  and  halved.  Values  reported 
under  time-of -minimum  and  time-of -maximum  represent  the 
sum  of  the  absolute  values  of  the  change  compared  to 
standard  run  output  (as  percentage  of  total  simulation 
time)  when  the  parameter  is  doubled  and  halved.  A  +  sign 
indicates  that  doubling  the  parameter  increases,  and 
halving  decreases  the  value  of  the  output  variable.   A  - 
sign  indicates  the  opposite.   Values  in  italics  indicate 
that  both  parameter  changes  cause  the  output  to  move  in 
the  same  direction.   Only  changes  >25%  are  reported. 
See  Table  5-1  for  parameter  description. 


OUTPUT  VARIABLE 

DEPTH  (m) 

PRE-SPILLWAY 

POST-SPILLWAY 

Parameter 

Avg 

Min 

Tmin 

Max 

Tmax 

Avg 

Min 

Tmin 

Max 

Tmax 

Output  Standard  Run 

1.50 

0.94 

14.7 

2.06 

0.30 

1.59 

1.17 

39.7 

2.01 

16.2 

AMPNODAM 

AMPDAM 

MAXBENTHICBACT 

BENTHICBACTFRAC 

PBENTHICBACTCONT 

PBENTHICEUKCONT 

Q 

SUBMLITAREA 

N 

BENTHOSSPECREGRT 

-85 

-32 

+35 
+35 
-32 

+41 

+27 
+27 

-47 
+31 
+32 
+32 
-32 

-32 

+29 

Table   5-4--continued. 
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OUTPUT  VARIABLE 

SESTON  (g  m-3) 

PRE- 

SPILLWAY 

POST-SPILLWAY 

Parameter 

Avg 

Min 

Tmin 

Max 

Tmax 

Avg 

Min 

Tmin 

Max 

Tmax 

Output  Standard  Run 

+28.8 

20.0 

0.0 

31.3 

4.7 

34.8 

27.4 

15.3 

38.3 

31.8 

ALPHA 

+  174 

+43 

+192 

+37 

+264 

+230 

+292 

BENTHOSSPECREGRT 

+32 

90 

+35 

♦155 

+103 

+203 

♦33 

PBENTHICEUKCONT 

+32 

90 

+35 

+155 

+103 

+203 

+33 

PBENTHICBACTCONT 

-32 

89 

-35 

-153 

-102 

-201 

-33 

BENTHICBACTFRAC 

-32 

88 

-35 

-153 

-102 

-202 

-33 

MAXBENTHICBACT 

-31 

82 

-35 

-147 

-97 

-193 

-33 

SUBMLITAREA 

+116 

+29 

-36 

+109 

+27 

+134 

+  138 

+132 

SVELP 

-54 

+36 

-78 

-33 

-132 

-98 

-164 

-33 

Q 

+91 

-36 

+92 

+35 

+130 

+124 

+131 

PLANKSPECREGRT 

63 

+33 

+114 

♦74 

+152 

♦33 

PEUKCONT 

62 

+33 

+  108 

♦70 

+144 

+33 

PBACTCONT 

-47 

-36 

-62 

-33 

-108 

-83 

-135 

-33 

SDR 

-84 

-86 

-37 

-101 

-101 

-103 

PLUPLANKTONCONT 

+52 

+53 

+66 

+61 

+68 

SVELD 

-56 

-53 

-59 

-62 

-58 

PDETRITUSCONT 

-57 

-60 

-56 

-56 

-61 

-55 

BACTFRAC 

26 

-33 

-55 

-40 

-70 

-33 

PAUTOCONT 

-25 

-26 

-33 

-51 

-42 

-59 

-33 

AMPNODAM 

-47 

-47 

+33 

-71 

N 

-46 

-39 

+33 

-66 

PINORGCONT 

-28 

-30 

-26 

-36 

-27 

-38 

PDELLITPLANKTONCON 

+33 

+35 

♦29 

+39 

+26 

AMPDAM 

-27 

+30 

42 

KT 

-27 

-26 

+25 

SVELI 

-25 

-25 

M 

+27 

-27 

BLIFETM 

-30 

KDETSED 

-30 

BULKDENS 

-27 

KM1 

+27 

KR2 

+27 

Table   5-4--continued, 
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OUTPUT  VARIABLE 

SRP-WATER  (g  m-3) 

PRE- 

SPILLWAY 

POST-SPILLWAY 

Parameter 

Avg 

Min 

Tmin 

Max 

Tmax 

Avg 

Min 

Tmin 

Max 

Tmax 

Output  Standard  Run 

0.016 

0.006 

0.0 

0.022 

0.7 

0.023 

0.012 

15.2 

0.028 

31.7 

ALPHA 

♦764 

+69 

+945 

+37 

+1056 

+  1173 

+1255 

BENTHOSSPECREGRT 

♦566 

+58 

-37 

+794 

♦37 

+1010 

+994 

+1405 

+33 

PBENTHICEUKCONT 

+565 

+58 

-37 

793 

+37 

+1008 

+992 

+1402 

+33 

BENTHICBACTFRAC 

-561 

-61 

+37 

-784 

-37 

-1000 

-985 

1389 

-33 

PBENTHICBACTCONT 

-560 

-58 

+37 

-784 

-37 

-999 

-983 

-1388 

-33 

MAXBENTHICBACT 

-534 

-61 

+38 

-742 

-37 

-951 

-938 

-1312 

-33 

PLANKSPECREGRT 

+430 

+52 

-38 

+583 

+37 

+757 

+755 

+1010 

+33 

PEUKCONT 

+413 

+51 

-38 

+557 

+37 

+726 

+725 

+964 

+33 

SVELP 

-367 

-48 

+38 

-482 

-37 

-634 

-637 

-828 

-33 

PBACTCONT 

-301 

-42 

+35 

-379 

-37 

-519 

-514 

-673 

-33 

BACTFRAC 

-197 

+35 

-218 

-35 

-320 

-320 

-399 

-33 

PAUTOCONT 

-151 

+35 

-142 

-35 

-221 

-228 

-258 

-33 

KR1 

-106 

-136 

-89 

-103 

-85 

SUBMLrTAREA 

+88 

+57 

+99 

+  116 

+99 

-30 

KM1 

♦79 

+50 

+35 

+98 

+119 

+106 

Q 

+69 

+44 

+35 

+96 

♦  105 

+99 

BULKDENS 

-65 

-37 

-33 

-72 

-96 

-76 

KDETSED 

-54 

-32 

-35 

-69 

-82 

-83 

-30 

BLIFETM 

-54 

-33 

-35 

-69 

-82 

-83 

-28 

PLIFETM 

+45 

+63 

+59 

+71 

BEDDENS 

+36 

+58 

+57 

+63 

+28 

V1 

-55 

-81 

-43 

-45 

-49 

MAXBENTHICAUTO 

-40 

-53 

-53 

-59 

PLITINORGCONT 

+35 

+52 

+48 

+58 

SVELI 

-50 

-54 

-51 

-57 

-43 

KR2 

+31 

36 

+35 

+51 

+72 

+51 

RISEDAM 

+50 

+64 

AMPNODAM 

-50 

+88 

-25 

-104 

N 

-45 

+35 

-30 

-78 

PDETSEDCONT 

-40 

-34 

-58 

-34 

SVELD 

-38 

-34 

-44 

MAXPAR 

-27 

-37 

-35 

-42 

MAXPLANKAUTO 

-27 

-37 

-35 

-42 

OPTPAR 

+27 

+37 

+35 

+42 

SDR 

26 

62 

36 

56 

48 

AMPDAM 

-35 

-30 

PLITPLANKTONCONT 

+26 

+30 

PDELLITINORGCONT 

+30 

PINORGSEDCONT 

-28 



Table   5-4--continued. 
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OUTPUT  VARIABLE 


SED.ORG.  MATTER  (%) 


Parameter 


PRE-SPILLWAY 


Avg 


Min 


Tmin 


Max 


Tmax 


POST-SPILLWAY 


Avg 


Min 


Tmin     Max     Tmax 


Output  Standard  Run 


0.70 


0.61 


0.5 


0.74 


15.0 


0.64 


0.61 


31.5 


0.74 


BENTHICBACTFRAC 

PBENTHICEUKCONT 

BENTHOSSPECREGRT 

MAXBENTHICBACT 

PBENTHICBACTCONT 

PLANKSPECREGRT 

PEUKCONT 

ALPHA 

SVELP 

PBACTCONT 

BACTFRAC 

PAUTOCONT 

PINORGSEDCONT 

SDR 

KR2 

KM1 

BULKDENS 

PDETSEDCONT 

PLIFETM 

KDETSED 

BLIFETM 

BEDDENS 

SUBMLITAREA 

Q 

PLriPLANKTONCONT 

AMPNODAM 

MAXBENTHICAUTO 

AMPDAM 

RISEDAM 


+131 

-131 

-131 

+129 

+130 

-108 

-106 

-104 

+101 

+88 

+65 

+54 

+49 

-37 

-37 

-31 

+26 


+  103 

-103 

-103 

+102 

+  101 

-90 

-88 

-85 

+83 

+72 

+51 

+35 

+54 


-45 


-37 
+37 
+37 
-37 
-37 
+36 
+36 
+36 
-36 
-36 
-36 
-36 

+25 
+34 
+33 
-25 


-35 
-35 


+44 
-44 
-44 
+44 
+43 
-39 
-39 
-39 
+39 
+37 
+33 
+33 
+45 

-28 
-27 
+25 


73 
73 
73 
73 
73 
68 
68 
70 
68 
68 
+38 
+38 

-37 
-37 
-37 
+37 


+35 
+35 


+152 

+159 

-151 

-158 

-151 

-158 

+150 

+159 

+150 

+158 

-136 

-143 

-134 

-141 

-131 

-138 

+129 

+137 

+121 

+  129 

+96 

+104 

+80 

+87 

+58 

+62 

-49 

-54 

-49 

-52 

-40 

-42 

+32 

+32 

-30 

-34 

-29 

-33 

-27 


+27 


-59 
+59 
+59 
-59 
-59 
+58 
+58 
+58 
-58 
-58 
-55 
-53 

-40 
+25 


-38 
-35 
+43 
-40 
-40 
-35 
+30 

-27 

+25 


+129 

-128 

-128 

+127 

+  128 

-114 

-112 

-109 

+108 

+97 

+74 

+60 

+45 

-41 

-41 

-37 

+31 


15.0 


66 

66 
66 
66 
66 


Table  5-4--continued. 
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OUTPUT  VARIABLE 

SEDIMENTATION  RATE  (g 

m-2y-1) 

PRE- 

SPILLWAY 

POST-SPILLWAY 

Parameter 

Avg 

Min 

Tmin 

Max 

Tmax 

Avg 

Min 

Tmin 

Max 

Tmax 

Output  Standard  Run 

647 

595 

0.1 

669 

4.4 

823 

629 

15.0 

879 

31.4 

PBENTHICBACTCONT 

-127 

-26 

+38 

-197 

-36 

-314 

-241 

-389 

-34 

BENTHOSSPECREGRT 

+124 

-38 

+198 

+36 

+314 

+240 

+390 

+34 

PBENTHICEUKCONT 

+124 

25 

-38 

+197 

+36 

+313 

+240 

+389 

+34 

BENTHICBACTFRAC 

-124 

25 

+38 

-195 

-37 

-312 

-239 

-387 

-34 

MAXBENTHICBACT 

-118 

-38 

-185 

-37 

-297 

-227 

-370 

-34 

PLANKSPECREGRT 

+93 

-37 

+147 

+38 

+239 

+  182 

+297 

+34 

PEUKCONT 

+89 

-37 

+141 

+38 

+230 

+175 

+285 

+34 

SVELP 

-77 

+37 

-120 

-34 

-202 

-152 

-250 

-34 

PBACTCONT 

-66 

+37 

-96 

-34 

-168 

-126 

-208 

-34 

ALPHA 

+84 

+100 

+36 

+  160 

+136 

+  184 

SUBMLITAREA 

+  109 

+45 

-37 

+83 

+26 

+124 

+133 

+122 

Q 

+85 

+38 

-37 

+74 

+36 

+121 

+  119 

+120 

BACTFRAC 

-44 

+37 

-57 

-33 

-106 

-80 

-129 

-32 

PDETSEDCONT 

-101 

-95 

+37 

-95 

-94 

-113 

-91 

SDR 

+75 

+104 

+37 

+  117 

-36 

♦73 

+64 

+77 

PAUTOCONT 

-34 

♦37 

-37 

-33 

-74 

-59 

-86 

-32 

PINORGSEDCONT 

-53 

-55 

-54 

-64 

-50 

-68 

PLrTPLANKTONCONT 

+46 

-37 

+35 

+57 

+54 

+57 

AMPNODAM 

-49 

-59 

+35 

-29 

-75 

KDETSED 

-39 

-27 

+37 

-36 

-46 

-49 

-49 

KR2 

+28 

-37 

+36 

+46 

♦40 

+49 

BLIFETM 

-38 

-26 

♦37 

-36 

-46 

-48 

-47 

N 

-45 

-53 

+35 

-26 

-58 

KM1 

-37 

+30 

+34 

♦32 

+36 

BEDDENS 

-37 

+31 

+33 

AMPDAM 

-30 

+35 

-30 

-35 

PDELLITPLANKTONCON 

+29 

-37 

+35 

PLmNORGCONT 

-37 

+28 

+30 

KT 

-26 

BULKOENS 

+37 

-28 

B3 

-37 

MAXBENTHICAUTO 

+37 

MAXPLANKAUTO 

+37 

MAXPAR 

+37 

PLIFETM 

-37 

BWDNODAM 

+37 

RISEDAM 

+32 

M 

-27 
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depth.   Halving  or  doubling  the  amplitude  of  water-level 
fluctuations  resulted,  of  course,  in  considerable  changes 
in  the  minimum  and  maximum  water  depths.   With  the 
exception  of  the  effect  of  ALPHA  on  model  behavior 
(discussed  below) ,  the  remaining  four  output  variables  were 
all  most-strongly  affected  by  changes  in  the  following 
parameters:  (l)  the  regeneration  of  P  by  eukaryotic 
heterotroph  benthos,  (2)  the  P  content  of  benthic  bacteria, 
(3)  the  P  content  of  eukaryotic  heterotrophs ,  (4)  the 
bacterial  fraction  of  benthos,  and  (5)  the  maximum  P-uptake 
rate  for  benthic  bacteria.   The  feedback  loops  in  which 
these  parameters  are  involved  occur  in  the  surface 
sediments  and  affect  the  rates  of  P  exchange  between 
bacteria  and  benthic  heterotrophs  with  the  interstitial 
pore-fluid  loops  21,  22,  and  28;  Figure  5-5).   Research 
focused  on  improving  the  estimates  of  these  active 
parameters  may  enhance  confidence  in  the  model.   The  rate 
of  P  uptake  by  benthic  bacteria,  in  particular,  has  no 
strong  foundation  in  the  literature.   Changes  in  parameters 
involved  in  P  exchanges  with  detritus  in  the  form  of  rates 
of  death  and  decomposition  have  a  much  smaller  impact  on 
model  behavior. 

ALPHA,  a  model  coefficient  determining  the  specific 
rate  of  resuspension  of  surface  sediments  (loops  7-10; 
Figure  5-3),  is  the  most  active  parameter  for  seston  and 
water-column  SRP  concentrations.   Its  effects  on  organic 
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matter  content  of  the  sediments  and  sedimentation  rate, 
however,  are  considerably  lower  than  the  above-discussed 
coefficients.   Similarly,  changes  in  the  settling 
velocities  of  inorganic  and  detrital  suspended  matter  are 
of  intermediate  significance  to  seston  and  water-column  SRP 
concentrations,  but  have  no  noticeable  influence  on 
sediment  organic  matter  content  or  sedimentation  rate. 
Therefore,  the  modeled  rate  of  sedimentation  (i.e.  the  rate 
at  which  surface  sediments  are  converted  to  permanent 
sediments)  is  not  as  much  a  function  of  the  net  difference 
between  settling  and  resuspension,  but  more  a  function  of 
benthic  processes  acting  on  sediments  once  they  have 
settled  to  the  bottom. 

To  test  whether  the  five  statistics  (used  in  Table  5- 
4)  were  adeguate  in  describing  changes  in  model  behavior 
due  to  halving  and  doubling  of  parameter  values,  dynamics 
for  the  two  most  active  parameters  were  analyzed.   Doubling 
ALPHA  (Figure  5-8)  and  halving  PBENTHICBACTCONT  (Figure  5- 
9)  showed  that  no  major  pattern  changes  (additional 
inflection  points  in  the  plots)  occurred  that  were 
different  than  those  from  the  standard  run  output  (Figure 
5-7) .   The  relative  stability  of  sedimentary  organic  matter 
content  in  the  standard  run  was  repeated  after  the  model 
output  passed  through  an  starting  transient  caused  by  a 
high  initial  value.   Pore-fluid  SRP  is  absent  from  these 
two  graphs  (compared  with  the  standard  run)  because  it  was 
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MODEL  OUTPUT  WHEN  ALPHA  IS  DOUBLED 


15         20         25 

Time  (years) 


Figure  5-8.   Model  output  when  ALPHA  is  doubled. 
Spillway  installation  initiated  at  15  years. 


MODEL  OUTPUT  WHEN  P-CONTENT  OF  BENTHIC  BACTERIA  IS  HALVED        14° 
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Figure  5-9.   Model  output  when  the  P  content  of 
benthic  bacteria  is  halved.   Spillway  installation 
initiated  at  15  years. 
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not  included  in  the  sensitivity  analysis.   Water-column 
depth,  which  was  part  of  the  analysis,  is  not  shown  because 
it  was  not  significantly  affected  by  the  two  parameters 
involved  in  this  test. 

The  plausibility  of  the  model  response  to  a  change  in 
a  parameter  may  be  analyzed  by  following  this  change 
through  the  equation  structure.   For  example,  doubling 
ALPHA  (Figure  5-8)  doubles  the  specific  rate  of 
resuspension.   This  results  in  an  increase  in  seston  and 
water-column  SRP  (hence  the  "+"  signs  for  the  influence  of 
ALPHA  on  these  two  model  variables  in  Table  5-4).   The 
magnitude  of  this  increase  depends  on  water  depth.   More 
seston,  in  turn,  produces  a  higher  rate  of  settling.   This 
rate,  however,  is  not  the  same  for  all  fractions  of  seston. 
Inorganic  seston  settles  faster  than  planktonic  seston  by 
an  order  of  magnitude.   This  produces  a  rapid  decline  in 
organic  matter  content  of  the  surface  sediments  (hence  the 
"-"  sign  for  the  influence  of  ALPHA  on  SEDORGMAT  in  Table 
5-4) .   Increased  settling  leads  to  more  surface  sediments. 
This  positive  resuspension-settling  loop  is  limited  by  the 
effect  of  negative  loops  between  resuspension  and  surface 
sediments  and  between  seston  and  settling  (see  also  Figure 
5-3) .   The  increase  in  water-column  SRP,  resulting  from 
doubling  of  ALPHA,  may  elevate  primary  production.   This 
produces  more  seston  and  increased  settling,  which  leads  to 
more  surface  sediments  and  higher  rates  of  sedimentation. 
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Again,  model  behavior  is  plausible  as  shown  by  the  "+•• 
signs  for  the  influence  of  ALPHA  on  sedimentation  rate 
(Table  5-4) . 

Shifts  in  the  time  of  peak  occurrences  were  similar 
among  changes  in  different  parameters.   An  example  may 
clarify  this.   While  the  values  of  prespillway  minima  were 
generally  stable,  their  occurrence  was  often  35  to  38% 
later  than  in  the  standard  run.   This  change  always 
occurred  when  either  halving  or  doubling  the  parameter 
value  changed  the  model  variable  from  gradually  increasing 
to  steady  state  (i.e.  time-of -minimum  is  zero)  to  gradually 
decreasing  to  steady  state  (i.e.  time-of -minimum  is  15 
years  or  37%  of  simulated  time  later) .   The  variability  of 
time-of-minimum  values  between  33  and  38%  was  caused  by  the 
difference  in  phase  of  the  sinusoid  time  pattern  at  t=15 
years.   This  likely  resulted  from  inherent  inaccuracies  of 
calculations  performed  on  a  digital  computer  which  are 
compounded  during  the  80,000  iterations  performed  in  one 
model  run. 

In  general,  the  parameters  associated  with  activity  in 
the  littoral  zone  were  of  intermediate  or  low  activity. 
Changes  in  the  efficiency  of  the  wetland  fringe  to  trap 
particulate  matter  and  P  from  runoff  from  the  watershed  had 
no  noticeable  effect  on  model  output.   Since  this  trapping 
efficiency  was  proportional  to  the  extent  of  the  wetland 
fringe,  the  latter  had  only  a  minor  impact  on  P  and  seston 
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loading  of  the  profundal  zone  in  this  model.   Although  this 
loading  increased  after  spillway  installation,  it  remained 
generally  less  than  1%  of  the  total  input  to  the  profundal 
water-column  dominated  by  internal  loading  through 
resuspension  of  the  surface  sediments  and  pore-fluid. 

The  effect  of  light  on  model  output  was  also  of  low 
activity.   A  doubling  of  MAXPAR  (maximum  available 
photosynthetically  active  radiation  at  average  lake  depth) 
doubles  PAR.   The  effect  of  this  on  P-uptake  rate  is, 
however,  reduced  by  the  way  in  which  PAR  and  its  effect  on 
specific  P-uptake  rate  (PAREFF)  relate  (parabolic,  concave 
to  the  x-axis;  see  Figure  5-4b)  and  by  the  manner  in  which 
PAREFF  affects  the  specific  P-uptake  rate  (as  the  square 
root  of  the  product  of  PAREFF  and  PEFF,  the  effect  of  SRP 
on  the  specific  P-uptake  rate) . 

Individual  changes  in  parameters  involved  in  the 
effect  of  the  dam  on  amplitude  of  water-level  fluctuation 
(AMP) ,  concentration  of  seston  in  the  outflow  (BWD) ,  and 
hydraulic  residence  time  (HYDRESTM)  were  also  of  low 
activity. 
Analysis  of  Alternative  Management  Strategies 

One  of  the  model  objectives  was  to  evaluate  the 
consequences  of  alternative  management  actions  in  Newnan's 
Lake.   The  dam  may  accelerate  sediment  and  detritus 
accumulation  by  altering  biological  processes  (primary 
production) ,  changing  sediment  dynamics  (settling  and 


144 
resuspension) ,  and  by  hydrological  modifications  (increased 
water  residence  time)  (see  Figure  5-1) .   Each  of  these  can 
be  evaluated  separately  with  the  model  by  manipulating  the 
set  of  parameters  involved  in  each  pathway.   By  ranking  the 
relative  significance  of  each  pathway  on  model  output, 
risks  and  benefits  of  management  strategies  to  intended 
uses  of  the  lake  may  be  compared.   The  differences  in  model 
output  (compared  to  the  standard  run)  for  these  three 
management  scenarios  are  given  in  Table  5-5. 


Table  5-5.   Change  in  model  output  variables  (as  percent  of 
standard  run  output)  for  simulations  addressing  each 
modeled  impact  of  the  dam  separately  and  comparing  this 
with  results  from  a  simulation  with  no  dam  installed.   See 
text  for  a  description  of  the  simulations.   The  statistic 
reported  is  the  average  value  of  the  model  output  variable 
postspillway. 


Output     Seston 
variable 

Water 
Depth 

Water- 
Column  SRP 

Sed . org . 

matter 

Sediment . 
rate 

Littoral     -6.3 

area 
Resus-        1.2 

pension 
Water  resi-  -10.4 

dence  time 

0.3 
-0.1 
-7.4 

-8.5 

2.4 

-31.4 

0.1 
-1.6 
13.2 

-6.3 

1.1 

-16.1 

L  Impacts    -15.5 

Model  run   -19.2 
without  dam 

-7.2 
-7.0 

-37.5 
-39.2 

11.7 

14.5 

-21.3 
-27.2 

In  the  "littoral  area"  run,  the  wetland  fringe  around 
the  lake  was  kept  at  20%  of  lake  surface  area,  while  other 
impacts  of  the  spillway  were  maintained.   Such  a  scenario 
may  be  accomplished  by  hand-planting  of  vegetation  in  the 
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littoral  area  of  the  lake.   This  labor-intensive  program 
has  been  applied  to  large  areas  of  the  Newnan's  Lake 
shoreline  since  1990  (Krummrich,  Florida  Game  and  Fresh 
Water  Fish  Commission,  pers.  comm.).   Its  effect  on  model 
output  is  minor,  although  it  may  assist  in  creating  refuge, 
spawning,  nursery,  and  feeding  grounds  for  many  species  of 
fish  (Durant  1980;  Schramm  et  al.  1983;  Conrow  1984). 

In  the  "resuspension"  run,  the  rate  of  resuspension  of 
surface  sediments  after  spillway  installation  was 
maintained  at  prespillway  levels.   Its  effect  on  model 
output  was  also  minor  (Table  5-5) .   In  other  words,  the 
reduction  in  sediment  resuspension  produced  by  water-level 
stabilization  and  the  effect  of  this  on  sediment  flushing 
was  of  minor  importance  to  the  rate  of  sediment  and 
detritus  accumulation  in  the  lake.   In  a  separate  run,  the 
model  also  predicted  that  modifying  the  dam  to  allow 
bottom-water  drainage  will  have  little  effect.   Maintaining 
the  bottom-water  drainage  at  prespillway  levels  reduced 
sedimentation  rate  by  only  3%. 

The  impact  of  the  dam  on  lake  hydrology  is  of  greater 
significance  to  the  model  behavior  (Table  5-5) . 
Maintaining  the  water  residence  time  at  prespillway  levels 
and  eliminating  the  rise  in  water  level  after  spillway 
installation  reduced  sedimentation  rate  by  16%  and  water- 
column  SRP  by  31%.   Management  strategies  designed  to 
reduce  water  residence  time,  such  as  gravity  drawdowns, 
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appear  most  effective  in  reducing  sedimentation  rate.   A 
more  extreme  water-level  reduction,  such  as  may  be 
accomplished  under  a  mechanical  pumpdown,  was  simulated  by 
lowering  the  water  level  to  the  historic  low  of  19.65  m 
(msl)  for  a  period  of  6  months.   Shallow  water  depth  caused 
increased  resuspension  of  bottom  material,  which  produced 
15-20%  higher  levels  of  seston  and  water-column  SRP  during 
the  period  of  pumpdown.   Increased  resuspension  also 
lowered  the  amount  of  surface  sediments  by  8%  during  3-4 
years . 

The  combined  effects  of  dam  installation  on  model 
output  are  illustrated  by  a  model  run  without  the  dam. 
Without  the  dam,  the  average  rate  of  sedimentation  was  27% 
lower,  and  water-column  SRP  and  seston  were  reduced  by  39 
and  19%,  respectively  (Table  5-5).   It  is  interesting  to 
note  the  difference  between  the  sum  of  the  separate  effects 
of  the  spillway  installation  and  the  combined  effects  in 
the  model  run  without  the  dam  (Table  5-5) .   Some  non- 
linearity  appears  to  exist  resulting  from  an  interaction 
among  the  ways  in  which  the  dam  affects  sediment  and  P 
dynamics  in  the  lake.   Such  an  interaction  occurs,  for 
instance,  between  the  ways  the  dam  may  influence 
phytoplankton  production  (see  Figure  5-1) .   After  spillway 
installation,  production  may  be  enhanced  by  both  increased 
light  availability  due  to  reduced  sediment  resuspension  and 
by  increased  external  nutrient  loading  caused  by  a 
reduction  in  littoral  vegetated  fringe. 
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The  spillway  was  built  in  1967  with  the  management 
intent  of  raising  the  lake  stage  and  eliminating  extreme 
low-water  levels  (Smith  and  McGriff  1958) .   Hydrologic 
records  (Adkins  1991)  indicated  that  its  benefit  to  lake 
access  and  use  for  outboard,  recreational  boats  has  been 
minor.   The  model  demonstrated  a  7%  increase  in  average 
water  depth  postspillway.  Much  of  this  increase  has  been 
eliminated  after  25  years  by  higher  rates  of  sedimentation. 
Increased  seston  concentrations  may  decrease  the  aesthetic 
value  of  the  lake  and  oppose  lake  management  goals  to 
promote  the  lake  for  water  recreation.   A  39%  increase  in 
simulated  water-column  SRP  may  increase  plant  and/or  algal 
biomass  and  further  interfere  with  recreation  and 
navigation.   While  the  extent  of  the  vegetated  littoral 
zone  and  its  ability  to  filter  particulate  matter  and  P 
from  runoff  from  the  surrounding  watershed  likely  has  been 
reduced  by  water-level  stabilization,  this  reduction  has 
not  impacted  sediment  accumulation  rates  in  the  model. 
Short-term  removal  of  the  dam  produced  plausible 
results  in  the  model  (Figure  5-10) .   Conditions  in  this 
simulation  were  set  to  mimic  the  short-term  drawdown  and 
drought  of  the  spring  of  1989,  documented  in  Chapters  3  and 
4.   This  was  done  by  simulating  removal  of  the  spillway  at 
t=30  years  and  lowering  the  water  level  to  19.85  m  msl  for 
3  months.   Low  water  levels  increased  resuspension  of 
surface  sediments  and  pore-fluid  SRP.   More  suspended 
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Figure  5-10.   Model  output  with  spillway  installation 
at  15  years  and  removal  of  the  spillway  at  3  0  years 
for  a  period  of  3  months. 
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matter  in  the  water  column  depleted  water-column  SRP,  and 
produced,  after  a  settling  delay,  a  higher  organic  matter 
content  in  the  surface  sediments.   While  the  plausibility 
of  this  generated  time  pattern  provides  support  for  the 
feedback  hypothesis,  a  more  rigorous  validation  is 
desirable.   An  example  of  such  a  test  is  to  measure  seston 
or  water-column  SRP  in  the  field  before,  during,  and  after 
the  short-term  drawdown  and  compare  these  measurements  with 
the  model  output. 

The  current  long-term  dam  removal  (initiated  in  1991 
based  on  the  results  of  my  field  studies)  provides  future 
opportunities  for  model  validation.   Monthly  water  quality 
sampling  is  now  carried  out  in  Newnan's  Lake  (Ware,  St 
Johns  River  Water  Management  District,  pers.  comm.)  and  may 
be  helpful  in  this  respect.   The  effect  of  long-term 
spillway  removal  on  model  behavior  is  shown  in  Figure  (5- 
11) .   Removal  was  simulated  by  resetting  the  parameters 
directly  involved  with  the  dam  to  prespillway  conditions  at 
t=30  years.   The  model  output  variables  showed  a  gradual 
return  to  prespillway  levels  over  a  time  period  in  excess 
of  10  years. 

Summary  and  Conclusions 
The  model  has  synthesized  the  best,  current 
understanding  of  the  impact  of  water-level  stabilization  in 


Newnan's  Lake.   This  was  done  in  the  form  of  a  feedback 
loop  hypothesis  based  on  field  data,  literature,  opinion, 
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Figure  5-11.  Model  output  with  spillway  installation 
at  15  years  and  removal  of  the  spillway  at  3  0  years. 
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and  estimates.  Since  the  general  behavior  of  the  model  was 
quite  similar  to  available  real-system  dynamics,  and  since 
the  simple  response  pattern  of  the  model  was  not  altered  by 
systematic  parameter  changes,  support  has  been  gained  for 
this  feedback  loop  hypothesis  as  an  adequate  explanation  of 
sediment  and  P  dynamics  in  this  aquatic  ecosystem. 

The  model  is  most  sensitive  to  changes  in  feedback 
loops  involved  in  sediment  resuspension  and  in  actions  of 
bacteria  and  benthic  heterotrophs  on  the  production  and 
movement  of  P  in  the  surface  sediments.   Hence,  parameters 
affecting  P  flows  between  benthos  and  sedimentary  pore- 
fluid  and  those  affecting  wind-wave  induced  resuspension 
require  additional  investigation  to  enhance  confidence  in 
the  model.   Field  measurements  of  these  processes  are 
scarce.   Future  work  should  focus  on  (1)  determining  the 
flow  of  (labeled)  P  in  surface  sediments  and  (2) 
quantifying  the  relationship  between  wind -wave  action, 
water  depth  and  resuspension  of  bottom  material. 

Each  effect  of  the  spillway  separately  had  minor 
implications  for  the  behavior  of  the  model.   In 
combination,  however,  they  resulted  in  elevated  seston, 
water-column  SRP  and  sedimentation  rate.   Over  a  period  of 
25  years,  the  intended  benefit  of  the  dam  (lake  access, 
navigation)  has  been  eliminated  and  replaced  with  common 
lake  management  "costs"  such  as  increased  water-column 
turbidity,  primary  production  and  material  deposition  in 
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the  lake.   Spillway  removal,  recently  initiated,  provides 
an  excellent  opportunity  to  test  further  whether  the  major 
causal  forces  between  the  spillway  and  sediment-P  dynamics 
in  Newnan's  Lake  have  been  identified  in  this  model. 


CHAPTER  6 
SUMMARY  AND  DISCUSSION 

The  long-term  significance  of  a  spillway  in  the  outlet 
of  Newnan's  Lake  for  the  transfer  of  bulk  sediment,  organic 
matter,  and  nutrients  between  water  and  sediment  was 
evaluated  using  a  paleolimnological  approach.   The  effect 
of  a  short-term  removal  of  the  spillway  on  the  transfer  of 
this  material  out  of  the  lake  (e.g.  flushing,  oxidation)  or 
within  the  system  (e.g.  erosion,  redistribution)  was 
determined  in  a  number  of  field  and  laboratory  experiments. 
The  findings  of  these  studies  were  integrated  with 
independently  derived  information  from  literature  and 
theory  in  a  feedback  dynamics  model,  and  their  implications 
for  lake  management  were  tested  in  a  series  of  simulations. 

Paleolimnological  analysis  of  the  sedimentary  record 
in  Newnan's  Lake  allowed  interpretations  about  past 
conditions  of  the  ecosystem.  The  sedimentary  signal  at  the 
time  of  spillway  installation  and  the  agreement  between 
210Pb  and  137Cs  derived  dates  were  remarkable,  particularly 
since  the  shallow  depth  of  the  lake  and  the  flocculent 
nature  of  the  sediments  likely  produced  a  more  disturbed 
stratigraphy  than  that  found  in  the  traditional 
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paleolimnology  of  deeper  lakes  (cf.  Shapiro  et  al.  1971; 
Robbins  and  Edgington  1975;  Pilskaln  and  Johnson  1991). 
This  corresponds  with  earlier,  promising  paleolimnological 
data  from  other  shallow  productive  lakes  in  Florida 
(Brenner  and  Binford  1986;  Deevey  et  al.  1986)  and  is 
encouraging  for  paleolimnological  work  in  such  systems. 

The  finding  in  this  study  that  net  accumulation  rate 
of  bulk  sediment  increased  by  45%  since  spillway 
construction,  while  rates  for  organic  matter  and  nutrients 
tripled  compared  to  prespillway  conditions,  has 
implications  for  lake  depth  and  productivity.   More 
sediment  deposition  increases  the  rate  of  filling  of  the 
lake.  A  higher  nutrient  content  of  the  surficial  sediments 
will  enhance  internal  loading  of  nutrients  through  mixing 
of  this  material  in  the  water  column.   As  such,  the 
spillway  may  accelerate  the  successional  development  of  the 
lake  towards  a  freshwater  marsh.   Further  paleolimnological 
work  may  supply  additional  lines  of  evidence  to  support 
this  conclusion.   For  instance,  changes  in  lake 
productivity  may  be  reflected  in  relative  abundance  of 
subfossil  Chironomus  sp.  and  Tanvtarsus  sp.  in  the  sediment 
profile.   These  highly  ubiguitous  midges  leave  identifiable 
morphological  remains  in  lake  sediments  and  Chironomus  sp. 
is  generally  considered  to  be  indicative  of  a  higher 
trophic  state  (Crisman  1978) .  Alternatively,  the  relative 
abundance  of  chydorid  Cladocera  remains  in  the  sediments 
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may  illustrate  historical  changes  in  the  lake's  littoral 
productivity.   Chydoridae  are  generally  non-planktonic  and 
occur  in  association  with  littoral  and  benthic  substrate 
(Goulden  1969) . 

The  short-term  drawdown  demonstrated  the  improbability 
of  lowering  the  water  level  in  the  lake  by  a  gravity 
drawdown  to  a  stage  less  than  that  accomplished  (19.70  m 
MSL) .   Small  elevation  gradients  in  the  lake  basin  and  the 
gradual  build-up  of  obstructions  in  the  lake  outlet 
upstream  and  downstream  from  the  spillway  since  its 
construction  (Crider  1972;  personal  observation)  likely 
prevented  a  more  dramatic  drawdown  of  water  level. 
However,  dredging  of  the  lake  outlet  to  enhance  outflow 
conditions  and  increase  the  amplitude  of  a  drawdown  may  not 
be  necessary.   Periodic  removal  of  the  spillway  will  likely 
scour  the  outflow  to  accommodate  larger  discharges  and 
return  this  channel  to  its  historic  configuration.   During 
this  drawdown,  the  water  level  was  barely  lakeward  of  the 
cypress  tree  fringe  and  drying  and  consolidation  of  lake 
bottom  was  limited  to  a  narrow  stretch  of  littoral  zone.  In 
other  Florida  lakes,  significant  consolidation  only 
occurred  when  substrate  was  exposed  for  a  minimum  of 
several  months  and  water  levels  were  significantly  below 
the  level  of  the  unconsolidated  sediments  (Greening  and 
Doyon  1990) .   A  shorter  exposure  time  may  result  in  the 
formation  of  a  "cap"  of  dried  organic  material  over 
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unconsolidated  wet  sediments.   This  cap  may  then  become 
dislodged  and  float  to  the  water  surface  during  and  after 
refill  of  the  lake. 

Flushing  rates  of  particulate  organic  matter  and 
nutrients  from  the  lake  during  the  drawdown  were  low.   This 
removal,  however,  was  accomplished  at  low  cost  since 
drawdowns  are  unquestionably  among  the  least  expensive  lake 
management  techniques  (Cooke  et  al.  1986).   Drawdowns  may 
be  an  appropriate  management  strategy  in  a  shallow,  exposed 
lake  such  as  Newnan's  Lake,  where  downstream  areas  are  able 
to  absorb  particulate  matter  loadings  during  the  period  of 
flushing. 

Maximum  flushing  of  organic  matter  and  nutrients  from 
Newnan's  Lake  by  temporary  removal  of  the  dam  occurred  at  a 
high  lake  stage  and  during  periods  of  high  average  wind 
speed.  Maximum  monthly  lake  stage  in  Newnan's  Lake  occurs 
either  from  February  to  April  or  from  July  through 
September  (Adkins  1991) .   No  wind  data  were  analyzed  for 
north-central  Florida,  but  data  for  Lake  Okeechobee 
indicated  that  maximum  wind  velocities  between  October  and 
May  were  approximately  30  percent  greater  than  between  June 
and  September  (Maceina  and  Soballe  1990) .   Assuming  that 
these  data  can  be  applied  to  Newnan's  Lake,  it  appears  that 
short-term  drawdowns  during  February-April  have  the  highest 
probability  of  leading  to  improved  lake  conditions. 
Drawdowns  during  this  time  period,  however,  are  not  optimal 
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from  a  fisheries  perspective.   Fish  generally  spawn  in 
February  in  Florida  and  need  nursery  areas  in  spring  and 
early  summer.   Frequent  drawdowns  (i.e.  every  year)  during 
February-April  are  therefore  not  recommended  if  the  lake  is 
managed  for  sportfishing.   Spillway  closure  in  May  will 
allow  re-filling  of  the  lake  by  convective  summer  rains. 
In  order  to  more  firmly  establish  a  recommended  time  span 
for  board  removal,  research  should  focus  on  (1)  analysis  of 
average  hourly  wind  speed  throughout  the  year  for  the 
Newnan's  Lake  watershed,  and  (2)  hydrologic  modeling  to 
evaluate  the  refill  potential  of  the  lake  following  a 
February-April  drawdown  compared  to  a  July-September 
drawdown . 

Erosion  of  flocculent  littoral  substrate  was  recorded 
during  the  drawdown.  Although  variability  was  high,  an 
average  of  6.08  g  m"2day_1  of  organic  matter  (0.58 
gTKN  m"2day_1  and  0.03  gTp  m"2day"1)  was  removed  from  the 
littoral  zone.   Sediment  depth  to  sand  decreased  by  an 
average  of  42%  after  drawdown,  and  the  bulk  density  of  the 
remaining  substrate  increased  by  an  average  of  250%. 
Removal  of  soft  substrate  from  the  littoral  zone  may 
promote  rooted  plant  growth  which  can  serve  as  nursery  area 
and  refuge  for  fish.   In  addition,  such  removal  may  improve 
littoral  habitat  for  future  sportfish  spawning. 

The  quantitative  record  of  the  transfer  of  this 
material  to  deeper  areas  of  the  lake  was  not  clear.   The 
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small  drop  in  lake  level  during  the  drawdown  did  not  create 
large  areas  of  erosion  in  the  littoral  zone.   This  reduced 
the  potential  redistribution  of  material  in  the  lake.   In 
addition,  the  use  of  field  markers  would  have  eliminated 
error  due  to  imprecise  matching  of  pre-  and  postdrawdown 
stations.   These  markers  (e.g.  buoys,  flags)  were  not  used 
in  the  study  because  of  safety  concerns  for  people  using 
the  lake  for  recreation  and  sportfishing. 

A  sediment  and  phosphorus  model  was  developed  to 
synthesize  the  best  current  understanding  of  the  effects  of 
the  spillway  on  Newnan's  Lake.   It  combined  the  results 
from  laboratory  and  field  study  with  information  from 
literature,  opinion,  and,  where  necessary,  guesses. 
Conclusions  from  the  model  were  given  in  the  previous 
chapter.   Here,  the  emphasis  is  on  an  evaluation  of  the 
interaction  between  the  experimental  components  of  this 
study  and  the  model.   Results  from  field  and  laboratory 
work  assisted  in  model  development  and  analysis,  especially 
in  quantifying  feedback  relationships,  evaluating  model 
behavior  during  the  sensitivity  analysis,  and  field  testing 
the  model.   Conversely,  the  model  was  used  to  improve  the 
understanding  of  the  system  and  to  identify  specific 
research  areas  needing  further  study.   In  this  case,  model 
analysis  suggested  research  needs  in  quantifying  the 
relationship  between  wind-wave  action,  water  depth,  and 
resuspension  of  bottom  material,  as  well  as  measuring  the 
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flow  of  phosphorus  in  surface  sediments  between  benthos  and 
interstitial  pore-fluid.   The  actions  of  bacteria  and 
benthic  heterotrophs  showed  a  considerable  impact  on  the 
production  and  movement  of  phosphorus  and,  thus,  on  the 
productivity  of  the  overlying  water.   This  corresponds  with 
the  key  areas  of  research  needs  identified  in  a  comparative 
study  of  a  large  number  of  lake  ecosystem  models  (Scavia 
1979) . 

Modeling  has  advanced  the  understanding  of  this 
aquatic  system  by:  (1)  pointing-out  errors  in  thinking  when 
the  model  did  not  adequately  reproduce  what  was  known  to  be 
true  about  the  behavior  of  the  real  system,  (2)  requiring 
explicit  statements  about  separate  processes  and  the  way  in 
which  they  interact,  (3)  encouraging  continued  feedback 
between  the  model  and  the  real  system,  and  (4)  by  forcing 
the  integration  of  theory  from  disciplines  (e.g.  hydrology, 
rheology,  ecology,  lake  management)  not  routinely  combined. 


LITERATURE  CITED 

Adkins,  M.W.  1991.   Investigation  of  the  water  resources  in 
the  Orange  Lake/ Florida  Ridge  watershed.   St.  Johns 
River  Water  Management  District,  Palatka,  Florida. 

American  Public  Health  Association.  1985.   Standard 

methods  for  the  examination  of  water  and  wastewater. 
Washington,  D.C.:  1193  pp. 

Anderson,  J.M.  1973.   The  eutrophication  of  lakes,  p.  117- 
140.   In  D.L.  Meadows  and  D.H.  Meadows  (eds.),  Toward 
global  equilibrium:   Collected  papers,  Wright -Allen 
Press,  Cambridge,  MA. 

Anderson,  R.F.,  S.L.  Schiff,  and  R.H.Hesslein.  1987. 

Determining  sediment  accumulation  and  mixing  rates 
using  ^•LUPb,  137Cs,  and  other  tracers:  Problems  due  to 
postdepositional  mobility  or  coring  artifacts.   Can. 
J.  Fish.  Aquat.  Sci.  44:  231-250. 

Appleby,  P.G.,  P.J.  Nolan,  D.W.  Gifford,  M.J.  Godfrey,  P. 
Oldfield,  N.J.  Anderson,  and  R.W.  Battarbee.  1986. 

Pb  dating  by  low  background  gamma  counting. 
Hydrobiologia  143:  21-27. 

Appleby,  P.G.  and  F.  Oldfield.  1978.   The  calculation  of 
lead-210  dates  assuming  a  constant  rate  of  supply  of 
unsupported  210Pb  to  the  sediment.   Catena  5:  1-8. 

Appleby,  p.g.  and  F.  Oldfield.  1983.  The  assessment  of 

Pb  data  from  sites  with  varying  sediment 

accumulation  rates.   Hydrobiologia  103:  29-35. 

Baker,  L.A. ,  P.L.  Brezonik,  and  C.R.  Kratzer.  1981. 
Nutrient  loading-trophic  state  relationships  in 
Florida  lakes.   Florida  Water  Resources  Research 
Center  publ.  56,  Gainesville,  Florida. 

Baxter,  R.M.  1977.   Environmental  effects  of  dams  and 
impoundments.   Ann.  Rev.  Ecol.  Syst.  8:  255-283. 


160 


161 

Bays,  J.S.  and  T.L.  Crisman.  1983.   Zooplankton-trophic 
state  relationships  in  Florida  lakes.   Can.  J.  Fish. 
Aquat.  Sci.  40:  1813-1819. 

Bedinger,  M.S.  1979.   Relation  between  forest  species  and 

flooding.   In  P.E.  Greeson,  J.R.  Clark,  and  J.E.  Clark 
(eds.),  Wetland  functions  and  values:  The  state  of  our 
understanding.  Amer.  Water  Res.  Assoc,  Minneapolis, 
MN:  427-435. 

Beeton,  A.M.  1958.   Relationship  between  Secchi  disk 

readings  and  light  penetration  in  Lake  Huron.  Trans. 
Amer.  Fish.  Soc.  87:  73-79. 

Bengtsson,  L.  1975.   Phosphorus  release  from  a  highly 
eutrophic  lake  sediment.   Verh.  Internat.  Verein. 
Limnol.  19:  1107-1116. 

Bengtsson,  L. ,  T.  Hellstr6m,  and  L.  Rakoczi.  1990. 

Redistribution  of  sediments  in  three  Swedish  lakes. 
Hydrobiologia  192:  167-181. 

Binford,  M.W.  1990.  Calculation  and  uncertainty  analysis 
of  210Pb  dates  for  PIRLA  project  lake  sediment  cores. 
J.  Paleolimnol.  3:  253-267. 

Binford,  M.W.  and  M.  Brenner.  1986.   Dilution  of  210Pb  by 
organic  sedimentation  in  lakes  of  different  trophic 
states,  and  application  to  studies  of  sediment-water 
interactions.   Limnol.  Oceanogr.  31(3):  584-595. 

Bishop,  E.W.  1967.  Florida  lakes — Part  I.  Study  of  the 
high-water  lines  of  some  Florida  lakes.  Division  of 
Water  Resources,  Florida  Board  of  Conservation, 
Tallahassee. 

Bowman,  G.T.  and  J.J.  Delfino.  1980.   Sediment  oxygen 
demand  techniques:  A  review  and  comparison  of 
laboratory  and  in  situ  systems.   Water  Res.  14:  491- 
499. 

Bremner,  J.M.  and  C.S.  Mulvaney.  1982.   Nitrogen- total.   In 
A.L.  Page,  R.H.  Miller,  and  D.R.  Keeney  [eds.], 
Methods  of  soil  analysis,  part  II,  chemical  and 
microbiological  properties.   Amer.  Soc.  Agro.  &  Soil 
Sci.  Soc.  Amer.,  Madison,  WI:  599-609. 

Brenner,  M.  and  M.W.  Binford.  1986.   Material  transfer  from 
water  to  sediment  in  Florida  lakes.   Hydrobiologia 
143:  55-61. 


162 

Bruno,  N.A.   1984.   Nest  site  selection  by  and  spawning 

season  vegetation  associations  of  largemouth  bass  in 
Orange  Lake,  Florida.   M.S.  thesis,  University  of 
Florida,  Gainesville. 

Byrkit,  D.R.   1975.   Elements  of  statistics,  2nd  edition. 
D.  van  Nostrand  Publ.  Co.,  New  York,  NY:  431  pp. 

Canfield,  D.E.  Jr.  1981.   Chemical  and  trophic  state 
characteristics  of  lakes  in  relation  to  regional 
geology.   Final  Report  to  Coop.  Fish  and  Wildlife  Res. 
Unit,  University  of  Florida,  Gainesville. 

Canfield,  D.E.  Jr.,  J.V.  Shireman,  D.E.  Colle,  W.T.  Haller, 
C.E.  Watkins,  and  M.J.  Maceina.   1984.   Prediction  of 
chlorophyll  a  concentrations  in  Florida  lakes: 
Importance  of  aquatic  macrophytes.   Can.  J.  Fish. 
Aguat.  Sci.  41:  497-501. 

Chant,  L.A.  and  R.J.  Cornett.  1991.  Smearing  of  gravity 
core  profiles  in  soft  sediments.  Limnol.  Oceanogr. 
36(7):  1492-1498. 

Clark,  W.,  R.  Musgrove,  C.  Menke,  and  J.  Cagle.  1964. 
Interim  Report  on  the  water  resources  of  Alachua, 
Bradford,  Clay,  and  Union  counties,  Florida.   Florida 
Geological  Inform.  Circular  35,  State  Board  of 
Conservation,  Division  of  Geology,  Tallahassee, 
Florida. 

Conrow,  R.  1984.   Habitat  preferences  and  seasonal 

succession  of  early  life  stages  of  fishes  in  Orange 
Lake,  Florida,  with  an  evaluation  of  sampling  methods. 
M.S.  thesis,  University  of  Florida,  Gainesville. 

Cooke,  G.D.  1980.   Lake  level  drawdown  as  a  macrophyte 
control  technique.   Water  Res.  Bull.  16:317-322. 

Cooke,  G.D.,  M.R.  McComas,  D.W.  Waller,  AND  R.H.  Kennedy. 
1977.   The  occurrence  of  internal  phosphorus  loading 
in  two  small  eutrophic  lakes  in  northeastern  Ohio. 
Hydrobiologia  56:  129-135. 

Cooke,  G.D.,  E.B.  Welch,  S.A.  Peterson,  and  P.R.  Newroth. 
1986.   Lake  and  reservoir  restoration.   Butterworth 
Publ.,  Stoneham,  Massachusetts. 

Crider,  D.E.  1972.   Newnan's  Lake;  its  history,  current 
environmental  problems  and  their  solution.   Florida 
Game  and  Fresh  Water  Fish  Commission,  open  file 
report,  Lake  City:  6  pp. 


163 

Crisman,  T.L.  1978.   Reconstruction  of  past  lacustrine 
environments  based  on  the  remains  of  aquatic 
invertebrates.   In  D.  Walker  and  J.C.  Guppy  (eds.), 
Biology  and  quaternary  environments.   Acad.  Sci., 
Canberra,  Australia:  69-101. 

Crisman,  T.L.  1986a.   Algal  control  through  trophic-level 
interactions:  Investigations  at  Lakes  Wauburg  and 
Newnan's,  Florida.   Florida  Department  of  Natural 
Resources,  Tallahassee,  Florida. 

Crisman,  T.L.   1986b.   Eutrophication  control  with  an 
emphasis  on  macrophytes  and  algae.   In  N.  Polunin 
(ed.),  Ecosystem  theory  and  application.   John  Wiley  & 
Sons,  Ltd.,  London:  200-239. 

Davis,  J.H.  Jr.  1973.   Establishment  of  mean  high-water 
lines  in  Florida  lakes.   Univ.  of  Florida  Water 
Resources  Research  Center,  publ.  24,  Gainesville. 

Davis,  M.B.  1968.   Pollen  grains  in  lake  sediments: 

Redeposition  caused  by  seasonal  water  circulation. 
Science  162:  796-799. 

Davis,  M.B.  and  M.S.  Ford.  1982.   Sediment  focusing  in 

Mirror  Lake,  New  Hampshire.   Limnol.  Oceanogr.  27(1): 
137-150. 

Day,  J.W.  Jr.,  C.A.S.  Hall,  W.M.  Kemp,  and  A.Yaftez- 

Arancibia.  1989.   Estuarine  ecology.  J.  Wiley  &  Sons, 
New  York. 

Deevey,  E.S.,  M.W.  Binford,  M.  Brenner,  T.J.  Whitmore. 
1986.   Sedimentary  records  of  accelerated  nutrient 
loading  in  Florida  lakes.   Hydrobiologia  143:  49-53. 

Dillon,  P.J.  and  F.H.  Rigler.  1974.   A  test  of  a  simple 
nutrient  budget  model  predicting  the  phosphorus 
concentration  in  lake  water.   J.  Fish.  Res.  Board  Can. 
31:  1771-1778. 

Durant,  D.F.  1980.   Fish  distribution  among  habitats  in 

hydrilla  infested  Orange  Lake,  Florida.   M.S.  thesis, 
University  of  Florida,  Gainesville. 

Eakins,  J.D.  and  R.T  Morrison.  1978.   A  new  procedure  for 
the  determination  of  lead-210  in  lake  and  marine 
sediments.   Int.  J.  Applied  Rad  Isotopes  29:  531-536. 

Eddy,  S.  and  J.C.  Underhill.  1978.   How  to  know  the 

freshwater  fishes.   Brown  Co.  Publ.,  Dubuque,  Iowa. 


164 

Engstrom,  D.R.  and  H.E.  Wright,  Jr.  1984.   Chemical 
stratigraphy  of  lake  sediments  as  a  record  of 
environmental  change.   In  E.Y.  Haworth  and  J.W.G.  Lund 
[eds.],  Lake  sediments  and  environmental  history. 
Leicester  Univ.  Press,  Leicester,  U.K.:  11-67. 

Environmental  Protection  Agency.   1979.   Methods  for  the 
chemical  analysis  of  water  and  wastes.   Method 
351.2-Storet  00625,  E.P.A.,  Corvalis,  OR. 

Evans,  R.D.  and  F.H.  Rigler.  1980.   Measurement  of  whole- 
lake  sediment  accumulation  and  phosphorus  retention 
using  lead-210  dating.   Can.  J.  Fish.  Aquat.  Sci.  37: 
817-822. 

Ferrante,  J.G.  1976.   The  role  of  zooplankton  in  the 
intrabiocoenotic  phosphorus  cycle  and  factors 
affecting  phosphorus  excretion  in  a  lake. 
Hydrobiologia  49:  203-214. 

Florida  Game  and  Fresh  Water  Fish  Commission.  1982-1989. 
Fish  management  reports  for  the  northeast  Florida 
region.   F.G.F.W.F.C. ,  Lake  City,  Florida. 

Fontaine,  T.D.  1978.   Community  metabolism  patterns  and  a 
simulation  model  of  a  lake  in  central  Florida.   Ph.D. 
dissertation,  University  of  Florida,  Gainesville. 

Forrester,  J.W.  1961.   Industrial  dynamics.   MIT  Press, 
Cambridge,  MA:  464  pp. 

Forrester,  J.W.  1968.   Principles  of  systems.  Wright -Allen 
Press,  Cambridge,  MA. 

Fox,  J.L.,  P.L.  Brezonik,  and  M.A.  Keirn.   1977.   Lake 
drawdown  as  a  method  of  improving  water  quality. 
Ecological  Research  Series,  E. P. A. -600/3-77005, 
Environmental  Protection  Agency,  Atlanta,  GA. 

Fuhs,  G.W.,  S.D.  Demmerle,  E.  Canelli,  and  M.  Chen.  1972. 
Characterization  of  phosphorus-limited  plankton  algae 
(with  reflections  on  the  limiting  nutrient  concept) . 
In  G,E'  Likens  (ed.),  Nurients  and  eutrophication:  The 
limiting  nutrient  controversy.   Special  symposium, 
Amer.  Soc.  Limnol.  Oceanogr. ,  Michigan  State 
University,  Hickory  Corners,  1:  113-133. 

Fukuda,  M.K.  and  W.  Lick.  1980.   The  entrainment  of 

cohesive  sediments  in  freshwater.   J.  Geophys.  Res. 
85:  2813-2824. 


165 

Gale,  P.M.,  K.R.  Reddy,  and  D.A.  Graetz.  1992. 

Mineralization  of  sediment  organic  matter  under  anoxic 
conditions.   J.  Env.  Qual.  21:  (in  press) 

Gates,  D.M.  1963.   The  energy  environment  in  which  we  live. 
Amer.  Scientist  51(4):  327-347. 

Gayle,  T.L.  1975.   Systems  models  for  understanding 
eutrophication  in  Lake  Okeechobee.   M.S.  thesis, 
University  of  Florida,  Gainesville. 

Goldberg,  E.D.  1963.   Geochronology  with  210Pb.  In 
Radioactive  Dating,  IAEA,  Vienna:  121-131.  " 

Goodrick,  R.L.  and  J.F.  Milleson.   1974.   Studies  of 

floodplain  vegetation  and  water  level  fluctuation  in 
the  Kissimmee  River  Valley.   Central  and  Southern 
Flood  Control  District,  West  Palm  Beach,  Florida:  60 
pp. 

Gottgens,  J.F.   1987.   Water  resources  of  the  Prairie  Creek 
area.   In  M.T.  Brown  and  R.G.  Hamann  (eds.),  The 
Prairie  Creek  CARL  project.   Alachua  County 
Conservation  and  Recreation  Task  Force,  Gainesville, 
Florida:  Dl-4. 

Gottgens,  J.F.  and  T.L.  Crisman.  1991.   Newnan's  Lake, 
Florida:  Removal  of  particulate  organic  matter  and 
nutrients  using  a  short-term  partial  drawdown.   Lake 
and  Reserv.  Manage.  7(1):  53-60. 

Gottgens,  J.F.  and  T.L.  Crisman.  1992.   Sediments  of 
Newnan's  Lake:  Characteristics  and  patterns  of 
redistribution  following  a  short-term  drawdown.   St. 
Johns  River  Water  Management  District,  Palatka, 
Florida. 

Gottgens,  J.F.  and  C.L.  Montague.  1987a.   Orange, 
Lochloosa,  and  Newnan's  lakes:  A  survey  and 
preliminary  interpretation  of  environmental  research 
data.   St.  Johns  River  Water  Management  District, 
Palatka,  Florida. 


Gottgens,  J.F.  and  C.L.  Montague.  1987b.  A  categorized 
bibliography  of  the  Orange  Creek  Basin.  St.  Johns 
River  Water  Management  District,  Palatka,  Florida. 


Gottgens,  J.F.  and  C.L.  Montague.  1988.   Comprehensive 

reconnaissance  profile  of  the  Payne's  Prairie  Basin, 
Florida.  St.  Johns  River  Water  Management  District, 
Palatka,  Florida. 


166 

Gosselink,  J.G.  and  P.E.  Turner.  1978.  The  role  of 

hydrology  in  freshwater  wetland  ecosystems.   In  R.E. 
Rood,  D.F.  Whigham,  and  R.L.  Simpson  (eds.), 
Freshwater  wetlands:  Ecological  processes  and 
management  potential.   Academic  Press,  New  York:  63- 
78. 

Goulden,  C.E.  1969.   Interpretative  studies  of  cladoceran 
microfossils  in  lake  sediments.   Mitt.  Int.  Ver. 
Limnol.  17:  43-55. 

Greening,  H.  and  S.  Doyon.  1990.   Environmental  and 
ecological  effects  of  drawdown  and  enhanced 
fluctuation  for  Lake  Apopka,  Florida.   St.  Johns  River 
Water  Management  District,  Palatka,  Florida. 


Griffiths,  M.  and  W.T.  Edmondson.  1975.   Burial  of 

oscillaxanthin  in  the  sediment  of  Lake  Washington. 
Limnol.  Oceanogr.  20  (6):  945-952. 


Gunderson,  L.H.  1989.  Historical  hydropatterns  in  wetland 
communities  of  Everglades  National  Park.   In  R.R. 
Sharitz  and  J.W.  Gibbons  (eds.),  Freshwater  wetlands 
and  wildlife.   DOE  symposium  series  61,  U.S. 
Department  of  Energy.  Oak  Ridge,  Tennessee:  1099-1111. 

Gutierrez,  L.T.  and  W.R.  Fey.  1975.   Simulation  of 

secondary  autogenic  succession  in  the  shortgrass 
prairie  ecosystem.   Simulation  24:  113-125. 

Hakanson,  L.  1982.   Bottom  dynamics  in  lakes. 
Hydrobiologia  91:  9-22. 

Hamilton,  D.P.  and  S.F.  Mitchell.  1988.   Effects  of  wind  on 
nitrogen,  phosphorus,  and  chlorophyll  a  in  a  shallow 
New  Zealand  lake.  Verh.  Int.  Verein.  Limnol.  23:  624- 
628. 

Hargrave,  B.T.  1969.   Similarity  of  oxygen  uptake  by 

benthic  communities.   Limnol.  Oceanogr.  14:  801-805. 

Harris,  S.W.  and  W.H.  Marshall.  1963.   Ecology  of  water- 
level  manipulations  on  a  northern  marsh.   Ecology  44: 
331-343. 


Hestand,  R.S.  and  C.C.  Carter.   1975.   Succession  of 
aguatic  vegetation  in  Lake  Ocklawaha  two  growing 
seasons  following  a  winter  drawdown.   Hyacinth  Contr. 
J.  13:  43-47. 


167 

Holcomb,  D.  and  W.  Wegener.   1971.   Hydrophytic  changes 
related  to  lake  fluctuation  as  measured  by  point 
transects.   Proc.  25th  Ann.  Conf .  Southeastern  Assoc. 
Game  Fish  Comm. :  570-583. 

Holcomb,  D.E.,  W.L.  Wegener  and  V.P.  Williams.  1975.   Lake 
level  fluctuation  for  habitat  management:  A  case  in 
point.   In  P.L.  Brezonik  and  J.L.  Fox  (eds.),  Water 
quality  management  through  biological  control. 
University  of  Florida,  Gainesville:  15-24. 

Holdren,  G.C.  Jr.  and  D.E.  Armstrong.  1980.   Factors 

affecting  phosphorus  release  from  intact  lake  sediment 
cores.   Environ.  Sci.  Technol.  14  (1):  79-87. 

Holly,  J.B.  1976.   Stratigraphy  and  sedimentary  history  of 
Newnan's  Lake.   M.S.  thesis,  University  of  Florida, 
Gainesville. 

Huber,  W.C.,  P.L.  Brezonik,  J. P.  Heaney,  R.E.  Dickinson, 

S.D.  Preston,  D.S.  Dwornik,  and  M.A.  Demaio.  1982.   A 
classification  of  Florida  lakes.   Florida  Water 
Resources  Research  Center  Publ.  72,  Gainesville. 

Hwang,  K.N.  and  A.J.  Mehta.  1989.   Behavior  of  fine 

sediment  in  generating  wave-induced  turbidity.  Proc. 
Beach  Preserv.  Technol.,  Florida  Shore  Beach  Preserv. 
Assoc,  Tallahassee,  Florida:  207-218. 

James,  A.  1974.   The  measurement  of  benthal  respiration. 
Water  Res.  8:  955-959. 

Johansen,  D.  1970.  Air  breathing  in  fishes.   In  w.s.  Hoar 
and  D.J.  Randall  (eds.),  Fish  physiology,  vol.  IV. 
Academic  Press,  New  York:  361-411. 

KBN  Engineering  and  Applied  Sciences.  1991.   Assessing  the 
feasibility  of  restoring  Newnan's  Lake.   KBN, 
Gainesville,  Florida. 

Keddy,  P. A.  and  A. A.  Reznicek.  1985.   Vegetation  dynamics, 
buried  seeds,  and  water  level  fluctuations  on  the 
shorelines  of  the  Great  Lakes.   In  Coastal  wetlands, 
Lewis  Publ.,  Chelsea,  Michigan:  33-58. 

Knoll,  G.F.  1979.   Radiation  detection  and  measurement.   J. 
Wiley  &  Sons,  New  York,  NY:  762  pp. 

Kramer,  K.J.M.,  R.  Misdorp,  G.  Berger,  and  R  Duijts.  1991. 
Maximum  pollutant  concentrations  at  the  wrong  depth:  A 
misleading  pollution  history  in  a  sediment  core. 
Marine  Chemistry  36:  183-198. 


168 

Krishnaswami ,  S  and  D.  Lai.  1978.   Radionuclide 

limnochronology.   In  A.  Lerman  (ed.),  Lakes: 
Chemistry,  geology,  Physics.   Springer-Verlag,  New 
York:  153-177. 

Kuenzler,  E.J.  and  B.H.  Ketchum.  1962.   Rate  of  phosphurus 
uptake  by  Phaeodactylum  tricornutum.   Biol.  Bull.  123: 
134-145. 

Kuhl,  A.  1974.  Phosphorus.  In  W.D.P.  Stewart  (ed.),  Algal 
physiology  and  biochemistry.  University  of  California 
Press,  Berkeley:  636-654. 

Kushlan,  J. A.  1974.   Effects  of  a  natural  fish  kill  on  the 
water  quality,  plankton,  and  fish  population  of  a  pond 
in  the  Big  Cypress  Swamp,  Florida.   Trans.  Amer.  Fish. 
SOC.  103:  235-243. 

Lam,  D.C.L.  and  J.M.  Jaquet.  1976.   Computations  of 

physical  transport  and  regeneration  of  phosphorus  in 
Lake  Erie,  Fall  1970.   J.  Fish.  Res.  Board  Can.  33: 
550-563. 

Lantz,  K.E.,  J.T.  Davis,  J.S.  Hughes,  and  H.E.  Schafer. 
1964.   Water  level  fluctuation — Its  effects  on 
vegetation  control  and  fish  population  management. 
Proc.  Southeast.  Assoc.  Game  Fish  Comm.  18:  483-494. 

Lee,  CM.  and  D.L.  Macalady.  1989.   Towards  a  standard 
method  for  the  measurement  of  organic  carbon  in 
sediments.   Intern  J.  Environ.  Anal.  Chem.  35:  219- 
225. 

Lee,  D.Y.,  W.  Lick,  and  S.W.  Kang.  1981.   The  entrainment 

and  deposition  of  fine-grained  sediments  in  Lake  Erie. 
J.  Great  Lakes  Res.  7(3):  224-233. 

Lewis,  W.M.,  Jr.  1970.  Morphological  adaptations  of 
cyprinodontoids  for  inhabiting  oxygen  deficient 
waters.   Copeia  1970:  319-326. 

Li,  Y.H.,  D.E.  Armstrong,  J.D.H.  Williams,  R.F.  Harris,  and 
J.K.  Seyers.  1972.   Rate  and  extent  of  inorganic 
phosphate  exchange  in  sediments.   Soil  Sci.  Soc.  Amer. 
Proc.  36:  279-285. 

Lick,  W.  1982.   Entrainment,  deposition,  and  transport  of 
fine-grained  sediments  in  lakes.   Hydrobiologia  91: 
31-40. 

Livingstone,  D.A.  1955.   A  lightweight  piston  sampler  for 
lake  deposits.   Ecology  36:  137-139. 


169 

Loftus,  W.F.  and  J. A.  Kushlan.  1987.   Freshwater  fishes  of 
southern  Florida.   Bull.  Florida  State  Museum,  Biol. 
Sci.  31(4) :  344  pp. 

Luettich,  R.A.Jr.,  D.R.F.  Harleman,  and  L.  Somlyody.  1990. 
Dynamic  behavior  of  suspended  sediment  concentrations 
in  a  shallow  lake  perturbed  by  episodic  wind  events. 
Limnol.  Oceanogr.  35(5):  1050-1067. 

Maceina,  M.J.  and  D.M.  Soballe.   1990.   Wind  related 

limnological  variation  in  Lake  Okeechobee,  Florida. 
Lake  Reserv.  Manag.  6(1):  93-100 

McCormack,  B.  1967.   Aerial  respiration  in  the  Florida 

spotted  gar.   Quart.  J.  Florida  Acad.  Sci.  30:  68-72. 

McDonnell,  A.J.  and  S.D.  Hall.  1969.   Effects  of 

environmental  factors  on  benthal  oxygen  uptake.   J. 
Wat.  Pollut.  Control  Fed.  41:  R353-R363. 


McKinney,  S.P.  and  W.S.  Coleman.  1980.   Hydrilla  control 
and  vegetation  responses  with  multiple  dewaterings. 
Florida  Game  and  Fresh  Water  Fish  Commission, 
Kissimmee,  Florida. 

Mickle,  A.M.  and  R.G.  Wetzel.  1978.   Effectiveness  of 

submersed  angiosperm-epiphyte  complexes  on  exchange  on 
nutrients  and  organic  carbon  in  littoral  systems.   I. 
Inorganic  nutrients.   Aguatic  Bot.  4:  303-316. 

Montague,  C.L.,  W.R.  Fey  and  D.M.  Gillespie.  1982.   A 

causal  hypothesis  explaining  predator-prey  dynamics  in 
Great-Salt  lake,  Utah.  Ecological  Modelling  17:  243- 
270. 

Mortimer,  C.H.  1971.   Chemical  changes  between  sediments 
and  water  in  the  Great  Lakes — Speculations  on 
probable  regulatory  mechanisms.   Limnol.  Oceanogr.  16: 
387-404. 

Nagy,  J.W.  1988.   Simultaneous  determination  of  supported 
and  unsupported  210Pb  concentrations  in  lake-bottom 
sediments  using  low-energy,  high-purity  germanium 
gamma-ray  spectroscopy.   Department  of  Environmental 
Engineering  Sciences,  University  of  Florida, 
Gainesville. 

Nordlie,  F.G.  1976.   Plankton  communities  of  three  central 
Florida  lakes.   Hydrobiologia  48  (1):  65-78. 

Nozaki,  Y.,  D.J.  DeMaster,  D.M.  Lewis,  and  K.K.  Turekian. 
1978.  Atmospheric  2IOPb  fluxes  determined  from  soil 
profiles.   J.  Geophys.  Res.  83:  4047-4051. 


170 

Odum,  E.P.  1971.   Fundamentals  of  ecology.   3rd  edition. 
W.B.  Saunders  Co.,  Philadelphia,  Pennsylvania. 

Odum,  H.T.  1983.   Systems  ecology.   Wiley  &  Sons,  New  York, 
NY:  644  pp. 

Odum,  H.T.  1988.   Energy,  environment,  and  public  policy:  A 
guide  to  the  analysis  of  systems.   UNEP  Regional  Seas 
Report  and  Studies  95,  New  York,  NY:  109  pp. 

Oldfield,  F.,  P.G.  Appleby,  and  R.  Thompson.  1980. 

Palaeoecological  studies  of  the  lakes  in  the  highlands 
of  Papua  New  Guinea:  I.   The  chronology  of 
sedimentation.   J.  Ecol.  68:  457-477. 

Opper,  S.C.  1982.   The  hydrogeology  of  Lake  Wauburg  and 
vicinity,  Alachua  County,  Florida.   M.S.  thesis, 
University  of  Florida,  Gainesville. 

Paul,  J.F.,  R.  Kasprzyk,  and  W.  Lick.  1982.   Turbidity  in 
the  western  basin  of  Lake  Erie.   J.  Geophys.  Res.  87: 
5779-5784. 

Pennington,  W.  1973.   The  recent  sediments  of  Windermere. 
Freshwat.  Biol.  3:  363-382. 

Pennington,  W.  1981.  Records  of  a  lake's  life  in  time:  The 
sediments.  Hydrobiologia  79:  197-219. 

Peters,  R.H.  1975.   Orthophosphate  turnover  in  central 

European  lakes.   Mem.  1st.  Ital.  Idrobiol.  32:  297- 
311. 

Pilskaln,  C.H.  and  T.C.  Johnson.  1991.   Seasonal  signals  in 
Lake  Malawi  sediments.   Limnol.  Oceanogr.  36(3):  544- 
557. 

Pirkle,  E.C.  and  H.K.  Brooks.  1959.   Origin  and  hydrology 
of  Orange  Lake,  Santa  Fe  Lake,  and  Levy  Prairie  Lakes 
of  north-central  peninsular  Florida.   J.  of  Geology  67 
(3):  302-317. 

Pollman,  CD.  1983.   Internal  loading  in  shallow  lakes. 
Ph.D.  dissertation,  University  of  Florida, 
Gainesville. 

Reddy,  K.R.  and  W.H.  Patrick.  1984.   Nitrogen 

transformations  and  loss  in  flooded  soils  and 
sediments.   CRC  Critical  Rev.  Environ.  Contr.  13(4): 
273-309. 


171 

Redfield,  A.C.,  B.H.  Ketchum,  and  F.A.  Richards.  1963.  The 
influence  of  organisms  on  the  composition  of  seawater. 
In  M.N.  Hill  (ed.),  The  sea.   Interscience,  New  York, 
Vol.  2:  26-77. 

Richardson,  L.V.  1975  Water  level  manipulation:  A  tool  for 
aquatic  weed  control.   Hyacinth  Control  J.  13:  8-11. 

Ritchie,  J.C.,  J.R.  McHenry,  and  A.C.  Gill.  1973.   Dating 
recent  reservoir  sediments.   Limnol.  Oceanogr.  18(2): 
254-263. 

Robbins,  J. A.  and  D.N.  Edgington.  1975.   Determination  of 
recent  sedimentation  rates  in  Lake  Michigan  using  Pb- 
210  and  Cs-137.   Geochim.  Cosmochim.  Acta  39:  285-304. 

Robbins,  J. A. ,  D.N.  Edgington,  and  A.L.W.  Kemp.  1978. 

Comparative  210Pb,  13'Cs,  and  pollen  geochronologies  of 
sediments  from  lakes  Ontario  and  Erie.  Quat.  Res.  10: 
256-278. 

Saunders,  G.W.,  K.W.  Cummins,  D.Z.  Gak,  E.  Pieczynska,  V. 
Straskrabova ,  and  R.G.  Wetzel.  1980.   Organic  matter 
and  decomposers.  In  E.D.  LeCren  and  R.H.  Lowe- 
McConnell  (eds.),  The  functioning  of  freshwater 
ecosystems.   Cambridge  Univ.  Press,  Cambridge:  341- 
392. 

Scavia,  D.  1979.   The  use  of  ecological  models  of  lakes  in 
synthesizing  available  information  and  identifyig 
research  needs.   In  D.  Scavia  and  A.  Robertson  (eds.), 
Perspectives  on  lake  ecosystem  modeling.   Ann  Arbor 
Science  Publ.,  Ann  Arbor,  Michigan:  109-169. 

Schramm,  H.L.,  M.V.  Hoyer,  and  K.J.  Jirka.  1983.   Relative 
ecological  value  of  common  aquatic  plants.   Final 
Report  to  Florida  Dept.  of  Natural  Resources, 
Tallahassee,  Florida. 

Sederholm,  H. ,  A.  Mauranen,  and  L.  Montonen.   1973.   Some 
observations  on  the  microbial  degradation  of  humous 
substances  in  water.   Verh.  Int.  Verein.  Limnol.  18: 
1301-1305. 

Serruya,  C.  and  U.  Pollingher.  1977.   Lowering  of  water 
level  and  algal  biomass  in  Lake  Kinneret. 
Hydrobiologia  54:  73-80. 

Shannon,  E.E.  1970.   Eutrophication-trophic  state 

relationships  in  north  and  central  Florida  lakes. 
Ph.D.  dissertation,  University  of  Florida, 
Gainesville. 


172 

Shapiro,  J.,  W.T.  Edmondson,  and  D.E.  Allison.  1971. 

Changes  in  the  chemical  composition  of  sediments  of 
Lake  Washington,  1958-1970.   Limnol.  Oceanogr.  16  (2): 
437-452. 

Sheng,  P.Y.  and  W.  Lick.   1979.   The  transport  and 

resuspension  of  sediments  in  a  shallow  lake.   J. 
Geophys.  Res.  84:  1809-1826. 

Skoglund.  1990.   Newnans  Lake  profile.   Florida  Game  and 
Fresh  Water  Fish  Commission,  Lakeland:  41  pp. 
(+  appendices) 

Smeltzer,  E.  and  E.B.  Swain.  1985.   Answering  lake 

management  questions  with  paleolimnology.   Proc.  4th 
Ann.  Conf.  North  Amer.  Lake  Manag.  Soc. ,  Williamsburg, 
VA:  268-274. 

Smith,  D.B.  and  P  McGriff.  1958.   Water  control  program  for 
the  Newnans,  Lochloosa,  and  Orange  lakes  drainage 
basin.   Alachua  County  Recreation  and  Water 
Conservation  and  Control  Authority,  Gainesville, 
Florida. 

Somlyody,  L.  1982.   Water  quality  modelling:  A  comparison 
of  transport-oriented  and  ecology-oriented  approaches. 
Ecol.  Modelling  17:  183-207. 

Stauffer,  R.E.  1981.   Sampling  strategies  for  estimating 
the  magnitude  and  importance  of  internal  phosphorus 
supplies  in  lakes.   Report  #  600/3-81-015.   U.S. 
Environmental  Protection  Agency,  Corvalis,  Oregon. 

Steemann  Nielsen,  E.   1962.   Inactivation  of  the 

photochemical  mechanism  in  photosynthesis  as  a  means 
to  protect  the  cells  against  too  high  light 
intensities.  Physiol.  Plant  15:  161-171. 

Symoens,  J.J.,  M.  De  Clercq,  R.A.  Osafo,  and  U.  Sansen. 

1988.   The  lake  of  Virelles  (Prov.  Hainaut,  Belgium): 
Management  of  a  lake  for  nature  conservation  and 
recreation.   Verh.  Internat.  Verein.  Limnol.  23:  1057- 
1062. 

Tarver,  D.P.  1980.  Water  fluctuation  and  the  aquatic  flora 
of  Lake  Miccosukee.  J.  of  Aquat.  Plant  Manag.  18:  19- 
23 . 

Terwindt,  J.H.J.  1977.   Deposition,  transportation,  and 

erosion  of  mud.   In  H.L.  Golterman  (ed.),  Interactions 
between  sediments  and  freshwater.   Junk,  Dordrecht, 
The  Netherlands:  19-24. 


173 

Theis,  T.L.  and  P.J.  McCabe,  1978.   Phosphorus  dynamics  in 
hypereutrophic  lake  sediments.   Water  Res.  12:  677- 
685. 

U.S.  Army  Coastal  Engineering  Research  Center.  1977.   Shore 
protection  manual,  v.  I.   U.S.  Army  Coastal 
Engineering  Research  Center,  Vicksburg,  MI. 

U.S.  Geological  Survey.  1966,  1988.   Topographic  maps,  7.5 
minute  series,  scale  1:24,000.   Gainesville  East 
(29082-F3-TF-024) ,  Micanopy  (N2930-W8215/7 . 5) ,  Orange 
Heights  (29082-F2-TF-024) ,  and  Rochelle  (N2930- 
W8207.5/7.5)  quadrangles,  Reston,  VA. 


Vallentyne,  J.R.  1970.   Phosphorus  and  the  control  of 

eutrophication.   Can.  Res.  Development  1970:  36-43. 

Vollenweider,  R.A.  1976.   Advances  in  defining  critical 

loading  levels  for  phosphorus  in  lake  eutrophication. 
Mem.  1st.  Ital.  Idrobiol.  33:  53-83. 

Walters,  C.J.,  Gunderson,  L.H.  and  C.S.  Holling.  1992. 
Experimental  policies  for  water  management  in  the 
Everglades.  Ecol.  Applications  2:  145-160. 

Walton,  W.C.  1970.   Groundwater  resources  evaluation. 
McGraw-Hill  Book  Co.,  New  York. 

Wang,  C.H. ,  D.L.  Willis,  and  W.D.  Loveland.  1975. 
Radiotracer  methodology  in  the  biological, 
environmental,  and  physical  sciences.   Prentice-Hall, 
Englewood  Cliffs,  NJ:  480  pp. 

Wegener,  W.  and  V.  Williams.   1974.   Fish  population 
responses  to  improved  lake  habitat  utilizing  an 
extreme  drawdown.   Proc.  28th  Ann.  Conf .  Southeast. 
Assoc.  Game  and  Fish  Comm. ,  Tallahassee,  FL:  144-161. 

Wegener,  W. ,  V.  Williams,  and  T.D.  McCall.   1974.   Aquatic 
macro invertebrate  responses  to  an  extreme  drawdown. 
Proc.  28th  Ann.  Conf.  Southeast.  Assoc.  Game  and  Fish 
Comm.,  Tallahassee,  FL:  126-143. 

Wetzel,  R.G.  1983.   Limnology.   2nd  ed.   CBS  College 
Publishing,  Philadelphia,  PA:  767  pp. 


APPENDIX  A 

LEAD-210  (210PB)  AND  CESIUM-137  (137CS)  DATING  OF  SEDIMENTS 

USING  LOW-ENERGY,  INTRINSIC-GERMANIUM, 

y-RAY  SPECTROSCOPY 

Application 
Direct  vray  assay  is  used  to  record  changing 
concentrations  of  210Pb  (unsupported  and  supported)  and 
±J/Cs  in  sediment  profiles.   This  forms  a  reliable  and 
precise  basis  for  age/depth  and  dry-sedimentation-rate 
calculations. 

Counting  system 
Location 

The  system  used  for  spectral  analysis  is  located  at 
the  University  of  Florida's  Department  of  Environmental 
Engineering  Sciences'  Low-Background  Counting  Room.   This 
room  was  designed  to  reduce  the  level  of  background 
radiation  interference  during  sample  counting. 
Detector  specifications 

Manufacturer  Princeton  Gamma  Tech,  Inc. 

Model  Intrinsic  Germanium  Well 

12225-14 
Preamplifier  RG11B/C 

Crystal  geometry        P-Type  Well 
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Detector  Bias 
Crystal  Size 

Well  Size 
Well  Material 
Bias  Polarity 
Efficiency  ( counts/ y) 


3000  V 

61  mm  (0),  60  mm  (length), 
139  cm3  (active  volume) 
14  mm  (0) ,  40  mm  (length) 
Vespel  (0.51  mm) 
Positive 

0.57  (@  46.5  keV) 
0.27  (§  295.2  keV) 
0.22  (@  351.9  keV) 
0.08  (§  609.3  keV) 
0.07  (@  661.6  keV) 
Standard  Peak/Compton  24:1  -  90:1 
FWTM/FWHM  1.9  -  2.0 

Lowest  Detectable  Energy  9  keV 
Background  Suppression 

Background  suppression  is  accomplished  with  a  10.1  cm 
thick  lead  shield  (low  activity).   In  addition,  evaporated 
liquid  nitrogen  is  piped  into  the  detector  shield  to  create 
excess  pressure  and  displace  ambient  air  during  operation. 
Operating  Conditions 

The  detector  operates  at  a  temperature  near  that  of 
liquid  nitrogen  (-196°C) ,  which  is  provided  by  a  high 
vacuum  cryostat-dewar  system  (capacity  30  1) .   Liquid  N2 
use  is  approximately  2  1/day. 
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Electronics 

System  electronics  include: 

Amplifier  EG  &  G  Ortec,  672 

Bias  Supply  Princeton  Gamma  Tech,  5  kV, 

315A 
Power  Supply  EG  &  G  Ortec,  402M 

MCB  EG  &  G  Ortec,  ACE  4K  PC  Card 

Computer  and  Software 

A  computer  with  multi-channel  analyzer  ("Maestro" 

ADCAM  100,  EG  &  G  Ortec,  0.5  keV/channel)  and  mathematical 

spreadsheet  (Quattro,  Borland  Inc.)  are  used  for  data 

conversion  and  analysis. 

Procedure 

1.  Sediment  cores  ,  kept  at  4°C,  are  sectioned  in  1  cm 
intervals  within  48  hours  after  collection.   Separate 
sections  are  stored  in  plastic  at  4°C. 

2.  Bulk  density  (g/cm3)  is  determined  by  drying  1  cm3 
subsamples  at  95°C  for  24  hours  followed  by  cooling  and 
weighing  on  an  analytical  balance. 

3.  Samples  for  isotope  analysis  are  dried  at  95 °C  for  24 
hours,  pulverized  by  mortar  and  pestle,  re-dried  (95°C,  1 
hour) ,  weighed,  and  placed  in  small,  low-density 
polypropylene  tubes  (capacity  4  ml) .   The  sample  volume  and 
density  should  correspond  to  the  volume  and  density  of  the 
standard  to  insure  the  same  counting  efficiencies  for  both. 
The  samples  and  standard  are  covered  with  an  identical  seal 
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(epoxy)  and  archived  for  14  days  to  equilibrate  radon 
(222Rn)  with  radium  (226Ra)  . 

4.  Counting  times  vary  from  7  to  14  hours  depending  on  the 
weight  of  the  sample;  low  weight  samples  need  longer 
counting  times  to  reduce  error.  A  blank  is  counted  for 
every  three  samples  to  determine  background  levels  of 
radiation.   Standards  are  run  with  the  same  frequency  to 
track  efficiency  (counts/y)  and  calculate  a  226Ra 
conversion  factor  (pCi/cps) . 

5.  Sample  spectra  are  analyzed  for  activity  in  the  46.5 
keV  (210Pb)  and  662  keV  (137Cs)  peaks.  Activities  at  295 
keV  (214Pb)  ,  352  keV  (214Pb) ,  and  609  keV  (214Bi) 
representing  uranium-series  peaks  are  used  to  compute 
supported  levels  of  210Pb. 

Calculations  and  Accuracy 
Radio-nuclide  Concentrations 

The  concentration  (pCi/g)  for  each  y-energy  of 
interest  is  determined  with 

C  =  Pc  "  Bc  (A-l) 

Exlxtxgx  0.037 

where  Pc  =  Peak  counts 

Bc  =  Background  counts 

E  =  Efficiency  (counts/y) 

I  =  Intensity  (y/disintegration) 
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t    =  Count  time  (sec) 
g    =  Sample  weight  (gram) 
and   0.037=  Conversion  factor  from 

disintegrations/sec  to  pCi. 

Counting  efficiency  for  each  peak  is  given  by 


E  =      ?c  I  Bc  (A-2) 

A  x  I  x  t  x  0.037 

where  A    =  Activity  of  the  standard  (pCi) 

The  total  210Pb  concentration  is  calculated  using  Eguation 
(A-l)  applied  to  the  46.52  keV  y  peak.   The  supported  210Pb 
concentration  is  determined  by  Equation  (A-l)  applied  to 
the  295,  352  and  609  keV  peaks,  and  averaging.  The 
unsupported  210Pb  concentration  is  then  calculated  by 

Unsupported  210Pb  =  Total  210Pb  -  Supported  210Pb  (A-3) 


137Cs  concentration  is  determined  by  applying  Equation 

(A-l)  to  the  662  keV  peak. 

Background 

The  background  characteristics  of  the  system  for 
selected  isotopes  were  optimized  by  a  series  (N=13)  of 
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blank  counts  during  a  10  week  period.   Counts  per  minute 
(over  840  minutes)  were  as  follows: 


210pb        2Kpb        2Kpb        2KBi         137Cs 
@  46.5  keV  @  295  keV   @  352  keV   @  609  keV    @  662  keV 


0.02  ±  0.02  0.01  t  0.01  0.04  ±  0.02  0.02  ±  0.02  0.04  ±  0.02 

These  levels  are  sufficiently  low  and  constant  to  achieve 
acceptable  precision  and  counting  times.   Background  counts 
in  the  energy  ranges  of  interest  need  to  be  monitored  after 
every  3-4  sample  counts. 
Error  Prediction 

Analytical  error  prediction  follows  nuclear  statistics 
(Knoll  1979)  assuming  a  Poisson  distribution  for  the 
recorded  counts.   Uncertainty  (pCi/g)  depends  linearly  on 
counting  time  and  sample  weight.   Error  bars  (±  la)  for 
Total  210Pb  were  0.51  pCi/g  (6.3%)  for  a  very  small  sample 
(0.723  g)  with  an  activity  of  8.14  pCi/g  and  a  14  hr.  count 
time.   A  2.087  g  sample  with  an  activity  of  6.41  pCi/g 
counted  for  7  hr.  produced  0.53  pCi/g  (8%)  error  bars. 
These  errors  represent  only  internal  uncertainty  (i.e. 
associated  with  the  accuracy  of  y-ray  detection) . 
Calculation  of  210Pb  Dates 

Calculation  of  ZJ-uPb  dates  follows  the  Constant  Rate 
of  Supply  model  (Goldberg  1963;  Appleby  and  Oldfield  1978; 
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Robbins  1978) .   As  such,  the  cumulative  residual 
unsupported  210Pb,  A^.,  beneath  sediments  of  age  t  varies 
according  to 

At  =  A0  e"kt  (A-4) 


k  At 
r  =  __E  (a-6) 


where  C    =  Concentration  of  unsupported  210Pb 
(pCi/g)  in  layer  of  interest. 


where  A0   =  Total  residual  unsupported  210Pb 


J.UUUJ.      unou^/puj.  ucu 

(pCi/cm2) 
k    =  2iopb  radioactive  decay  constant 

At  and  A0  are  calculated  by  direct  numerical  integration  of 
the  210Pb  profile.  The  age  of  sediments  of  depth  x  is  then 
given  by: 


t  -  1  i»£ 


The  sedimentation  rate  (r)  can  then  be  calculated  directly 
(Appleby  &  Oldfield  1978): 
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Cesium-137 

First  occurrence  of  137Cs  in  the  profile  generally 
coincides  with  the  onset  of  widespread  atmospheric  testing 
of  nuclear  weapons  (i.e.  1954).   1959-1960  and  1963-1964 
are  reported  as  periods  of  maximum  fallout  (Ritchie  et  al. 
1973)  . 
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APPENDIX  B 
MODEL  PROGRAM  LISTING 


Newnan's  Lake  model  -  Simulation  of  the  effect  of  spillway 
installation  on  in-lake  sediment  and  phosphorus  dynamics. 
Time  scale  is  40  years  with  spillway  construction  initiated 
at  15  years. 


DEFDBL  A-Z 


STRING  VARIABLES 


Double  precision 


TITLE$  =  "Effect  of  Spillway  on  Sediment  and  P-Dynamics" 
SUBTITLE$  =  "Newnan's  Lake,  FL" 
T$  =  "Years" 


TIME  PARAMETERS 

TO  =  0 
THALF  =  20 
TEND  =40 
DT  =  .0005 
PLTINTVL  =  .05 


Start  time  (y) 

End  time  (y) 

Integration  interval  (y) 
Plotting  interval  (y) 


INITIAL  VALUES 

PLOTCOUNT  =  0 
PERMSED  =  60000 
SRPWATER  =  .006 
PDETRITUS  =  .016 
PPLANKTON  =  .048 
PINORG  =  .016 
SRPPORE  =  .43 


Permanent  sediments  (g  m~2) 

Water  column  soluble  reactive 
phosphorus  cone.  (gSRP  m"3) 
P-concentration  of  suspended 
dead  organic  matter  (gp  m"3) 
P-concentration  of  suspended 
live  organic  matter  (gp  m-3) 
P-concentration  of  suspended 
inorganic  matter  (gp  m~3) 
SRP  in  Pore  fluids  in  surf sed 


(g 


SRP 


m-<) 
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PDETSED  =  1.64 
PBENTHOS  =  .4 
PINORGSED  =  1.32 
VEGLITAREA  =  20 
LITPART  =200 
DELLITPART  =5000 
AVGLEVEL  =  20.24 


P-concentration  in  dead 

organic  surfsed  (gp  m~2) 

P-concentration  in  live 

organic  surfsed  (gp  m"2) 

P-concentration  in  inorganic 

surfsed  (gp  m"2) 

Percent  vegetated  littoral 

area  in  lake 

Seston  +  surfsed  in  littoral 

(g  *T2) 

Seston  +  surfsed  temp,  trapped 

in  littoral  (g  m-2) 

Average  lake  stage  (m  msl) 


CONSTANTS 
AMPNODAM  =  .55 
AMPDAM  =  .39 
FREQLEVEL  =50 
RISEDAM  =  .13 
RAMPDURATION  =  .5 
BWDNODAM  =  .5 
BWDDAM  =  .3 
HYDRESTMNODAM  =  .6 
HYDRESTMDAM  =  .8 
LAKEBED  =  18.57 
BULKDENS  =  30000 
BEDDENS  =  500000 
SDR  =  .5 
SVELP  =300 
SVELD  =2523 
SVELI  =2523 
ALPHA  =  36.85 
GAMMA  =  -.084 


Amplitude  of  lake  level 
fluctuation  without  dam  (m) 
Amplitude  of  lake  level 
fluctuation  with  dam  (m) 
Freguency  of  water  level 
fluctuation 

Water-level  rise  following  dam 
installation  (m) 
Ramp  coefficient  for 
water-level  rise  (y) 
Concentration  of  seston  in 
outflow  without  dam  (g  m-3) 
Concentration  of  seston  in 
outflow  with  dam  (g  m"3) 
Hydraulic  residence  time 
without  dam  (y) 
Hydraulic  residence  time  with 
dam  (y) 

Initial  elevation  of  sediment- 
water  interface  (m  msl) 
Surfsed  bulk  density  (g  m"3) 

Permsed  bulk  density  (g  m"3) 

Specific  rate  of  deposition 

(%  Y"1) 

Settling  velocity  of  plankton 

(m  y"1) 

Settling  velocity  of  detritus 

(m  y"1) 

Settling  velocity  of  inorganic 

particles (m  y"1) 

Model  coefficient  (for  Srr) 

(%  y"1) 

Model  coefficient  (for  Srr) 


184 


MAXPAR  =16 

OPTPAR  =  8 
BETA  =  .069 
MAXBENTHICBACT  =15.5 

MAXPLANKBACT  =  15.5 

MAXBENTHICAUTO  =  3 
MAXPLANKAUTO  =10 

KM1  =  .31 

KM2  =  .036 

KM3  =  .31 

KM4  =  .036 


Maximum  photosynthetically 

active  radiation  (Par)  at  mean 

lake  depth  (cal  cm"2hr_1) 

Optimum  Par  at  mean  lake  depth 

(cal  cm~2hr_1) 

Extinction  coefficient 

(%  per  g  m-3  of  seston) 

Max.  P-uptake  rate  for  benthic 

bacteria 

(9srp   Spbenthicbact     Y     ) 
Max.    P-uptake  rate  for 

planktonic  bacteria 

(9srp   9ppiankbact     Y     ) 

Max.    P-uptake  rate  for  benthic 

algae    (gSRp  gptenthicauto"1  Y'1) 
Max.    P-uptake  rate  for 
planktonic  algae 

(9srp   gPplankajito     Y     ) 
Half-saturatxon  constant  for 

specific  P-uptake  rate  by 

benthic  bacteria  (gSRP  m--3) 

Half-saturation  constant  for 

specific  P-uptake  rate  by 

benthic  algae  (gSRP  m"3) 

Half-saturation  constant  for 

specific  P-uptake  rate  by 

planktonic  bacteria  (gSRP  m"3) 

Half-saturation  constant  for 

P-effect    (gSRP  m"3) 


BACTFRAC  =  .5 
AUTOFRAC  =  .42 
EUKFRAC  =  .08 

PBACTCONT  =  .008 
PAUTOCONT  =  .004 
PEUKCONT  =  .004 

PDETRITUSCONT  =  .002 
PINORGCONT  =  .004 

BENTHICBACTFRAC  =  .75 


Bacterial  fraction  of  plankton 

(Ssact  ^pian^ton   ) 

Algal  fraction  of  plankton 

(^Algae   ^plankton      ) 
Eukaryotic  heterotroph 

fraction  of  plankton 

(Sgukhet  ^plankton   ) 

P-content  of  bacterio-plankton 

(9p   ^Bacteria      ) 
P-content  of  algae 

(gP    ^Algae"1) 

P-content  of  eukaryotic 

heterotroph-plankton 

(9p    ^Eukhet      ) 
P-content  of  suspended 

detritus    (gp  goetritus"1) 
P-content  of   inorganic  seston 

(9p   9inorgseston      ) 

Bacterial   fraction  of  benthos 
(^Bact    ^Benthos      ) 
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BENTHICAUTOFRAC  =  .2 
BENTHICEUKFRAC  =  .05 

PBENTHICBACTCONT  =  .012 
PBENTHICAUTOCONT  =  .004 
PBENTHICEUKCONT  =  .004 

PDETSEDCONT  =  .002 
PINORGSEDCONT  =  .003 


Algal  fraction  of  benthos 

(^Algae  ^Benthos   ) 
Eukaryotic  heterotroph 

fraction  of  benthos 

(^Eukhet  ^Benthos   ) 

P-content  of  bacterial-benthos 

(9p   ^Bacteria      ) 

P-content  of  benthic  algae 

(9p   ^Algae      ) 

P-content  of  eukaryotic 

heterotroph-benthos 
(9p    ^Eukhet      ) 

P-content  of  detrital  surfsed 

(9p   ^Detsed      ) 

P-content  of   inorganic  surfsed 

(9p   ^inorgsed      ' 


KF   =    .045 

PLOAD   =    .05 

Q   =    1350 

KT   =    .02 

LITRESTM  =    .15 

M  =    .003 

B2    =    .007 

SUBMLITAREA   =    6000000 

LAKEVOLUME   =    40500000 

SLSR   =    .005 

HALFLIFE   =    5 

B3    =    64.42 

N   =    -5.07 


Model  coefficient  (for  filter) 

Srp-contribution  from  runoff 

(9§rp  m'V1) 

Gain  in  Litpart  from  influx 
and  production  (g  m"2y_1) 
Model  coefficient  (for 
trapping  of  Litpart) 
Litpart  residence  time  (y) 

Slope  coefficient  for  specific 

export  rate 

Y-intercept  for  specific 

(m2) 


export  rate  (g  g_-Ly"x) 
Submerged  littoral  area 


r-ll 


Lake  volume  (m3) 

Specific  littoral 
sedimentation  rate  (%  y"x) 
Littoral  area  residence  time 

(y) 

Model  coefficient  for  gain  in 

veglitarea 

Model  coefficient  for  gain  in 

veglitarea 


PLITDETRITUSCONT   =    .002 
PLITINORGCONT   =    .004 
PLITPLANKTONCONT   =    .006 


P-content  of  detrital  Litpart 

(9p   ^Litdetritus      ) 

P-content  of   inorganic  Litpart 

(9p    "?Litinorg      ) 
P-content  of  planktonic 

Litpart    (gp  gLitpiankton-1) 
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DETRITUSLITFRAC  =  .2 
INORGLITFRAC  =  .2 
PLANKTONLITFRAC  =  .6 


Detrital  fraction  of  Litpart 

(^Detritys    9i,itpai;t      ) 
Inorganic  fraction  of  Litpart 

(Sinorg   litpart      ) 

Planktonic  fraction  of  Litpart 


-li 


(^plankton  9l it part   ) 


PDELLITDETRITUSCONT  =  .002 
PDELLITINORGCONT  =  .004 
PDELLITPLANKTONCONT  =  .006 
DETRITUSDELLITFRAC  =  .2 
INORGDELLITFRAC  =  .2 
PLANKTONDELLITFRAC  =  .6 


P-content  of  detrital 

Dellitpart  (gp  goeii^detr"1) 
P-content  of  inorganic 
Dellitpart  (gp  g^mti^g-1) 
P-content  of  planktonic 

Dellitpart  (gp  goeiiitpiank"1) 
Detrital  fraction  of 
Dellitpart  (gUetj.  gJ)BllLtpart-^ 
Inorganic  fraction  of 

Dellitpart  (gInorg  gDeiiitpart"1) 
Planktonic  fraction  of 
Dellitpart  (gPlank  gDeiiitpart"1) 


VI  =  500 


KR1  =2.67 


KR2  =3.06 


V2  =  500 


Specific  reaction  rate  for 
Pinorg  -*  Srpwater 

(^Srowater,  ^Pinorg  Y     ) 
Eguiiibnum  constant  for 

Pinorg  ■*  Srpwater 

'^p^norg  ^Srpwater   » 
Equilibrium  constant  for 

Pinorgsed  -»  Srppore 

'^Pinprgsed   ^Srppore      ' 
Specific  reaction  rate  for 

Pinorgsed  -»  Srppore 
^Srppore   ^Pinorgsed     Y      ) 


PLANKSPECREGRT  =500 

BENTHOS SPECREGRT  =500 

PLIFETM  =  .05 
BLIFETM  =  .025 
KDETRITUS  =20 

KDETSED  =20 


Specific  regeneration  rate  of 
Srp  by  euk. he t. plankton 

(9srp  9pukhet  Y     ) 

Specific  regeneration  rate  of 

Srp  by  euk . het . benthos 

(9srp  fyukhat  Y     ) 

Mean  lifetime  of  Pplankton  (y) 


Mean  lifetime  of  Pbenthos  (y) 

Specific  rate  of  plankton 
decomposition 

(^Pplankton  ^Pdetritus  Y     ) 
Specific  rate  of  benthos 

decomposition 
(^Pbenthos  Spdetsed  Y     ) 
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COMPUTED  CONSTANTS 

AMP  =  AMPNODAM 

HYDRESTM  =  HYDRESTMNODAM 

BWD  =  BWDNODAM 

KIM  =  SUBMLITAREA  /  LAKEVOLUME 

Conversion  factor  from  m2 
littoral  to  m3  profunda 1  zone 

PPLANKTONCONT  =  (PBACTCONT  *  BACTFRAC)  +  (PAUTOCONT  * 
AUTOFRAC)  +  (PEUKCONT  *  EUKFRAC) 

P-content  of  sestonic  plankton 

PAUTOFRAC  =  (PAUTOCONT  *  AUTOFRAC ) ^PPLANKTONCONT 

Algal-P  fraction  of  Pplankton 
PBACTFRAC  =  (PBACTCONT  *  BACTFRAC)  /  PPLANKTONCONT 

Bacterial-P  fraction  of 

Pplankton 
PEUKFRAC  =  (PEUKCONT  *  EUKFRAC)  /  PPLANKTONCONT 

Euk.Het.-P  fraction  of 

Pplankton 
PBENTHOSCONT  =  (PBENTHICBACTCONT  *  BENTHICBACTFRAC)  + 

(PBENTHICAUTOCONT  *  BENTHICAUTOFRAC)  +  (PBENTHICEUKCONT 
*  BENTHICEUKFRAC) 

P-content  of  benthos 

PAUTOBENTHICFRAC  =  (PBENTHICAUTOCONTen*hBENTHICAUTOFRACl  / 
PBENTHOSCONT  ' 

Algal-P  fraction  of  Pbenthos 
PBACTBENTHICFRAC  =  (PBENTHICBACTCONT  *  BENTHICBACTFRAC )     I 
PBENTHOSCONT 

Bacterial-P  fraction  of 
Pbenthos 

PEUKBENTHICFRAC  =  (PBENTHICEUKCONT  *  BENTHICEUKFRAC)  / 
PBENTHOSCONT 

Euk.Het.-P  fraction  of 
Pbenthos 
Al  =  l  /  (-OPTPAR2)  Model  coefficient  for  Pareff 

Bl  =  -Al  *  MAXPAR  Model  coefficient  for  Pareff 

GOSUB  5000  Call  the  plot-initialization 

subroutine 

FOR  T  =  TO  TO  TEND  STEP  DT 


188 
AUXILIARY  EQUATIONS 

IF  T  >=  15  THEN 

AMP  =  AMPDAM 
HYDRESTM  =  HYDRESTMDAM 
BWD  =  BWDDAM 
RISE  =  RISEDAM 

END  IF 

IF  T  >=  15  +  RAMPDURATION  THEN  RISE  =  0 

LEVEL  =  AVGLEVEL  +  AMP  *  SIN((T  /  7.957)  *  (FREQLEVEL) ) 

Annual  lake  level  fluctuation 
(m  msl) 
BOTTOM  =  (PERMSED  /  BEDDENS)  +  (SURFSED  /  BULKDENS)  + 
LAKEBED 

Mean  elevation  (m  msl)  of 

„„„„,„  sediment-water  interface 

DEPTH  =  LEVEL  -  BOTTOM 

Water  column  depth  (m) 

SRR  =  ALPHA  *  EXP (GAMMA  *  DEPTH) 

Specific  rate  of  resuspension 
by  wind-generated  waves  and 
currents  (g  g_1v_1) 

SURFSED  =  PINORGSED  /  PINORGSEDCONT  +  PBENTHOS  / 
PBENTHOSCONT  +  PDETSED  /  PDETSEDCONT 

Flocculent  Surfsed  subject  to 
resuspension  (a   m~2) 

VSRPPORE  =  SRPPORE  *  BULKDENS  /  SURFSED 

Srppore  (g-^  m"3) 
SPECBENTHICBACTPUR  =  (MAXBENTHICBACT  *  VSF&PORE)  /  (KM1  + 
VSRPPORE)  '  '     K 

Spec .  dark  P-uptake  rate  for 
benthic  bacteria 

*  —1—1 

BENTHICBACTPUR  =  PBENTHOS  *  SPECBE^raiCBACTPUR  * Y  ) 
PBACTBENTHICFRAC 

P-uptake  rate  by  benthic 
bacteria  (ge   m~2v-1} 
SPECBENTHICAUTOPUR  =  (MAXBENTHICAUTO  *  VSRPPORE)  /  (KM2  + 
VSRPPORE)  '  '  l 

Spec,  dark  P-uptake  rate  for 
benthic  algae 

BENTHICAUTOPUR  =  PBENTHOS  *  SPECBENraiCAOTOPUR  * ^  ^ 
PAUTOBENTHICFRAC 

P-uptake  rate  by  benthic  algae 
(gSrp  m^y"1) 
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SESTON  =  PINORG  /  PINORGCONT  +  PPLANKTON  /  PPLANKTONCONT  + 
PDETRITUS  /  PDETRITUSCONT 

Inorganic,  dead,  and  live 
seston  combined  (gSeBton  m~3) 

PAR  =  MAXPAR  *  EXP (-BETA  *  SESTON) 

Photosynthetically  active 
radiation  (cal  cm^hr-1) 

PAREFF  =  Al  *  (PAR2)  +  Bl  *  PAR 

Effect  of  Par  on 
Specplankautopur 

PEFF  =  SRPWATER  /  KM4  +  SRPWATER 

Effect  of  P  on 
Specplankautopur  (0<Peff<l) 


SPECPLANKAUTOPUR  =  MAXPLANKAUTO  *  SQR (PAREFF  *  PEFF) 

Spec.  P-uptake  rate  for 
planktonic  algae 

(9srp  9pPiankauto  Y   ) 
PLANKAUTOPUR  =  PAUTOFRAC  *  PPLANKTON  *  SPECPLANKAUTOPUR 

P-uptake  rate  by  planktonic 
algae  (gSRP  m^y-1) 
SPECPLANKBACTPUR  =  (MAXPLANKBACT  *  SRPWATER)  /  (KM3  + 
SRPWATER) 

Spec.  P-uptake  rate  for 
planktonic  bacteria 

(9srp  9pPiankbact  Y     ) 
PLANKBACTPUR  =  PBACTFRAC  *  PPLANKTON  *  SPECPLANKBACTPUR 

P-uptake  rate  by  planktonic 

bacteria  (gSRP  nT^y-1) 

FILTER  =  KF  *  VEGLITAREA      Percent  P  filtered  from  runoff 

by  littoral 
TRAPEFF  =  KT  *  VEGLITAREA     Percent  of  part,  matter 

trapped  by  littoral 
SER  =  M  *  SRR  +  B2  Spec,  export  rate  (g  g"1y"1) 

RATE-EQUATIONS 

PPLANKTONFLUSH  =  PPLANKTON  *  BWD  /  HYDRESTM 

Flushing  of  Pplankton  through 

outlet  (gpplankton  nT3y-1) 
PINORGFLUSH  =  PINORG  *  BWD  /  HYDRESTM 

Flushing  of  Pinorg  through 

outlet  (gPinorg  nrV1) 

PDETRITUSFLUSH  =  PDETRITUS  *  BWD  /  HYDRESTM 

Flushing  of  Pdetritus  through 
outlet  (gPdetritus  hTV1) 

SRPFLUSH  =  (1  /  HYDRESTM)  *  SRPWATER 

Flushing  of  Srpwater  through 
outlet  (gSrp  hTV1) 
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SEDRT  =  SDR  *  SURFSED 

Sedimentation  rate 

(9surfs*»d  m~2y_1) 
PBENTHOSSEDRT  =  SDR  *  PBENTHOS 

Pbenthos-sedimentation  rate 
a  in™  v™  \ 

PDETSEDRT  =  SDR  *  PDETSED       Pbenthos 

Pdetsed-sedimentation  rate 

(Spdetsed  m~2y_1) 
PINORGSEDSEDRT  =  SDR  *  PINORGSED  ae"ea 

Pinorgsed-sedimentation  rate 

(Spinorgsed  ^Y'1) 

APBENTHOSRESUS  =  SRR  *  PBENTHOS 

Areal  rate  of  Pbenthos 
resuspension  (gDh«„*-v„,=  m~2y_1) 

VPBENTHOSRESUS  =  SRR  *  PBENTHOS  /  DEPTH       "»ntnos 

Volumetric  rate  of  Pbenthos 
resuspension  (gPbenthoa  m_3y-1) 
APDETSEDRESUS  =  SRR  *  PDETSED  roentnos 

Areal  rate  of  Pdetsed 

resuspension  (gPdetSed  m"2y_1) 
VPDETSEDRESUS  =  SRR  *  PDETSED  /  DEPTH         *aetsea 

Volumetric  rate  of  Pdetsed 

resuspension  (gPdetsed  m~3y_1) 
APINORGSEDRESUS  =  SRR  *  PINORGSED  raersea 

Areal  rate  of  Pinorgsed 
resuspension  (gPinoraaed  m"2y_1) 
VPINORGSEDRESUS  =  SRR  *  PINORGSED  /  DEPTH     Pinorgsed 

Volumetric  rate  of  Pinorgsed 

resuspension  (gPinoraaed  m"3y-1) 
AENTR  =  SRR  *  SRPPORE  Pinorgsed    i      J 

Areal  entrainment  of  Srppore 
by  shear  (gSrMOPe  m"2y_1) 
VENTR  =  SRR  *  SRPPORE  /  DEPTH  "ppore 

Vol. entrainment  of  Srppore  by 
shear  (gSrppore  m'V1) 

APPLANKTONSETTL  =  SVELP  *  PPLANKTON 

Areal  rate  of  settling  of 
Pplankton  (gPr,ia „!,<-«„  m~2y-1) 

VPPLANKTONSETTL  =  SVELP  *  PPLANKTON  /  DEPTIT 

Volumetric  rate  of  settling  of 
Pplankton  (gPn,. _,,«.,„,  m"3y-1) 
APDETRITUSSETTL  =  SVELD  *  PDETRITUS        fpianjcton 

Areal  rate  of  settling  of 
Pdetritus  (gPdetritus  m-2y-1) 
VPDETRITUSSETTL  =  SVELD  *  PDETRITUS  /  DEPTH 

Volumetric  rate  of  settling  of 
Pdetritus  (gPdetritu8  m-V1) 
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APINORGSETTL  =  SVELI  *  PINORG 

Areal  rate  of  settling  of 
Pinorg  (gPinorg  m~2y-1) 

VPINORGSETTL  =  SVELI  *  PINORG  /  DEPTH 

Volumetric  rate  of  settling  of 
Pinorg  (gPinorg  m_3y-1) 

PINORGSORP  =  VI  *  (PINORG  -  KR1  *  SRPWATER) 

Adsorption/desorption  rate 

between  Pinorg  and  Srpwater 

(gSr pv/ater  m~  Y~*) 
PINORGSEDSORP  =  V2  *  (PINORGSED  -  10*2  *  SRPPORE) 

Adsorption/desorption  rate 

between  Pinorgsed  and  Srppore 

(9srppore   m'V1) 

PREG  =   PEUKFRAC    *    PPLANKTON   *   PLANKSPECREGRT 

Srp-regeneration  rate  by  euk. 
het.    plankton    (gSrpwater  nrV1) 

PLANKPUR   =   PLANKBACTPUR  +   PLANKAUTOPUR 

Phosporus  uptake  rate 

(9srP  m"3y_1) 
PPLANKTONDR  =  PPLANKTON  /  PLIFETM 

Average  death  rate  of 
Pplankton  (gPpiankton  m_3y-1) 
PDETRITUSDECRT  =  PDETRITUS  *  KDETRITUS 

Average  detritus  decomposition 

rate  (gPdetritus  m_3y"1) 
PPOREREG  =  PEUKBENTHICFRAC  *  PBENTHOS  *  BENTHOS SPECREGRT 

Srp-regeneration  rate  by  euk. 
het.  benthos  (gSrp  m~2y~*) 

BENTHOSPUR  =  BENTHICBACTPUR  +  BENTHICAUTOPUR 

Phosporus  uptake  rate 
(9srP  m^y-1) 

PBENTHOSDR  =  PBENTHOS  /  BLIFETM 

Average  death  rate  of  Pbenthos 

(Spbenthos  m~2y_1) 
PDETSEDDECRT  =  PDETSED  *  KDETSED 

Average  detsed  decomposition 

rate  (gPdetsed  ^'V1) 

INFLOW  =  Q 

Gain  of  Litpart  from  runoff 

and  prod.  (gLitMrt  m"2y"1) 
TRAP  =  LITPART  *  TRAPEFF  /  LITRESTM 

Accrual  of  Litpart 

(^Litpart  m"2y_1) 
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THROUGH  =  (1  -  TRAPEFF)  *  LITPART  /  LITRESTM 

Particles  not  trapped  by 


LITSEDRT  =  SLSR  *  DELLITPART 
DELTHROUGH  =  DELLITPART  *  SER 


Litarea  (gLitpart  m  V1) 

Littoral  sedimentation  rate 
(9Dellitpart  ^'^'^ 


Delayed  particles  flow  through 
Litarea  (gnoi  m-na^  m"2y_1) 
PEX  =  PLOAD  *  (1  -  FILTER)  /  DEPTH       "eixicparr 

Srp-export  rate  from  littoral 
to  pelagic  (gs   m^y-1) 

GAIN  =  (100  /  (1  +  B3  *  EXP(N  *  AMP))  /  (HALFLIFE  /  .69)) 

Veglitarea  gain  with  lake 
level  fluctuation  (%) 

LOSS  =  VEGLITAREA  /  (HALFLIFE  /  .69) 

Loss  in  Veglitarea  with  lake 
level  fluctuation  (%) 

PPLANKTONIMPORT  =  (LITPART  *  KIM  *  (1  -  TRAPEFF)  /  LITRESTM) 

*  PLANKTONLITFRAC  *  PLITPLANKTONCONT 

Import  of  Pplankton  from 
Litarea  (gp  ■,     tt.   m~^y-^) 
PDETRITUSIMPORT  =  (LITPART  *  KIM  *  (1  -  TRAPEFF)  /  LITRESTM) 

*  DETRITUSLITFRAC  *  PLITDETRITUSCONT 

Import  of  Pdetritus  from 

Litarea  (gPdetritua  m'3Y-1) 
PINORGIMPORT  =  (LITPART  *  KIM  *  (1  -  TRAPEFF)  /  LITRESTM)  * 
INORGLITFRAC  *  PLITINORGCONT 

Import  of  Pinorg  from  Litarea 
?inorg 


(9pinorq  "'"V1) 


PPLANKTONDELIMPORT  =  KIM  *  DELLITPART  *  SER  * 
PLANKTONDELLITFRAC  *  PDELLITPLANKTONCONT 

Delayed  import  of  Pplankton 
from  Litarea  (gp  ■.  ,.   m~^y_1) 
PDETRITUSDELIMPORT  =  KIM  *  DELLITPART  *  SER  *ppianJcton 
DETRITUSDELLITFRAC  *  PDELLITDETRITUSCONT 

Delayed  import  of  Pdetritus 
from  Litarea  (gPdetritus  m_3y_:L) 
PINORGDELIMPORT  =  KIM  *  DELLITPART  *  SER  *  INORGDELLITFRAC  * 
PDELLITINORGCONT 

Delayed  import  of  Pinorg  from 
Litarea  (gPinorg  nrV1) 
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PVAR(l) 
PVAR(2) 
PVAR(3) 
PVAR(4) 
PVAR(5) 
PVAR(6) 


SESTON 

SURFSED 

SRPWATER 

VEGLITAREA 

PPLANKTON 

SEDRT 


1st  plot 
2nd  plot 
3rd  plot 
4th  plot 
5th  plot 
6th  plot 


variable 
variable 
variable 
variable 
variable 
variable 


IF  T  =  0  THEN  GOSUB  5100 

IF  PLOTCOUNT  >=  PLTINTVL  AND  PLOTCOUNT  <  PLTINTVL  +  DT  THEN 


GOSUB  5100 

PLOTCOUNT  =  PLOTCOUNT  +  DT 

LOCATE  1,  1 

PRINT  USING  "####.## 

FOR  1%  =  1  TO  NPVAR 
COLOR  COLR(I%) 
PRINT  USING  "####.## 

NEXT  1% 


Call  the  plotting  subroutine 


II  •   rp  . 


";  PVAR(I%) ; 


IF  T  >  15  AND  T  <  15.001  THEN 
LOCATE  30,  1 

PRINT  USING  ••#####.###   ■•;  T; 
FOR  1%  =  1  TO  NPVAR 

COLOR  COLR(I%) 

PRINT  USING  »#####.##   ••;  PVAR(I%) ; 
NEXT  1% 

END  IF 


ACCUMULATIONS 

VEGLITAREA  =  VEGLITAREA  +  DT  *  (GAIN  -  LOSS) 

LITPART  =  LITPART  +  DT  *  (INFLOW  -  TRAP  -  THROUGH) 

DELLITPART  =  DELLITPART  +  DT  *  (TRAP  -  LITSEDRT  - 
DELTHROUGH) 

SRPWATER  =  SRPWATER  +  DT  *  (PEX  +  VENTR  +  PINORGSORP  +  PREG 
-  PLANKPUR  -  SRPFLUSH) 

PPLANKTON  =  PPLANKTON  +  DT  *  (PLANKPUR  +  VPBENTHOSRESUS  + 

PDETRITUSDECRT  +  PPLANKTONIMPORT  +  PPLANKTONDELIMPORT  - 
VPPLANKTONSETTL  -  PREG  -  PPLANKTONDR  -  PPLANKTONFLUSH) 


194 

PDETRITUS  =  PDETRITUS  +  DT  *  (PPLANKTONDR  +  VPDETSEDRESUS  + 
PDETRITUSIMPORT  +  PDETRITUSDELIMPORT  -  PDETRITUSDECRT  - 
VPDETRITUSSETTL  -  PDETRITUSFLUSH) 

PINORG  =  PINORG  +  DT  *  (VPINORGSEDRESUS  +  PINORGIMPORT  + 
PINORGDELIMPORT  -  VPINORGSETTL  -  PINORGSORP  - 
PINORGFLUSH) 

SRPPORE  =  SRPPORE  +  DT  *  (PINORGSEDSORP  +  PPOREREG  - 
BENTHOSPUR  -  AENTR) 

PBENTHOS  =  PBENTHOS  +  DT  *  (BENTHOSPUR  +  APPLANKTONSETTL  + 
PDETSEDDECRT  -  PBENTHOSSEDRT  -  APBENTHOSRESUS  - 
PPOREREG  -  PBENTHOSDR) 

PDETSED  =  PDETSED  +  DT  *  (APDETRITUSSETTL  +  PBENTHOSDR  - 
PDETSEDDECRT  -  PDETSEDRT  -  APDETSEDRESUS) 

PINORGSED  =  PINORGSED  +  DT  *  (APINORGSETTL  -  PINORGSEDRT  - 
APINORGSEDRESUS  -  PINORGSEDSORP) 

PERMSED  =  PERMSED  +  DT  *  SEDRT 

AVGLEVEL  =  AVGLEVEL  +  DT  *  (RISE  /  RAMPDURATION) 


NEXT  T 


Data  below: 


a)  number  of  plot  variables  (maximum  of  6) 

b)  for  each  plot  variable: 

1)  a  legend  symbol 

2)  a  color  number 

3)  the  minimum  scale  value 

4)  the  maximum  scale  value 


DATA  6 

DATA  "SEST",12,0,60 
DATA  "SURF" ,13,1000,4000 
DATA  "SRPW",10,0.0,0.2 
DATA  "LITA",15,0,20 
DATA  "PPLA",11,0,0.2 
DATA  "SEDR", 14, 0,1000 

SLEEP 


STOP 
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5000  Color  plot  initialization 

CLS 

SCREEN  12 
COLOR  14 
PRINT 

PRINT  SPC(40  -  (LEN(TITLE$)  /  2));  TITLE$ 
PRINT  SPC(40  -  (LEN(SUBTITLE$)  /  2));  SUBTITLE$; 
VIEW  (138,  90)-(500,  405),  ,  14 
WINDOW  (0,  10) -(10,  0) 

FOR  1%  =  1  TO  9 

LINE  (1%,  0)-(I%,  .1),  14 

LINE  (1%,  10)-(I%,  9.9),  14 

LINE  (0,  I%)-(.1,  1%),  14 

LINE  (10,  I%)-(9.9,  1%),  14 
NEXT  1% 

LOCATE  27,  17:  PRINT  TO; 

LOCATE  27,  39:  PRINT  THALF; 

LOCATE  27,  62:  PRINT  TEND; 

LOCATE  29,  38:  PRINT  T$; 

READ  NPVAR 

FOR  1%  =  1  TO  NPVAR 

READ  SYMBL$(I%),  COLR(I%) ,  MIN(I%) ,  MAX(I%) 
COLOR  COLR(I%) 

IF  1%  >  3  THEN  ZZ1  =  43  +  5.2  *  1%  ELSE  ZZ1  =  19  -  5.5  *  1% 

LOCATE  5,  ZZ1:  PRINT  SYMBL$(I%); 
LOCATE  6,  (ZZ1) :  PRINT  MAX(I%)  ; 

LOCATE  16,  (ZZ1) :  PRINT  INT((MAX(I%)  +  MIN(I%) )  /  2) 
LOCATE  26,  (ZZ1) :  PRINT  MIN(I%) ; 
NEXT  1% 

COLOR  10 

LOCATE  16,  2.5:  PRINT  "0.1n 

COLOR  11 

LOCATE  16,  69.5:  PRINT  ".1" 

RETURN 
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5100  Color  plotting  subroutine 

PLOTCOUNT  =  0 

FOR  1%  =  1  TO  NPVAR 

WINDOW  (TO,  MAX(I%))-(TEND,  MIN(I%)) 

PSET  (T,  PVAR(I%)),  COLR(I%) 
NEXT  1% 

IF  T  >  15  AND  T  <  15.001  THEN 
LOCATE  30,  1 
PRINT  USING  ••#####.###   «;  T; 

FOR  1%  =  1  TO  NPVAR 

COLOR  COLR(I%) 

PRINT  USING  »#####.##   »;  PVAR(I%) ; 
NEXT  1% 

END  IF 

RETURN 

END 
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